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Abstract 



The article reviews the physics related to kaon and antikaon production in heavy ion reactions 
, at intermediate energies. Chiral dynamics predicts substantial modifications of the kaon properties 

^> I in a dense nuclear environment. The status of the theoretical predictions as well as experimental 

' evidences for medium effects such as repulsive/attractive mass shifts for /K~ are reviewed. 

. In the vicinity of the thresholds, and even more pronounced below threshold, the production of 

\ strangeness is a highly collective process. Starting from elementary reaction channels the phe- 

nomenology of and K~ production, i.e. freeze-out densities, time scales etc. as derived from 
. experiment and theoretical transport calculations is presented. Below threshold kaon production 

shows a high sensitivity on the nuclear compression reached in heavy ion reactions. This allows to 
put constraints on the nuclear equation-of-state which are finally discussed. 

Q , 
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1 Introduction 



For more than two decades strangeness production has been one of the major research topics of heayy 
physics at intermediate energies. First experiments on strangeness production were already perforr 



in the early SOties at the LBL in Berkeley p]j. These were, however, restricted to collisions of light 
nuclei. In the middle of the QOties the whole field obtained a significant push when the KaoS [2^ and the 
FOPI j3] Collaborations at the Gesellschaft fiir Schwerionenforschung (GSI) in Darmstadt/Germany 
started to deliver high precision data for kaon and antikaon production in heavy systems. The kaon 
spectrometer KaoS has stopped data acquisition in 2004 and will be dismantled at the GSI. Hence it is 
a proper time to draw some resume on what has been achieved during the last twenty years and which 
problems still have to be settled. The present article tries to summarize the status of the field and to 
draw conclusions from the various experimental and theoretical efforts. 

Strangeness production in heavy ion reactions at intermediate energies is of high interest since it 
opens the possibility to attack several fundamental questions of nuclear and hadron physics which 
are not only interesting by themselves but have also astrophysical implications and/or are related to 
fundamental aspects of Quantum Chromo Dynamics (QCD). 

What makes the production of strange hadrons special is the fact that strangeness is exactly con- 
served in hadronic reactions. Open strangeness can be produced by the creation of kaon {K^{us), Js)) 
- antikaon {K~{us), K^{ds)) pairs or by kaon-hyperon pairs. The hyperons carrying one strange quark 
are A{uds) and S {Ti~{dds), Ti^{uds), hyperons. The production of hidden strangeness 

through (j){ss) mesons is possible but suppressed according to the Okubo-Zweig-Iizuka selection rule 
j3]. Since data on (j) production are extremely scarce at intermediate energies, the present article will 
mainly focus on kaon and antikaon production. 

A first consequence of strangeness conservation is the fact that mesons, once produced, cannot 
be absorbed by the surrounding nucleons. This results in a rather long mean free path of about 7 fm 
of mesons in nuclear matter and makes them a suitable 'penetrating' probe for the dense fireball 
produced in heavy ion reactions. Final state interactions such as elastic kaon-nucleon scattering or 
the propagation in potentials influence the dynamics but do not change the total yields. A second 
consequence of strangeness conservation are the high production thresholds. The cheapest way to 
produce an ss pair is the reaction NN — > NAK~^ which has a threshold of £^iab = 1-58 GeV for the 
incident nucleon. When the incident energy per nucleon in a heavy ion reaction is below this value one 
speaks about subthreshold kaon production. Subthreshold kaon production is in particular interesting 
since it ensures that the kaons originate from the high density phase of the reaction. The missing energy 
has to be provided either by the Fermi motion of the nucleons or by energy accumulating multi-step 
reactions. Both processes exclude significant distortions from surface effects if one goes sufficiently far 
below threshold. In combination with the long mean free path subthreshold production is an ideal 
tool to probe compressed nuclear matter in relativistic heavy ion reactions. Indeed, one of the major 
motivations to start the kaon project at the GSI was to explore the nuclear equation-of-state (EOS) 
at supra-normal densities, i.e. significantly above nuclear saturation density po — 0.16 fm^^. A better 
knowledge of the high density behavior of the nuclear EOS is relevant for astrophysical purposes, e.g. 
for the understanding of neutron stars and supernovae explosions. Heavy ion reactions provide the only 
possibility to attack this question experimentally and to constrain theoretical models above saturation 
density. After a more than thirty years quest for the nuclear EOS it seems that by studying subthreshold 
production substantial progress could be achieved in the recent years. 

A second push of more theoretical nature was given to the field when chiral models became more 
and more popular in the late SOties and early 90ties. Chiral symmetry is an exact symmetry of QCD 
in the limit of massless quarks. Since the up and down current quark masses are small, i.e. of the order 
of 5-10 MeV, this symmetry is still approximately fulfilled. In nature it is, however, spontaneously 
broken by the non-vanishing - and large - expectation value of the scalar chiral quark condensate 
{qq) of the QCD vacuum (for a pedagogical review see e.g. [5|)- The spontaneous symmetry breaking, 
similar to the spontaneous magnetization of a ferromagnet which breaks the symmetry of the underlying 
Hamiltonian, implies the existence of massless Goldstone bosons which are the pions. The small pion 
mass of 140 MeV ensures that the concept of chiral symmetry is a fundamental feature of low energy 
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hadron physics. By the extension to the full SU(3) sector the pseudoscalar meson octet of vr, rj, K and 
K mesons plays now the role of the Goldstone bosons. With a strange quark mass of about 150 MeV 
the explicit symmetry breaking is much larger in the strange sector than in the SU(2) sector but the 
general concepts of chiral symmetry and its spontaneous breaking are believed to be still valid. They 
do not only allow to understand the origin of the pseudoscalar meson masses. Chiral perturbation 
theory (ChPT) is also considered as the exact, QCD based, theory of the pion-nucleon interaction at 
low energies p. 

Hence it was a natural step to extend ChPT to the SU(3) flavor sector . However, the situation 
turned out to be more complicated than in pure SU(2). While the ttN and the KN interactions can 
be treated perturbatively, the KN interaction is already around threshold dominated by the presence 
of resonances. This makes a perturbative treatment of the antikaon-nucleon interaction impossible. 
Instead, non-perturbative approaches are required, which can be achieved within chiral dynamics, how- 
ever, by the price of loosing of a well defined expansion scheme. This essential difference between kaons 
and antikaons will be reflected at many points in the following discussions. It is the main reason why 
conclusions on in-medium modifications of kaon properties are much firmer than those for antikaons. 

Such medium modifications are closely connected to the conjecture of a partial restoration of chiral 
symmetry in dense matter. The chiral condensate (gg) is expected to be reduced at finite density 
and/or temperature [SI 13 CD] which should be reflected, e.g., in shifts of the corresponding meson 
masses. Heavy ion reactions open thus the possibility to test fundamental concepts of hadron physics 
and QCD. A strong reduction of the K~ mass, as also supported by the investigation of kaonic atoms 
[TTj . will have severe astrophysical consequences. It can lead to K~ condensation in neutron stars 
[3 1121 CSl mi which, due to additional negative charge, increases the proton fraction and softens thus 
the equation-of-state for the neutron star. The onset of K~ condensation is reached when the electron 
chemical potential starts to exceed that of the kaons /ig > which might happen at around 4-5 

times nuclear saturation density. This lowers the maximal neutron star mass to about 1.5 solar masses 
and the core of a supernova, if heavier than this value, will collapse into a black hole. This, on the 
other hand will lead to a large number of low mass black holes in the universe |l5j (for a recent review 
see dni). 

In the late SOties/ early 90ties first mean field calculations [H [TTj predicted already a moderately 

repulsive K+ potential {V ~ +(20-^30) MeV) and a strongly attractive K~ potential {V (100^200) 

MeV) at nuclear saturation density. Such a value of the potential is in accordance with the K^N 
scattering length and theoretical estimates for the potential remained fairly stable over the years. For 
antikaons the situation is more complex. The resonant character of the K~N interaction makes the 
mean field picture highly questionable and the size of the in-medium K~ potential is even not yet 
completely settled from the theoretical side. 

A major challenge studying strangeness production in heavy ion reactions is to verify (or falsify) the 
conjecture of the existence of these two potentials, to determine their size from experiment and to put 
constraints on theoretical models. However, heavy ion reactions are highly dynamical processes. At 
intermediate energies the phase space distributions of the colliding nuclei are far from global and even 
local equilibrium over most of the duration of the reaction ^H]- Since observables are not snapshots 
but the results of space-time integrals over the entire reaction dynamics it is unavoidable to account 
for this dynamical evolution. This means that the link between experiment and the underlying physics 
has to be provided by dynamical transport models. Semi-classical transport equations of a Boltzmann 
type can be derived from non-equilibrium quantum field theory [20] • Corresponding Boltzmann- 
Uehling-Uhlenbeck (BUU) |2I1 [22] niodels or, alternatively, the Quantum-Molecular-Dynamics (QMD) 
approach [23] are well established transport models and explain successfully a large variety of hadronic 
observables in heavy ion reactions, such as e.g. collective fiow pattern or spectra and abundances of 
newly produced particles [211 [23 • Most of the results presented in this review are based on the QMD 
approach. Details of the Tiibingen QMD model can be found in ^Hl [2Z] • 
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Such hadronic transport models propagate one-body phase space distributions in self-consistent 
potentials and account for elastic and inelastic two-body scattering processes. They contain basic 
quantum features such as final state Pauli-blocking for fermions. As soon as particle production comes 
into play the transport approach becomes a coupled- channel problem. The various hadron species 
are coupled via production and absorption processes and by their mean fields. Transport models 
are generally mean field approaches which rely on the quasi-particle approximation (QPA) for stable 
particles and narrow resonances. mesons can safely be treated within such a framework since 

they retain good quasi-particle properties also in a dense environment. Antikaons, in contrast, seem 
to develop complex structures in their spectral functions which would require to account for off-shell 
effects beyond the quasi-particle approximation. There exist first attempts to go in this direction [2H] 
but most available calculations on K~ production are based on the QPA. 

The review is now organized as follows: Chapter 2 gives an overview of the theoretical predictions for 
the in-medium modifications which the kaons should experience in a dense nuclear environment. Since 
chiral perturbation theory is considered as the most suitable tool to study the interactions of pseu- 
doscalar mesons with nucleons at low energies the chapter starts with a short outline of the derivation 
of effective kaon-nucleon models based on chiral dynamics. Mean field models as well as more elabo- 
rated coupled channel dynamics which are required by the resonant structure of the antikaon-nucleon 
interactions, are briefly discussed. For completeness also the expected scenario in the complementary 
case, i.e. in a high temperature, baryon dilute but pion dominated, environment is mentioned. Chapter 
3 summarizes the phenomenology of strangeness production in heavy ion collisions around threshold 
energies. Starting from elementary processes, i.e. strangeness production and exchange reactions, the 
presently accepted scenarios for kaon and antikaon production are developed as they can be deduced 
from experiments and corresponding dynamical simulations. This concerns the questions of system size 
dependences, freeze-out densities, time scales etc.. An important question is in this context the degree 
of equilibration which the mesons reach or do not reach in the expanding system. This question is 
intimately connected with the possibility to probe the early high density phase of a heavy ion reaction. 
If the kaons would have time to equilibrate they would loose their memory on the early reaction stages. 
Chapter 4 turns then to the treatment of strangeness production within dynamical transport models. 
Here we focus on the production and discuss exemplarily the realization within the Tiibingen QMD 
approach. The treatments within transport models used by other groups are similar in principle, but can 
differ in details. Finally a comparison of the predictions from various transport models for selected pion 
and kaon observables is given. The search for signatures of in-medium mass shifts, or more generally, 
in-medium potentials is discussed in Chapter 5. Mass shifts lead first of all to shifts of the production 
threshold which are reflected in the total particle yields. Since conclusions are only possible relative 
to a free mass scenario they have to be based on transport simulations. The data situation strongly 
supports the in-medium mass scenario, at least concerning the mesons. These observations are 
complemented by the study of dynamical observables such as the collective in-plane and out-of-plane 
pattern. Finally Chapter 6 turns to the original issue which motivated the kaon program at the GSI, 
namely to extract information on the nuclear equation-of-state. The measurement of excitation 
functions down to energies far below threshold can be considered as a breakthrough. Before coming 
to the interpretation of the data a brief summary on the present status of the theoretical prediction 
for the nuclear EOS is given. The dependence of the kaon production on the compression achieved in 
heavy ion reactions puts constraints on the EOS which are discussed in detail. Finally the consistency 
with information from other sources such as nucleon and neutron stars is outlined. 



The review closes with a summary of the major results. 
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2 Kaons in dense matter 



2.1 Chiral SU(3) Lagrangian 

The natural framework to study the interaction between pseudoscalar mesons and baryons at low 
energies is chiral perturbation theory (ChPT). Kaplan and Nelson were the first to apply the chiral La- 
grangian to the properties of kaons in nuclear matter [7] . Later on this framework has been used by many 
other authors |Hlinil2iEniEIlElE3EaESlEniEIlEHl- The corresponding chiral SU(3)LxSU(3)ij 
Lagrangian used by Kaplan and Nelson reads 



£ = + A[TrM,(S-l)+h.c. 



+ iTvBYK, B] + DTrBYl'{A^, B} + FTtB^^^'[A^, B] 
+ aiTiB{iMqi + h.c.)S + a2Tri?fi(^Mg^ + h.c.) 
+ a3[TrMgE + h.c.]Trfifi. 



mB)B 



The degrees of freedom in the Lagrangian 



B 



are the baryon octet B 



■^0 



P 

n 



with a degenerate mass m^, and the pseudoscalar meson octet 



^/2 



/ m I 71-° 



\ 



K- 



K+ 



(3) 



entering into the chiral pseudoscalar meson fields 

S = exp(2i0//^) and ^ 



exp(i0//^ 



(4) 



The pseudoscalar meson decay constants are equal in the SU(3)v limit and given by the weak pion 
decay constant /tt — 93 MeV. The current quark mass matrix which is responsible for explicit chiral 
symmetry breaking is given by 


















(5) 



if one neglects the small difference between the up and down quark masses (m^ ~ nid = niq ~ 5.5 MeV). 
The constants F and D are the SU(3) axial vector couplings with F + D = qa which determine the 
pseudo-vector meson-baryon coupling strengths through corresponding Goldberger-Treiman relations. 

Chiral symmetry of QCD is explicitely broken by the finite, but small quark masses (0). However, 
compared to the chiral symmetry breaking scale ~ 47r/^ ~ 1 GeV the up and down quark masses and 
also the strange quark mass [mg — 150 MeV) are small and the QCD Lagrangian is still approximately 
chirally invariant. The same holds for the mesonic sector of the chiral Lagrangian (P). This allows 
a systematic expansion in powers of hadron momenta and light quark masses over A-^, i.e. chiral 
perturbation theory [H]. In the baryonic sector chiral symmetry is broken due to the baryon mass m^. 
In the Lagrangian is degenerate for the baryon octet and the mass spectrum has to be fixed through 
the expansion coefficients. 
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The mesonic vector and axial vector currents are defined as 



V, = ^{eO,^ + ^d,e) and A, = '-{td,^ - ^d.t), (6) 

respectively. 

The treatment of the full Lagrangian (|T)) is complicated since it leads automatically to a coupled 
channel problem. This approach has been pursued by several authors [301 IHHl OH IMl IMl EH 
However, for many applications, in particular studying kaon properties at the mean field level, an 
effective chiral Lagrangian based on kaon and nucleon degrees of freedom can be used. We discuss the 
effective KN Lagrangian approach first and turn then to the more involved coupled channel problem. 



2.2 Effective chiral Lagrangian 

The Lagrangian (jT} can be reduced to an effective Lagrangian by the steps outlined below [71 [HHl • One 
should, however, be aware that the following expansion is not chiral perturbation theory since it mixes 
contributions of different order in ChPT. 

First the pseudoscalar meson field S is expanded up to order 1//^ and only the kaon field K is kept . 
The terms involving the axial vector current do not contribute to the kaon mass and can be ignored. 
The first two terms in Eq. are the kinetic and the mass term 

d^Kd^K - A{mq + ms)KK , (7) 

from where one can identify the kaon mass 

= A{mq + nis) ■ (8) 



The kaon field is given by 



K = ( ^0 ) and K = {K- K^) . (9) 



Keeping explicitly only nucleon and kaon degrees of freedom, the third and fourth terms in Eq. (Q) lead 
to the Dirac equation for the nucleon field 

N = and iV = (p n), (10) 

and the Weinberg- Tomozawa KN interaction term j4Uj 

3i - - ^ 

N{i-i^d^-mB)N -—N-i^'NK . (11) 

The last three terms in Eq. (Q) can be similarly worked out and lead to a scalar KN interaction, 
the so-called Kaplan-Nelson term [7j 



Tr5(^Mg^ + h.c.)5 = 2mqNN —{mg + m,)KK 



TTBB{^Mq^ + h.c.) = 2m,NN - —^{mg + ms)KK 

J TT 

2NN 

[TrMgE + h.c.]Tr55 = 2{2mg + m,)NN —{niq + nis)KK . (12) 

J n 



Combining these expressions, the full Lagrangian reads 

C = N{iYdf,-mB)N + d''Kd^K - A{mq + ms)KK 

- ^N^NK df,K+ [2m,ai + 2m,a2 + 2(2m, + m,)a3]NN 
NNKK 

- — TTF^ — (m^ + ms)(ai + 2a2 + 403) . (13) 

Now one can fix the remaining free parameters oi, 02, and 03 from the nucleon mass 

niiy = niB — 2[aimq + 02^^ + a^{2mq + m^)] (14) 
and the kaon-nucleon sigma term Skn can be identified 

Skn = -^("^q + "^s)(ai + 2a2 + 403) . (15) 

Eq. (|T^ is obtained using eq. for the nucleon mass and the definition of the sigma term 



1 



1 drri]^ drriN 



-2 dniq dm,^ 



{mq + nis) -— V— (16) 
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The effective chiral kaon-nucleon Lagrangian reads now up to order (1//^) 

J TT 

3i - - ^ 

- —NYNKd.K . (17) 

It contains a vector interaction, the Weinberg- Tomozawa term, which is repulsive for kaons and at- 
tractive for antikaons due to gf-parity. The attractive scalar interaction, the Kaplan-Nelson term, is 
equal for kaons and antikaons. The strength of the Kaplan-Nelson term depends thereby on the magni- 
tude of the kaon-nucleon sigma term Skn- The Weinberg- Tomozawa term is current algebra while the 
Kaplan-Nelson interaction is next to leading order in ChPT. 

In contrast to the pion-nucleon-sigma term which is experimentally well determined from pion- 
nucleon scattering (S^^r — 45 MeV), the kaon-nucleon-sigma term is a relatively uncertain quantity 
since it is related to the strangeness content of the nucleon. In addition to the explicit breaking from 
the quark masses (0) chiral symmetry is spontaneously broken by the large expectation values of the 
scalar quark condensate (qq) ~ ((230 ± 25) MeV)^ of the QCD vacuum. The pion, respectively the 
complete pseudoscalar meson octet assuming SU(3)y symmetry, plays the role of the Goldstone boson 
of chiral symmetry breaking. Thus, like for the pion, the kaon mass mx can be related to the vacuum 
quark condensates {qq) by the Gell-Mann-Oakes-Renner (GOR) relation [H] 

m| = ^(m„ + ms)(MM + ss) . (18) 

m„ in eq. (fTHj) is the up-quark current mass and is the strange quark mass. In the nucleon the right 
hand side of eq. (jl8p defines the kaon-nucleon sigma term 

^KN = ^im^ + ms){N\uu + ss\N) . (19) 

Kaon-nucleon scattering yields values for the isospin averaged sigma term of about Ekat — 400 MeV 
whereas lattice QCD predicts values between 300-450 MeV |H1 113113, heavy baryon ChPT gH Ekn = 
380 ± 40 MeV (1=1) and chiral quark model calculations Skn = 386 MeV. Thus Ekat can range 
from 2m^ up to 450 MeV, however, with the tendency of a value around 400 MeV to establish. 
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2.2.1 Mean field dynamics 

For estimates of kaon mass shifts in nuclear matter and kaon dynamics in heavy ion reactions the above 
Lagrangian (fTTj) is usually applied in mean field approximation. Already in the early QOties mean field 
calculations were carried out in the Nambu-Jona-Lasinio (NJL) model |17j. In the context of chiral 
SU(3) dynamics the mean field approximation means to treat the KN interaction at the tree level. The 
in-medium Klein-Gordon equation for the kaons follows from ()17|) via the Euler-Lagrange equations 



d,d^ ± + {ml - ^p. 



()K±(x) = . (20) 



Here = (N'jf^N) is the nucleon four-vector current and ps = (NN) the scalar baryon density. With 
the vector potential 

V, = (21) 

and an effective kaon mass mj^ defined as 



^K = \K-jfps + V,V'^ (22) 

the Klein-Gordon Eq. ()20|) can be written as 

[id^±tV^f + m*^](l)^±{x) = . (23) 

Thus the vector field is introduced by minimal coupling into the Klein-Gordon with opposite signs for 
and K~ while the effective mass is equal for both. Introducing effective momenta (/c* = {E*, k*)) 
as well 

k; = k^^ (24) 
the Klein-Gordon equation (j20l23p reads in momentum space 

[k*'-m*^]Mk) = . (25) 

Eq. is just the mass-shell constraint for quasi-particles inside the nuclear medium. There exists 
now a complete analogy to the quasi-particle picture for the nucleons in relativistic mean field theory, 
e.g. in the Walecka model of Quantum Hadron Dynamics (QHD) |17j where the nucleon obeys an 
effective Dirac equation 

[It* - m*] u{k) = 

for the in-medium nucleon spinors u. From Eqs. (f^3j) the dispersion relation follows 



E{k) = ko = •^k*^ + m*^ ± Vo . (26) 

In nuclear matter at rest where the space-like components of the vector potential vanish, i.e. V = 
and k* = k, Eq. reduces to 



Eik) = ^lk^ + mi~^p, + V,'±Vo . (27) 

Eq. ()26|) accounts for the full Lorentz structure, a fact which comes into play when heavy ion collisions 
are considered where one has to transform between different reference frames, e.g. the center-of-mass 
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frame of the colliding nuclei and the frame where a kaon is created. Like in electrodynamics the spatial 
components of the vector field give rise to a Lorentz force [IH] as discussed in detail in Chap. 4. 
Now one can also introduce the kaon optical potential through the in-medium dispersion relation 

= k;' - m*^ = kl - mi - 2mKUopt • (28) 

In mean field approximation the difference between the mass shell conditions for kinetic and canonical 
quantities is simply given by the optical or Schroedinger equivalent potential 

f/o tip k) = -S5 ± ^ + ' = - ^^ps (29) 

Here we introduced the total scalar kaon self-energy 



1 

2mK V 



j:s^mx-m*^^-—(:^p^-V^] . (30) 



From Eq. ^ follows 



where |koo| is the asymptotic momentum k^ = k^ — mi of the incoming particle in the Schroedinger 
equation. Thus [/opt of Eq. (jSHI) corresponds to the potential in the non-relativistic Schroedinger 
equation. Again the optical potential is of identical structure as the central part of the optical potential 
for nucleons |17j. However, in the latter case an additional spin-dependent part of the interaction can 
be obtained from a Fouldy-Wouthousen transformation of the Dirac equation. 

In this context it should be mentioned that in the literature it is often not distinguished between 
in-medium 'mass' and in-medium energy shifts, i.e. the energy (j27p at zero momentum is identified as 
an in-medium mass = £^(k = 0). In order to strengthen the analogy with the relativistic mean 
field picture for baryons and to distinguish clearly between the different Lorentz properties - m^ given 
by Eq. (j^^ is by definition a Lorentz scalar - we will distinguish between these two quantities in the 
discussion of medium effects. 

To call mj^ the in-medium mass is also consistent with the general picture of a reduction of meson 
masses in dense and hot hadronic matter. From the Hellmann-Feynman theorem follows a reduction 
of the non-strange quark condensate {uu + dd) which is linear to first order in the nuclear density* p 

ISEj 

{p\uu + dd\ p) ^ ^ _ S^Tv ^ _^ ^^2^ 

{uu + dd) flml 

Model calculations [1711211101, e.g. within the NJL model, predicted a similar reduction of the chiral 
condensate. Assuming an analogous behavior of the strange condensate yields 

M5!i±££ji!>^i_|!2L, + ... . (33) 

{uu + SS) tt'^^K 

Combining now the GOR relation ()18|) with ()33|). the effective kaon mass scales with density as 

< =<-^P + (4) (34) 

which is exactly the form of eq. ()22j) and consistent with the picture of partial restoration of chiral 
symmetry. 

*In isospin symmetric nuclear matter the density is related to the Fermi momentum hy p = 2fc|,/(3 7r^) 



10 



2.2.2 Higher order corrections 

The importance of higher order contributions beyond tree level can be estimated from the empirical 
kaon-nucleon scattering lengths. Such a comparison gives also a feeling how far the mean field model 
complies with kaon-nucleon scattering. 

To lowest order in density the mass or energy shift of a meson is generally given by the forward 
scattering length, in the case of kaons 

AEl{k = 0) = El-ml = -vr f 1 + ^) (a^^p + 3 a^p) p + O (4) . (35) 

The empirical values of the isospin 7 = and 1=1 fC^-nucleon scattering lengths are — 0.02 fm 

and a^^Tv — "0.32 fm which leads to a repulsive mass shift of about 28 MeV at nuclear saturation 
density ikp — 265 MeV). Higher order corrections in the density expansion of (j35|l were found to be 
small for . The k\ correction was found to increases the repulsive i^'"'"-mass shift only by about 
20% [SSj compared to expression This fact makes the density expansion useful in the sector. 
The empirical scattering lengths can now be compared to the tree level Weinberg- Tomozawa interaction 
which yields ci^+^ = fm and — —0.585 fm [SU]. Thus current algebra and the corresponding 

effective KN Lagrangian ()17|) is in rough qualitative agreement with the constraints from low energy 
nucleon scattering. In the effective model the too large vector repulsion 1)211) is compensated by 
the attractive scalar Kaplan- Nelson potential ()34)1 . Compared to the empirical values, the repulsion is 
now, however, overcompensated 

One way to overcome this problem at the mean field level has been suggested by Brown and Rho 
|31j . They subsumed higher correlations of chiral order and above into a medium dependence of 
the pion decay constant which should scale in matter similar as the chiral condensate: Assuming the 
Gell-Mann-Oakes-Renner relation ()18p to be still valid in the medium, one obtains the following relation 
for the in-medium pion decay constant /* 

_ K {P\m\ P) .ggX 

/2 mf {qq) 

From ChPT [50j and vr-mesonic atoms [HI] the (s-wave) pion mass is known to change only slightly with 
density. Using the empirical values of m*(po)/'"^7r ~ 1.05, one obtains 

ff{Po)/f:^0.6 (37) 

at nuclear saturation density po — 0.16 fm~^. A dropping pion decay constant enhances both, the vector 
repulsion and the scalar attraction. However, as found in one-loop ChPT the additional attraction is 
counterbalanced by the range term which is of the same order as the Kaplan-Nelson term. Similar 
results have been found for pions JET . Moreover, such a dropping of the pion decay constant seems to 
be supported by the potentials extracted from pionic atoms |^5^. At the mean field level these results 
can be incorporated by replacing H-i> f*'^ only in the vector potential (PT|) . 

Fig. n] shows the in-medium energy shift of and K~ and the in- medium mass defined by (|22p in 
nuclear matter. MFT denotes thereby the Lagrangian (|T7j) with a value of Skn = 350 MeV which has 
originally been used by Li and Ko MFT+corr. denotes the mean field model proposed by Brown 
et al. including the above mentioned higher order corrections with a value of Skn = 450 MeV. The 
MFT and MFT-(-corr. curves shwon in Fig. Q are obtained by Eqs. ()22I27|) with the corresponding 
values for Skn? however in the MFT+corr. case with the additional replacement of by in the 
vector field ()21|). The empirical energy shifts are shown as well in Fig. ^ For the value is obtained 
by Eq. (jHSj) from the empirical K^N scattering length. The K~ band corresponds to the empirical 
iso-spin averaged K~N scattering length of a^-N = 0.62 ±0.5 fm extracted from kaonic atom data |S3]. 
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The strength of the empirical i^'~-nucleus potential is in the analysis of Batty et al. [S3] quite robust 
against a variation of the nuclear wave functions and the inclusion of p-wave interactions. The s-wave 
potential shown in Fig. Ogives therefore a good impression of the values determined from kaonic atoms. 

Similar results have been obtained with slightly modified versions of the effective chiral Lagrangian 
PHj . in the quark- meson-coupling model [55 shown in Fig. Q and in relativistic mean field calculations 
where kaons are coupled to static (J,uj,p and 6 meson background fields |5(ij . 




Figure 1: In-medium kaon energy (left) and quasi-particle mass (right) in the chiral mean field theory 
(MFT ChPT) and including higher order corrections (MFT ChPT+corr.) Results from the 

mean field quark-meson-coupling (QMC) model [S3] are shown as well. The bands represent the values 
extracted from empirical K^N scattering and K~ atoms [54j. 

In this context it is worth to mention that the philosophy behind the KN mean field model is 
similar to that of effective relativistic nucleon-meson Lagrangiens of QHD |Tf[ IS7] . Both are designed 
to describe in-medium properties, in the latter case nuclear matter and finite nuclei. In both cases the 
structure of the interaction complies with the knowledge from free scattering, however, the models do 
not pretend to give a quantitative description of free scattering data. In the present case this is true 
for the sector. The situation becomes, however, much less satisfying turning to the K~ sector. 

Again one can use the low density theorem to estimate the medium effects to leading order in density. 
The empirical scattering lengths a^I^ — (—1.70 + iO.68) fm and a^l]^ — (0.37 -|- i0.60) fm [13 ISH] 
imply according to ()35|) a repulsive mass shift of 23 MeV and a width of F^^- ~ 147 MeV at saturation 
density. In contrast to K~^, the next order correction to the density expansion of Eq. (|35p is large, 
resulting in a total repulsive mass shift of 55 MeV and a width of F;^- ~ 195 MeV |35J. First of all, 
this questions the convergence of a density expansion for the K~-mode. Moreover, the leading terms 
suggest a repulsive K~ potential which stands in clear contradiction to the empirical knowledge from 
kaonic atom data {cik-n = 0.62 ± 0.5 fm) [S^ suggesting sizable attraction at small density . Finally 
the empirical K~N scattering lengths are in disagreement with the Weinberg- Tomozawa term which 
predicts an attractive mass shift. These facts imply that perturbation theory is not applicable in the 
K~ sector. The reason lies in the existence of a resonance, the A(1405) close the K~p threshold which 
makes the K~p interaction repulsive at threshold. The appearance of resonances requires generally a 
non-perturbative treatment of two-body scattering processes. 
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2.3 Non-perturbative coupled channel dynamics 

A non-perturbative calculation of the full two-body scattering amplitude T requires to solve the 
Lippmann-Schwinger equation, respectively its relativistic counterpart, the Bethe-Salpeter equation. 
For kaon-nucleon scattering the Bethe-Salpeter equation reads schematically 



-KN^KN 



Vkn^kn — i Vkn^mbGbDmTmb^kn 



(3J 



where Gb and Dm are baryon and meson propagators, respectively 

1 „ 1 



Gb{p) 



p - niB + ie 



Duik) 



k"^ — m\j + ie 



(39) 



The Bethe-Salpeter equation ()38p iterates the KN interaction kernel V to infinite order. It is a coupled 
channel equation since it contains not only kaon and nucleon degrees of freedom but involves the 
complete baryon {B = N, A, S) and pseudoscalar meson (M = vr, K) octet © and Q of the chiral 
Lagrangian. The coupling to S's and ?7's can be neglected. The A, E and n degrees of freedom are, 
however, essential for KN scattering. 




2 - 
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Figure 2: Real (left) and imaginary (right) part of the isospin zero s-wave i^'~-nucleon scattering 
amplitude from the iterated Weinberg- Tomozawa interaction in a coupled channel calculation. The 
figure is taken from 



Hence the strategy is different to the vrA^ sector where a perturbative expansion of the SU{2) chiral 
Lagrangian has been demonstrated to be very successful fBH- In the SU{3)v sector the interaction 
kernel V rather than directly the scattering amplitude has to be expanded. This kernel is then iterated 
to all orders in the Bethe-Salpeter equation. The leading order in the expansion of the KN interaction 
is current algebra, i.e. the Weinberg- Tomozawa term (fTT|) which is of chiral order Q. Fig. |21 shows 
the real and imaginary part of the isospin zero s-wave iir~-nucleon scattering amplitude from the 
iterated Weinberg- Tomozawa interaction in a coupled channel calculation from Ref. |3Bi- It is nicely 
demonstrated that using the physical kaon mass the A(1405) is dynamically generated as a pole in the 
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i^'~-proton scattering amplitude. A decrease of the K~ mass leads to a disappearance of the A(1405) 
which will be crucial for the discussion of in-medium effects. 

Expanding V beyond current algebra the corresponding coefficients have to be fixed by KN scat- 
tering data. Coupled channel calculations for s-wave scattering with the interaction kernel truncated 
at chiral order were first carried out by Kaiser et al. jSHI- p-wave contributions have been taken into 
account at the one-loop level by Kolomeitsev et al. inO] and in coupled channel calculations by Lutz 
et al. 0^. Higher partial waves have been taken into account in the G-Matrix calculations of Ramos 
and Oset jSI] and by Tolos et al. jHS EB] , in the latter case with an interaction kernel V based on the 
Jiilich meson-exchange potential |64j . 

After fixing the model parameters from free NK scattering one is now able to systematically incor- 
porate medium effects and to determine thus in-medium scattering amplitudes, mass shifts and spectral 
functions. Medium modifications of the Bethe-Salpeter-equation (jnH|) are the following ones: 

• Pauli-blocking of intermediate nucleon states: The Pauli principle is of course not active for 
hyperons and suppresses NK excitations compared to Ytt excitations. 

• Self-consistency: This means a self-consistent dressing of the K~ propagator 

by the in-medium kaon self-energy H-k-- Since K~ mesons receive a substantial width in the 
medium Hk- is generally complex 

3?nK-(A;) = 2E(k)sRf/opt(^,k) , ^IiK-{k) = ~2Eiy)VK-{E,\^) (41) 

• Dressing of the nucleon propagator: At finite nuclear density the nucleon propagator is 
dressed by the nucleon self-energy S^v due to the interaction with the surrounding nucleons 

Gn{p) = -^ \^ ^. • (42) 

Nucleons are still good quasi-particles and thus Stv = + 7^Sy is real. Scalar and vector 
contributions of S^v can e.g. be taken from the Walecka model of nuclear matter jUj. The same 
holds for the other baryons of the baryon octet (j21) where self-energy contributions can e.g. be 
derived from simple counting of non-strange quarks, e.g. Sa = 2/3S7V. 

• Dressing of the pion propagator: Analogous to the kaons the intermediate pion propagator 
Dj^ 1-^ D* is dressed by a pion self-energy 11^ due to A-hole or A^-hole excitations in the nuclear 
medium. 



The effect of Pauh blocking was first pointed out by Koch and later on studied in detail by Waas, 
Kaiser, Rho and Weise [HHl I34j . Pauli blocking effects were found to play a dominant role since the 
attractive K~N interaction is reduced at finite densities. This acts effectively as a repulsive force which 
shifts the A (1405) resonance above the K~p threshold and leads to a dissolution of this resonance 
at densities above 2 — 3po- Since the existence of the A (1405) was, on the other hand, the origin of 
the repulsive K~N scattering length at threshold, a shift or a dissolution of this resonance causes an 
in-medium K'^ potential which is now close to the tree-level result predicted by the attractive Weinberg- 
Tomozawa term. However, self-consistency, i.e. the dressing of the K~ propagator by the attractive 
potential counteracts the Pauli effect. As pointed out by Lutz [32] a decreasing K~ mass results in a 
negative shift of the A(1405) , regarded as a bound K~p state, and compensates the positive Pauli shift 
to large extent. The position of the A (1405) pole stays now fairly constant but the resonance is still 
substantially broadened and dissolves at high densities (as can be seen from the schematic calculation 
shown in Fig. |21) The influence of dressing of nucleon and hyperon propagator due to short-range A^A^ 
and NY correlations has been investigated in |j66|. 
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2.3.1 In-medium potentials 



The dressing of nucleoli and hyperon propagators is generally included in such type of calculations but 
the effects are of minor importance. A strong influence has, however, a dressing of the pion propagator 
by an attractive pion potential which arises due to nucleon- and A-hole excitations in the medium. As 
shown by Ramos and Oset [HI] and also found by Tolos et al. jH^ the thresholds of the vrA and vrE 
channels are lowered resulting in less attraction for the in-medium potential. 




Figure 3: In-medium K energy obtained in coupled channel calculations which includes Pauli-blo eking 
(Waas et al., [SSI), kaon dressing (Ramos et al., |6lj) and pion dressing (Ramos et al., [HI], Tolos et al. 

Fig. El shows the single particle energy or 'in-medium mass shift' i?(k = 0) = mx + WJopt{E, k = 0) 
for antikaons obtained in various coupled channel calculations. This quantity can be compared to the 
mean field picture although such a comparison has to be taken with care. At finite densities the antikaons 
aquire a substantial in-medium width and do no more behave like good quasi-particles, as assumed in 
a mean field picture. In particular at low momenta the spectral functions can be of complex structure 
without a well defined quasi-particle pole which makes the interpretation of the in-medium self-energy 
Hk in terms of on-shell potentials questionable. However, transport models are usually formulated in 
terms of quasi-particles. Hence, we do not want to refrain from this comparison. The microscopic 
coupled channel calculations deliver an attractive in-medium potential which is significantly smaller 
than in the mean field approaches, in particular when a self-consistent dressing of the kaon propagator 
is taken into account, and even smaller when pion dressing is included. This fact is also reflected in 
the optical potential (real part) shown in Fig. 0] We compare the momentum dependence of 3?f/opt 
at saturation density obtained in various approaches. Results are taken from the chiral mean field 
approach [21], denoted in Fig. [T]as MFT ChPT+corr., the coupled channel calculations of Tolos et al. 
[62j . with and without pion dressing, and a dispersion analysis of K^N and K~N scattering amplitudes 
by Sibirtsev and Gassing |B7|. Since the definition used to extract the optical potential from the self- 
energy H varies in the literature, this comparison is based on relation (PT|) which has been used in [^2] 
instead of Eq. (j2HI)- For the magnitudes of the potential are consistent, i.e. the dispersion analysis 
agrees with the mean field approach at zero momentum. It predicts an almost momentum independent 
potential while [/opt is slightly rising as a function of momentum in mean field models. For K~ all 
models predict a considerably reduced attraction at high kaon momenta, however, the potential depths 
strongly deviate. The self-consistent coupled channel calculations from Schaffner et al. [HH] predict an 
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Figure 4: Optical kaon potentials at nuclear saturation density. Results from the cliiral mean field 
approach (ChPT MFT), from coupled channel calculations (Tolos et al. [H^l) and from the dispersion 
analysis of ^7\) are compared. 



even smaller potential which is of the size of -32 MeV at saturation density. 

The dispersion analysis of Ref. comes close to the mean field result which is, however, not 
astonishing since the authors disregarded the repulsive contributions from the S(1385) and A(1405) 
resonances according to the argument that these resonances should dissolve at finite density. They 
found their parameterization of the K~ potential consistent with data from p + A reactions [HZ] . 

All the microscopic approaches predict K~ potentials of only moderate attraction and are thus in 
stark contrast to the mean field picture and the standard analysis of kaonic atoms JT]. The latter 
suggests a strongly attractive on-shell K~ potential of about 200 MeV at po- K is not yet fully clear if 
the microscopic approaches which comply with kaon-nucleon scattering data, can explain kaonic atoms 
as well. However, there are indications that kaonic atoms explore the antikaon potential only at the 
nuclear surface (HHl HOI and weak K~ potentials describe the available data as well [7T] . A final answer 
would require to account for the full off-shell behavior of self-energy and spectral properties of a bound 
K~ state. However, such calculations have not yet been performed. 

Figure El demonstrates finally the validity of the quasi-particle picture. It shows the kaon and an- 
tikaon spectral functions at saturation density obtained from coupled channel calculations from Ref. 
[7^ (including pion dressing) for different momenta. The kaons have still a clear quasi-particle peak 
which in the medium acquires a finite width. The latter is, however, quite small (less than 5 MeV) 
for small momenta and increases up to 15 MeV at a momentum above 400 MeV which is still moder- 
ate. Hence the quasi-particle picture and the mean field approximation are well justified. As already 
stressed several times, the situation for K~ mesons is quite different. In particular at low momenta 
the spectral functions are broad and of complex structure. At larger momenta a quasi-particle peak 
may still be visible but also here substantial strength is shifted to lower momenta. Thus the mean field 
approximation is questionable for the antikaons. 

2.3.2 In-medium cross sections 

Coupled channel calculations may predict sizeable in-medium modifications of the pion-induced K~ 
production cross sections and the corresponding absorption cross sections nY < — > NK~ . The fact that 
the s-wave A(1405) resonance lies only 27 MeV below the K~p threshold implies a strong coupling to 
this state and requires a non-perturbative treatment. The melting of the A(1405) and S(1385) bound 



16 




Energy [GeV] Energy [G^V] 



Figure 5: In-medium kaon (left) and antikaon (right) spectral functions from coupled channel calcula- 
tions of [72] at saturation density. Results are shown for different kaon momenta: pk = (solid line), 
200 MeV (dotted line), 400 MeV (dashed line), and 600 MeV (dash-dot-dot line). 

states due to Pauli blocking of the intermediate states in the BS-equation ()38|) leads to a dramatic 
increase in particular of the vrS — > NK~ cross section at threshold. In 68^ the enhancement factor 
was found to be more than one order of magnitude at po- However, self-consistency shifts the K~ mass 
below threshold and decreases the available phase space which counteracts the enhancement due to a 
melting A (1405). In the calculations of Schaffner et al. the vrS — ^ NK~ is then only enhanced by 
a factor of two and the vrA — > NK~ is hardly affected at all. In the self-consistent calculations of Lutz 
and Korpa jSZI the predicted in-medium modifications of these cross sections are practically opposite. 
They account in addition for the full in-medium modifications of the K~ spectral distributions and 
obtain a strong enhancement of the vrA — ^ NK~ cross section due to the coupling to the S(1385) 
but almost no changes for the vrS — > NK~ channel. The G-matrix calculations of Tolos et al. [02] 
came to opposite conclusions, namely an almost complete suppression of the pion induced reactions 
in the nuclear environment. Such strong modifications of the K~ production cross sections and the 
corresponding absorption cross sections would have severe consequences for the K~ dynamics in heavy 
ion reactions. 

2.4 Kaons in pion matter 

In heavy ion reactions at intermediate energies the matter is baryon dominated, at ultra-relativistic 
energies at CERN-SPS or at RHIC the matter is, however, pion dominated. For completeness it is 
thus instructive to consider this case as well. The problem of medium modifications experienced by 
kaons in a hot pion gas has been dressed in the early 90ties in Refs. ^] and [Tl] where the authors 
came, however, to opposite conclusions concerning the kaon mass shifts. Recently this problem has 
been picked up by Martemyanov et al. |75j|. In Ref. [73] the kaon self-energy has been determined in 
a model independent way to leading order in pion density, based on ChPT at low temperatures and 
experimental phase shifts at high temperatures. 

Analogous to nuclear matter (j35|) the kaon self-energy lixik'^iE) can be expressed in terms of the 
'kK forward scattering amplitudes for on-shell pions and off-shell kaons. The necessary on-shell ttK 
amplitudes have been evaluated in ChPT to order by several authors (see e.g. [7^ and references 
therein). Near the threshold, the isospin-even (-I-) and odd (— ) tiK scattering amplitudes can be 
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expressed in terms of scattering lengths and effective ranges al^"* and fol^^ 



+ cW(^'-0 (43) 



with p* = p*{y/s, m^r, ttt-k) the cm. momentum of the ttK system The i^'-meson self-energy is then 
obtained by integration over the corresponding pionic Bose distributions 

dP^ = ^J exp(^^^) - l) (44) 



(27r)3 V" T 
with the scalar pion density dpsn = dpT,/{2ET, 



-I{K{k\E) = j A^+\s,0,k^)idp,^+ + dp,^o + dp,^^) + J A^-\s,0,k^){-dp,^+ +dp. 



(45) 



To lowest order ChPT isospin symmetric pion matter does not change the kaon dispersion law. The 
leading order effect appears at the one loop level. The mean field, i.e. the scalar mass shift Suik and the 
vector potential Vk follow from the self-energy at the on-shell point which can be expressed in terms 
of s- and p-wave scattering lengths and s-wave effective ranges: 

]lK{mj^,mK) = -'inp^ ay . (46) 

With (5mK + Vk = n_ft:(m|;, mK)/2mK. From pUj) follows also the vector potential 

+ 2m^mK(6!,+^ + 3aS+^)l . (47) 



Vk - 



+ rriK 

Since ChPT is only valid at temperatures well below the pion mass in Ref. [7^] the high temperature 
behavior has been based on a more phenomenological approach which parameterizes the experimental 
phase shifts and matches smoothly with the one- loop ChPT low-temperature limit. The corresponding 
kaon self-energy at threshold, the mass shift and the vector potential are shown in Fig. El as a function 
of temperature. At T = 170 MeV one obtains a negative mass shift SMk = —33 MeV and a repulsive 
vector potential of Vk = 21 MeV. There exists a remarkable analogy to the nuclear matter case: the 
kaon mass shift at high temperatures is large and negative, the vector potential is large and positive, 
their sum is relatively small and negative. Kaons are therefore bound in pion matter similar to nucleons 
or antikaons in nuclear matter. The vector potential is, however, C-even, distinct from the case of 
nuclear matter. In addition both, kaons and antikaons aquire a substantial in medium width T]^ at 
finite temperature 

In-medium modifications of the kaons have direct implications on the meson properties. At ultra- 
relativistic energies an inconsistency has been observed between the yields measured through the 
dilepton and KK channels. In Pb-|-Pb collisions at CERN/SPS energies (-Eiab = 158 AGeV) the p~^p~ 
yield from decays was measured by the NA50 Collaboration TT and the KK channel by the NA49 
[75] . The leptonic channel was found to be enhanced by a factor of two to four compared to the mesonic 
channel. A similar observation has been made in Au+Au collisions at RHIC {^/sNN = 200 GeV) by 
PHENIX jini IHO] which measured simultaneously the — > e^e~ and KK channels. At RHIC energies 
this apparent enhacement of the leptonic channel can qualitatively be understood by an enhanced in- 
medium width r™*^*^ ~ 2 3r™'^ and the final state interaction of the corresponding kaons inside the 
fireball |73EI1- 
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Figure 6: Medium modifications of kaons in an isospin-symmetric liot pion gas: self-energy 3^11 /^/(2mK), 
mass sliift —SrriK, vector potential Vk, and kaon collision width F^^ versus temperature T. Results are 
taken from [73]. 



2.5 Chiral symmetry restoration? 

At the end of this chapter I want to shortly address the question of chiral symmetry restoration. The 
spontaneously broken chiral symmetry of QCD manifests itself in the large expectation value of the 
scalar quark condensate (qq) of the QCD vacuum while the small but finite current quark masses are 
responsible for the explicit chiral symmetry breaking of QCD. In the chiral limit of vanishing current 
quark masses, the Goldstone bosons of spontaneous chiral symmetry breaking, i.e. the pion, respectively 
the full pseudo-scalar obtect, become massless. This fact is guaranteed by the Gell-Mann-Oakes-Renner 
relation. However, through the GOR relation the pseudoscalar meson masses are also directly propor- 
tional to the scalar condensates. The GOR relation is of leading order in the quark masses. By 
in-medium chiral symmetry restoration one understands now the fact that dropping pseudoscalar me- 
son masses are caused by a reduction of the condensate at finite nuclear density and/or temperature. 
The meson masses approach the chiral limit although the explicit symmetry breaking is still valid. 

How such a scenario is connected with the medium properties obtained from chiral dynamics is, 
however, a non-trivial question. The leading term in the chiral pseudoscalar meson-nucleon interaction 
is the Weinberg- Tomozawa term. At finite nuclear density it is responsible for the splitting of the energy 
levels between the degenerate flavor eigenstates K^{us) and K~{us) due to SU(3) flavor symmetry 
breaking. As indicated in Fig. [71 SU(3) flavor symmetry is broken by the non- vanishing up and down 
quark densities pu/d while the strange quark density ps is still zero. Isospin symmetric pion matter, in 
contrast, is flavor symmetric and there occurs no mass splitting between kaons and antikaons. Charge 
symmetry breaking which occurs in isospin asymmetric nuclear matter leads to an additional mass 
splitting of the different isospin states and K^{ds), respectively K~ and K^{ds). 

In chiral coupled channel dynamics the leading order Weinberg- Tomozawa kaon-nucleon interaction 
is iterated to infinite order. The kaons obtain medium modifications, i.e. mass shifts and changes of 
their spectral distributions. However, in this framework their origin can not easily be traced back to 
a restoration of chiral symmetry. Such an interpretation is to some extent possible when the Kaplan- 
Nelson term in the KN- interaction is taken into account. The Kaplan- Nelson term is of order (qq) and 
therefore directly related to the in-medium condensates via the GOR relation. In the mean field picture 
the two effects, mass splitting (Weinberg- Tomozawa) and dropping mass (Kaplan-Nelson) are clearly 



19 



mass A 
[Mev] ^ 



vacuum 



495 - 



spontanous 
breaking 



explicit 
breaking 



kaons are 
Goldstone 
bosons 



nuclear matter 



symmetric 



Pu=Pd=Ps, p, = 



flavor 
asymmetric 

Pu=Pd'Ps = 




flavor & charge 
asymmetric 

Pd> Pu.Ps = 



K 



current quark masses of QCD 


mu = 





mu=5-10MeV 


md = 





md=5-10MeV 


ms = 





ms = 150 MeV 



Figure 7: Schematic representation of kaon and antikaon energy shifts in dense nuclear matter due to 
SU(3) symmetry breaking (p„ = pd, ps = 0). The Weinberg- Tomozawa term leads to the K~^/K~ mass 
splitting (thin line) while the Kaplan-Nelson term is responsible for the common mass shift (thick line). 



separated and the latter one can be interpreted in terms of a partial chiral symmetry restoration. Fig. 
[7| illustrates this scenario schematically. 

In coupled channel chiral dynamics it is no more straightforward to disentangle ordinary hadronic 
many-body effects from a change of the QCD vacuum. The principle connection of hadronic many-body 
effects with basic QCD quantities such as the chiral condensate is provided by the Hellman-Feynman 
theorem which relates (qq) with the derivative of the QCD-Hamiltonian with respect to the current 
quark masses. In nuclear matter the condensate is in the same way obtained from the total energy 
density £ 

{p\qq\p) = {qq) + . (48) 

Eq. ()48|) implies that many-body correlations, albeit based on hadronic degrees of freedom, provide con- 
tributions to the in-medium quark condensate. The practical use of this type of quark-hadron duality 
is, however, rather limited. While there exist sophisticated models for the treatment of hadronic many- 
body correlations, such as e.g. Brueckner theory, one would have in addition to know the dependence 
of these correlations on the current quark masses. For pionic correlations, lets call them C^^, such a pro- 
cedure may in principle be possible making use of the GOR relation, i.e. dCT^/drriq = dC,, / drnT^rriT, /2mq. 
The derivative of the nucleon mass is given by the pion-nucleon E term dmN/dmq = T^T^N/mq and 
analogous relations exist for the strange quark. However, these derivatives are only valid to leading 
order in ChPT which restricts their applications to the mean field level. Attempts to proceed in such 
a direction have been performed, e.g., by Cohen et al. 0Hj and more recently by Chanfray and Ericson 

IH21. 
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K- 


type 


BB — > BYK+ 


BB — > BBK+K' 


strangeness production 
primary/ secondary 


ttB — > YK+ 


ttB — > BK+K~ 


strangeness production 
secondary 




BY — > BBK- 


strangeness exchange 




TvY < — > BK- 


strangeness exchange 


YK+ — > ttB 




strangeness absorption 


BK+ < — > BK+/^ 


BK- < — > BK-{lt) 


elastic/charge exchange 



Table 1: Elementary hadronic reactions which are relevant for kaon dynamics at intermediate energies. 
B stands for nucleons or nucleon resonances {N, A, N*) and Y for hyperons (A, S). 

3 Phenomenology of strangeness production 
3.1 Elementary reaction channels 

We start with a short discussion of the various reaction channels which determine the kaon and antikaon 
production at intermediate energies, the role these channels play in the complex reaction dynamics and 
their experimental and theoretical knowledge. The various processes are summarized in Tab. ^ They 




B * B' 

Figure 8: Diagrammatic representation of the reaction BB — > BY within the one-boson-exchange 
-|- resonance model. 

can be classified in primary and secondary strangeness production reactions, strangeness and charge 
exchange reactions and strangeness absorption reactions. Strangeness production can be subdivided 
into two classes, namely baryon induced {BB — > BYK~^, BB — > BBK^K~) and pion induced 
{ttB — > YK'^, ttB — > BK^K~) reactions. B stands here for a baryon which can be either a 
nucleon or a nucleon resonance {N, A, N*) and Y for hyperons (A, S). The first class falls into two 
subclasses: primary reactions with two nucleons in the entrance channel and secondary reactions with 
at least one nucleon resonance in the entrance channel. Meson induced reactions are by definition of 
secondary type. Both involve at least a two step process with first the production of the intermediate 
resonance or meson. Primary reactions play the dominant role for kaon and antikaon production at 
higher energies and contribute dominantly to the high momentum part of the spectra. However, at 
subthreshold energies the production rates are dominated by the secondary type reaction mechanism. 
The same feature appears already in proton- nucleus collisions jHSl IHH EH ESI- Other processes like 
the production of multi-strange baryons (S, Q) or hidden strangeness (0) are subdominant and will be 
discussed separately. 
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BB — > BYK+: 

The reaction pp — > pAK~^ sets the threshold for production in free space {Tp = 1.58 GeV). 
Historically this cross section was overestimated for a long time jHEji in particular close to threshold. 
The situation improved when the COSY-11 Collaboration delivered data for pp — > pAK~^ close to 
threshold An open question is still the isospin dependence of this channel. Isotopic relations 

predict {pn — > nAK~^) = 2{pp — > pAK~^) while recent proton-deuteron data indicate that the pn 
cross section might be even larger |HH]- Already rather early it was noticed that secondary reactions 
(A^A — y NYK~^), where the A resonance plays the role of an energy storage, contribute essentially to 
the subthreshold kaon yield in HIC's [HHl IHOl IHH 1021 IHEl IHl] while multi-step processes of higher order 
with two A resonances in the entrance channel are suppressed jHH IHS HH] • Reactions with a A in the 
exit channel are energetically suppressed. Since no experimental data exist for cross sections involving 
nucleon resonances one has here to rely on model calculations. Commonly used parameterizations are 
based on one-boson-exchange models [HZl IHH] or on more microscopic one-boson-exchange + resonance 
models PHI IIUUI IIUH llU2j . Both models are usually applied at the tree level. The resonance model 
assumes thereby that the heavy meson production process runs over the excitation of an intermediate 
nucleon resonance R = N*,A as depicted in Fig. |H1 The model parameters, i.e. coupling constants 
can be fixed from data on vrA^ — > YK^ reactions. As already shown in |^ the resonant s-channel 
is dominant in pp — > pAK~^, non- resonant t-channel contributions can be parameterized through K* 
exchange. This type of resonance models are popular and generally quite successful in the description 
of heavy meson production (p, u, (p) close to threshold |103t I104| I1U5| I106[ I1U7| I108j . 
-kB — > YK+: 

The importance of this channel for HIC's was first pointed out by Fuchs et al. |lU9j and confirmed 
by the Giessen group ^94j. Depending on energy and system size the pion-induced reactions ttN — > 
YK^, IT A — > YK^ can contribute up to about 50% of the total kaon yield. Naturally the contribution 
is largest in heavy systems. The corresponding cross sections are relatively well constrained by pion- 
proton scattering data. The parameterizations of Tsushima et al. |ll(Jj obtained within an one-boson- 
exchange + resonance model are presently used in most transport calculations. 
YK+ — > TiB: 

At subthreshold energies strangeness absorption plays no significant role. That the probability for 
such an event is, however, not complete vanishing has been demonstrated in jlllj . Since strangeness 
absorption would be the only possibility to drive kaons towards chemical equilibrium the chemical kaon 
freeze-out occurs already at a non-equilibrium stage [112] . 
BB — > BBK+K-: 

The reaction pp — > ppK^K^ sets the threshold for production in free space (Tp = 2.50 GeV). 
The history of this reaction is similar to that of pp — ^ pAK^ . Also here early predictions |113j 
overestimated this cross section in the vicinity of the threshold by about two orders of magnitude, as 
it turned out when COSY-11 provided the first data point close to threshold jll4j . Tree level OBE 
calculations based on pion and kaon exchange are able to reproduce the near threshold behavior |115j . 
In contrast to pp — > pAK^ this reaction seems to be dominated by t-channel exchange. Hence the 
resonance model which is based on resonant s-channel exchange has not been applied. Reactions which 
involve nuclear resonances in the initial and/or final state have not been investigated theoretically and 
are neglected in corresponding transport calculations. 
TiB — > BK+K-: 

The vrA^ — > NK^K~~ and ixN — > NttK^K^ cross sections are experimentally well constrained. Cal- 
culations based on tree level one-boson-exchange have been performed in |115j . Transport calculations 
|1161 inni estimate this channel of about equal importance as baryon induced primary K~ production. 
BY i — > BBK-: 

As first pointed out by Ko |117j . strangeness exchange is the dominant mechanism for K~ production 
at subthreshold energies. The relative importance of the two possibilities, through nucleon-hyperon 



22 



scattering or pion absorption are not yet completely settled. More than 20 years ago Ko |117j estimated 
the contribution of this channel to about 10%. The NY — > NNK~ can in principle easily be deter- 
mined in a OBE picture since the cross section factorizes in the pion exchange and the experimentally 
known tiY < — > NK^ part. The importance of this channels is, however, still controversial. In |116j it 
was found to contribute only a few percent to the total antikaon yield in heavy ion reactions whereas 
a recent reanalysis of this channel by Barz et al. |118i stresses its importance (about 25%) for p + A 
reactions. 

TiY < — > BK-: 

This channel is experimentally well known from -proton scattering (for a recent partial wave anal- 
ysis see |119j ). The cross sections are large, e.g. 40-70 mb for NK^ — > ttY depending on relative 
momentum. Strangeness exchange reactions are essential for the K~ dynamics in heavy ion collisions. 
The cross sections are generally large, i.e. of the order of mbarn, whereas the strangeness production is 
of the order of /ibarn. In the standard scenarios strangeness exchange through pion absorption is one of 
the major sources for the K~ yield at subthreshold energies [Ml I120j . The reason lies in the fact that 
the threshold for strangeness production through hyperons in association with mesons is about 
1 GeV below the direct K^K~~ production threshold and these strange quarks can be transfered to 
antikaons by pion absorption. This drives the Y — K~ system towards chemical equilibrium whereas for 
anti-strangeness no comparable process exists. Hence, antikaons freeze out later than kaons |11H 1120"] . 
However, due to the strong coupling of the i^~-nucleon system to the A(1405) large medium effects 
have to be expected for the ttE < — > NK~ cross section and this reaction has theoretically been studied 
in great detail [HHl EZl IHH I121j . The medium dependence of the cross sections is, however, still 
controversial and ranges from moderate enhancement close to threshold to a strong suppression (see 
discussion in chapter 5). This, of course would change the complete scenario of K~ production and 
absorption. 

BK < — > BK, BK < — > BK : 
Because of strangeness conservation the kaon final state interaction is at low energies dominated by 
elastic kaon-nucleon scattering. This does in principle not affect the total yields. However, elastic 
scattering changes the shape of the momentum spectra significantly. Rescattering makes the spectra 
harder and changes the angular distributions |122[ I12c?l I124j . In systems with large isospin asymmetry 
charge exchange pK^ < — > nK^ may also slightly affect the final yields. The elastic antikaon cross 
section is large (~ 20 — 100 mb), in particular at low momenta, which is the result of the strong 

coupling to the A(1405) resonance. It competes, however, with K~ absorption. The corresponding 

_ — 

charge exchange cross section pK~ < — > nK is comparatively small and does not significantly influence 
the K~ dynamics. 

3.2 Hidden strangeness production 

In addition to open strangeness production through KK pairs or kaon-hyperon pairs hidden strangeness 
can be produced through (f){ss) mesons. The is a 0^(1 ) vector meson with the same quantum 
numbers as the uj and has a mass of 1020 MeV. However, compared to the lj the production of (p 
mesons is strongly suppressed due to the Okubo-Zweig-Iizuka (OZI) f3| selection rule which forbids 
the appearance of disconnected quark line diagrams. According to the OZI estimate mesons can 
only be produced due to a small admixture of non-strange light quarks in their wave function. The 
corresponding mixing angle 9mix is equal to 9mix ~ 3.7° |125j . According to this mixing angle the ratio 
of and u mesons cross sections should at comparable energies naively be equal to 

R^/^ = tan^ e^i, ■ F ^ 4.2 X 10"=^ ■ F , (49) 

where F is a correction due to the different phase space factors for and u mesons. Experimentally, 
the ratio R^/u, is in pp — > pp(j){uj) reactions, however, known to be one order of magnitude larger than 
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the naive expectation |126j . Elementary production has been calculated within several approaches, 
mostly within meson exchange models \V27\ 11281 1129j . In Refs. |lU7ij and the cross section for 

LJ and (f) meson production in pp collision have been calculated within the resonance model. The 
experimental data are well reproduced and the large violation of the OZI prediction for the over u 
production, observed experimentally, could be explained in |lU8j without the introduction of additional 
parameters. The reactions of the type NN — > NNcp and vrA^ — > N(f) should be OZI suppressed compared 
to open strangeness production as well, whereas the DISTO results [126 show that almost 50% of the 
pp — > ppK^K~ cross section is resonant, i.e. running over an intermediate (f). When parameterizations 
of the A^A^ — > NNK'^K~ are used in transport calculations the resonant parts of these cross sections 
are usually imphcitly included. 

The FOPI Collaboration measured the production in Ni+Ni reactions at 1.93 AGeV and came to 
the conclusion that at least about 20% of the observed K~ mesons originate from decays |130j . Due 
to the large coupling to the KK channel (the partial r(0 — > K^K^) decay width is about 50%) one 
can expect that the reacts sensitive on medium modifications of the kaons which will mainly result in 
a broadening of the width in hadronic matter (see e.g. [73] and refs. therein). The production in 
heavy ion reactions has been studied in |127t I131j . However, also with the inclusion of medium effects 
the corresponding yields are underestimated by a factor two to three by these transport calculations 
which indicates that production in heavy ion reactions is not yet fully understood. The situation will 
probably be improved when the HADES Collaboration will be able to measure the production through 
the e^e~ decay channel. The dilepton spectra should provide a very clear access to a broadening of the 
and the underlying in-medium modifications of the kaons. 

3.3 Densities and time scales 

Since the conjecture that the early and dense phase of heavy ion reactions can be probed by 
production has been put forward by Aichelin and Ko about two decades ago |132j strong efforts were 
undertaken to settle this question, both from the experimental and the theoretical side. The original 
motivation was to study the nuclear equation of state at supra-normal densities. This point will be 
addressed in detail in Chap. 6. The prerequisite for such an enterprise is of course the fact that 
mesons originate indeed from the early phase and from supra-normal densities. Both features were 
predicted by transport calculations in the early QOties |132t 1^ 1^ 1^ I133j . At that time experimental 
data on subthreshold kaon production were scarce and the same was true for the elementary cross 
sections which serve as input for the transport calculations. Hence these calculations were burdened 
with large uncertainties but the predictions on densities and time scales turned out to be stable. A 
more recent analysis is e.g. given in |lllj . 

3.3.1 Conditions for strangeness production 

mesons are produced in the early and dense phase of a heavy ion reaction. E.g. at an incident 
energy of 1 AGeV they reach their chemical freeze-out values already after 15 fm/c which is the phase 
were maximal compression is reached. One should, however, keep in mind that at such time scales the 
surrounding nuclear environment is still in a pronounced non-equilibrium state |11H ITH] . The same 
is in principle also true for K~ mesons. However, due to high absorption rates through strangeness 
exchange reactions K~N ^ Yn they have the tendency to reach chemical equilibrium much faster. We 
will come back to this point in subsection 3.4. 

Experimental evidence for an early freeze-out is provided by the fact that mesons and high- 
energy pions behave qualitatively similar. This fact is demonstrated by Fig. Elwhere the experimentally 
extracted inverse slope parameters for kaons and high-energy pions {E^^ > 300 MeV) are compared. 
Also rapidity distributions of and high pt pions are similar and clearly distinct from low pt pions 
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Figure 9: Measured inverse slope parameters for high-energy pions and mesons in Au+Au reactions 
at 1 AGeV. The figure is taken from |134j . 

|135j . At SIS energies both, kaons and pions with px > 0.5 GeV are produced at equivalent "subthresh- 
old" energies and identical slopes of the spectra imply that these particles freeze out under similar 
conditions. High-pT particles probe in general the early stage of a heavy ion reactions and this feature 
is independent of the considered energy range. For pions at SIS energies this behavior was predicted 
by transport models |136t I137| I77j and experimentally verified by the simultaneous measurement of 
pion in-plane and out-off- plane fiow |1381 1139j . The energy dependence of the pion fiow is caused by 
spectator matter shadowing which can be used as a time clock for the reaction 1381 1134j . To obtain 
a quantitative picture of the density range probed by production, in Fig. ITUl barvon densities are 
shown at the space-time coordinates where mesons are created. The figure shows the densities in 
central collisions (b=0 fm) for three typical mass systems, Au-f-Au, Ni-|-Ni and C+C at 1.0 AGeV. The 
density distribution dMK+ / dp is defined as 



where is the baryon density at which the kaon i was created and Pi is the corresponding production 
probability. For a better comparison of the different reaction systems the curves are normalized to 
the corresponding mass numbers. Fig. [TUl illustrates the fact that nuclear compression is probed by 
production. Most kaons are produced at supra- normal densities between one and three times po- 
The maximal densities are, on the other hand, correlated with the system size. A dependence on the 
nuclear EOS will therefore be most pronounced in large systems as discussed in detail in chapter 6. 
Closely connected to the densities is the average number of collisions Nc which particles encountered 
prior to the production of a meson. The distributions shown in Fig. ^2 are obtained in analogous 
way as the density distributions, i.e. dM/dNc = dPi/dN}j with Nq = \{N}j^ + N}.^) the average 
number of collisions which two hadrons (1,2 = iV, A, vr) experienced before they produced kaon i with 
production probability Pi. As before the curves are normalized to the mass numbers of the different 
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Figure 10: Nuclear density at the production of a meson. Results from transport calculations for 
central C+C, Ni+Ni and Au+Au collisions at 1 AGeV, normalized to the mass numbers, are shown. 





C+C 


Ni+Ni 


Au+Au 


(1.0 AGeV) 


58 ±6 


75 + 6 


85 + 6 


K' (1.8 AGeV) 


55 ±6 


90 + 5 





Table 2: Experimental inverse slope parameters for and K at comparable energies above threshold 
for different mass systems (from |1H4| I14()| I141j ) . 

reaction systems. This quantity is a suitable measure for the collectivity probed by production. 
Fig. [TT] (left panel) demonstrates that subthreshold kaon production is indeed a multi-step process. At 
1 AGeV the average number of hadron-hadron collisions prior to the production of a kaon is larger than 
one and increases with the system size. This feature is also reflected in the Apart dependence of the 
experimentally extracted inverse slope parameters of the kaon spectra shown in Fig. El In more central 
reactions kaons seem to originate from an environment with apparently higher temperature which must, 
however, not be thermally equilibrated. The inverse slope parameters rise also with system size. Tab. 
13 which is consistent with the collision history shown in Fig. ^2 The inverse slope parameters can 
also simply be considered as a measure of the collectivity experienced by the kaons. The question if 
equilibrium has been reached or not is a subtle one and cannot be decided exclusively from experiment. 
With increasing energy the number of multi-step processes involved in the kaon production decreases 
which can be seen from the right panel of Fig. ^2 There the number of collisions are shown at different 
energies 1 and 1.5 AGeV. For better comparison the curves are normalized to each other. In the light 
C+C system the decrease is more pronounced than in Au+Au and similar to p+A reactions most kaons 
are now produced in primary reactions. 

3.3.2 Energy and mass dependence 

The fact that subthreshold kaon production is a collective phenomenon is clearly reflected in the de- 
pendence of the kaon yield on the number of participating nucleons. Apart is defined as the number of 
nucleons contained in the overlapping region of the two colliding nuclei. For a given impact parameter 
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Figure 11: System size and energy dependence of the number of binary collisions which particles 
encountered prior to the production of a meson. Left: results from transport calculations for 
central Au+Au, Ni+Ni and C+C reactions at 1 AGeV. Right: Same for Au+Au and C+C at 1 and 1.5 
AGeV (arbitrary normalization of the curves). 



the overlap in a geometrical model is that of two hard spheres with radius tqAs [tq = 1.2 fm). In 
symmetric A+A collisions the mean Apart is given by (Apart) = A/2, in exclusive reactions Apart is 
directly related to the centrality. 

However, before turning to heavy ions it is useful to consider proton-nucleus reactions. A survey 
on the experimental situation of the mass dependence of production in p+A can be found in |143j . 
If the kaons were dominantly produced through primary reactions the scaling law of the aP~^^ *k+x 
cross section should follow that of the inelastic proton-nucleus cross section. For proton energies Tp < 2 
GeV the inelastic cross section is given by o"inei ~ j\om±om -^j-^ici^ is close to that of an opaque nucleus 
(cr ~ As). At Tp = 1.5 GeV, i.e. close to threshold, a scaling law aK+ ~ A"" with a = 0.73 ± 0.04 
has been extracted at SATURNE |144j which supports the primary reaction mechanism. Deep below 
threshold the missing mass of the kaons has to be provided by the Fermi motion inside the nucleus or 
the high momentum tails of the nucleon spectral functions. Such collective effects are not supposed to 
depend strongly on the system size. On the other hand, deep subthreshold energies favor the secondary 
production mechanism since, e.g., the Fermi motion can be utilized several times in multi-step processes. 
This would lead to a proportionality ax+ ~ A. Indeed, at PNPI a coefficient a = 1.04 ± 0.01 has been 
extracted at Tp = 0.8 1.0 GeV ^145 j. Such an exponent is in line with the predictions from transport 
models which indicate that subthreshold kaon production in p+A is governed by secondary reactions 
jHni- More recently ANKE extracted values of a = 0.74 ± 0.05 at Tp = 1.0 GeV [113] in fine 
with SATURNE. The scaling depends, however, strongly on the experimantal acceptence and thus 
experimental situation is not fully conclusive. 

In summary, p+A reactions reveal a scaling law of the kaon multiplicity Nk+ = 0"^+ /cinei ~ A" 
with 1.05 < a < 1.48. In a Glauber picture an A+A collision can be understood as the superposition of 
independent NN collisions or as Ax (N+A). Assuming the Glauber picture to hold, the kaon multiplicity 
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Figure 12: multiplicities per Apart for different colliding mass systems at 1 AGeV. Data from KaoS 
|1411 1142j are compared to theoretical QMD transport calculations. For C+C and Ni+Ni also K~ data 
|14H 1140] at equivalent energy above threshold (1.8 AGeV) are shown. 



Nk+ should be linear proportional to Apart- This is evidently not the case. Fig. I12l ' shows KaoS data 
|1341 1142j for different mass systems at 1 AGeV which can be fitted by Nk+ oc A" with a = 1.36. 
The observed non-linear increase as a function of system size is thus a clear indication for a collective 
production mechanism through multi-step scattering processes. The system size dependence of the 
measured multiplicities can be reproduced by theoretical transport calculations. In Fig. El we show 
results from QMD calculations which include a repulsive in-medium potential (see also chapters 4 and 
5). The same non-linear scaling with Apart occurs as a function of centrality, i.e. the yield in 
central reactions is strongly enhanced relative to peripheral reactions. The scaling of Nx+ oc A^.^^^ with 
a = 1.8 ± 0.15 is even more pronounced than the average A dependence of the kaon multiplicity of 
different A-l-A systems (Fig. IT^ . Since the bulk of kaons originates from semi-central collisions (e.g. 
in Au+Au from 6 ~ 5 fm which corresponds to Apart/Amax ~ 0.7 in Fig. IT^ the pronounced centrality 
dependence is to some extent washed out in minimal bias reactions. As a function of centrality a is 
significantly larger than the value of a = 1.48 which can maximally be deduced from p+A within a 
Glauber picture. This demonstrates that heavy ion collisions show qualitatively novel aspects compared 
to p+A. Fig. EJcompares the experimental Apart dependence in Ni+Ni and Au+Au reactions at 1 AGeV 
|14H I134j to QMD calculations. In both cases the data are well reproduced. The contributions of the 
different production mechanisms can be seen from Fig. As expected, the primary NN — > NY 
channels show a linear Apart dependence while the strong non-linear increase of the total yield is caused 
by secondary type reaction mechanisms. Here the processes AA — > NYK^ and ttN — > YK^ 
dominate while the higher order multi-step processes AA — >■ NYK^ and ttA — > YK^ are of minor 
importance. 

For Ni+Ni the scaling parameter a has also been measured as a function of the beam energy 
[TiT] . A slight decrease of a = 1.9 ± 0.25 at 0.8 AGeV to a = 1.65 ± 0.05 at 1.8 AGeV indicates 
the decreasing importance of collective effects with rising beam energy. The experimental facts are 
consistent with the theoretical findings discussed above, namely that the average number of collisions 

^Although in minimal bias reactions the mean ^part is given by A/2, we scale the yields by (Apart) = A/2 which 
accounts for the fact that kaons are dominantly produced in central and semi-central reactions. 
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Figure 13: multiplicities as a function of Apart in Au+Au and Ni+Ni reactions at 1 AGeV. Data 
from KaoS |134t I141j are compared to QMD transport calculations. 



which the nucleons encountered prior to the production of a meson increases a) with system size 
and b) going deeper below threshold. The particles have to accumulate the necessary energy by multiple 
scattering which is easier to achieve in central reactions and/or larger systems where higher densities 
are reached. The conditions are optimal in the early phase of the reaction where baryon as well as 
energy density are maximal. Corresponding data for antikaons are rare. However, existing K~ data 
indicate that the behavior of antikaons follows closely that of the kaons. In Ni+Ni at 1.8 AGeV a scaling 
law of oc Apart "with a = 1.8 ± 0.3 has been observed |14m I141j . At a first glance this may be 
astonishing since one could expect that the strong K~ absorption leads to a weaker Apart dependence 
than for K^. On the other hand, strangeness exchange is the dominant source for K~ production at 
subthreshold energies and primary hyperons are produced in association with the i^T+'s. Since these 
primary hyperons have the same Apart as the K~^^s the K~^s are forced to follow them. 

3.3.3 Chemical freeze-out 

The question of strangeness equilibration has a long history. Already in the classic paper by Koch, 
Miiller and Rafelski |147j it was argued that at CERN energies strangeness saturation requires equi- 
libration times Teq ~ 80 fm/c much longer than the reaction time and that the fact that particle 
spectra can nonetheless be described by a common Boltzmann factor T ~ 170 MeV hints towards a 
sudden hadronization where the system crosses the phase boundary from the quark-gluon plasma to the 
hadronic phase. At SIS energies there is of course no such phase transition but the question remains 
whether the observed particle abundances have purely dynamical origin or can be understood by means 
of statistical concepts. 

Transport calculations indicate that kaons as well as antikaons are dominantly produced in the early 
and dense non-equilibrium phase. Since kaons interact with the surrounding medium almost exclusively 
by elastic reactions KN — > NK (here we consider charge exchange as elastic) they have a long mean 
free path ~ 7 fm and chemical freeze-out occurs early. This picture is supported by the fact that 
subthreshold kaons show similar features like high energy pions, see Fig. El High energy pions were 
experimentally proven to originate from the early phase |135| 11381 1139j . The influence of strangeness 
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Figure 14: multiplicities as a function of Apart in Au+Au at 1 AGeV for the different production 
channels NN (primary), A^A + vrA^ (secondary) and A A + it A (secondary). 



absorption which might drive kaons to equilibrium has been studied in |lllj . This process has extremely 
low probability since at subthreshold energies no more than one YK^ pair per nucleus-nucleus collision 
is created. This means that the same pair of particles has to meet again in the course of the reaction 
for the kaon to get annihilated. As long as kaon production is active the annihilation (loss) rate dNi/dt 
is about one to two orders of magnitude smaller and can be safely neglected in corresponding 

transport calculations. Hence the kaon yield freezes out early and at supra-normal nuclear density 
The conditions are not much different from those where the kaons are produced, see Fig. E3 
The same conclusion can be drawn from the investigations of Bratkovskaya et al. jll2j where chemical 
equilibration times in "infinite" hadron matter, i.e. a box with periodic boundary conditions, were 
studied using a BUU cascade model. The resulting equilibration times for pions and kaons as functions 
of the baryon density and the energy density e are shown in Fig. ^1 In contrast to pions, kaons need 
at least Teq > 40 fm/c to reach chemical equilibrium. Equilibration takes even longer at lower baryon 
and/or energy density. At SIS conditions, corresponding to the lower limit of e in Fig. [TUl r^q is above 
100 fm/c which is at SIS energies about twice the duration of an heavy ion reaction like Au+Au. 

For antikaons the situation is more complex. Though primary antikaon production takes also 
place in the initial and dense phase, strangeness can be redistributed between hyperons and antikaons. 
Strangeness exchange leads to annihilation rates which exceed the production rates in the expansion 
phase of a heavy ion collision |11HI120] . When the fireball becomes dilute the endothermic K~ produc- 
tion process dies out while the exothermic inverse reaction can still continue. Due to such a cross over 
there exists a situation where the equilibrium condition, namely production rate = annihilation rate, is 
formally fulfilled. However, this is not a stationary state since the system is expanding. The situation is 
illustrated by Fig. [TBI where kaon and antikaon production (gain) rates dNc/dt and annihilation (loss) 
rates dN^/dt in a central (b=0 fm) Ni+Ni reaction at 1.0 AGeV are shown. The results were obtained 
by Pal et al. |lllj . Though closer, also K~^s seem not to reach chemical equilibrium in a thermody- 
namical sense. In thermal equilibrium the yields do exclusively depend on the chemical potentials and 
the temperature. As has been demonstrated by Hartnack at al. [120j this is not the case even when 
the equilibrium condition dNc/dt = dNi/dt is fulfilled: In the expanding system an artificial scaling of 
the strangeness exchange cross section ttY NK~ is directly reflected in the total yield. 
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Figure 15: Chemical equilibration time r^q for pions and kaons in infinite hadronic matter as a function 
of baryon density and energy density e. Results from a BUU cascade are shown. The figure is taken 
from [TT^ . 



In view of these facts it is astonishing that many features can also be understood by means of purely 
statistical concepts. The statistical model of an ideal hadron gas assumes that particles are emitted 
from a thermalized, eventually radially expanding fireball. It turned out to describe hadron yields and 
ratios successfully from AGS to RHIC energies (for a recent review see jl48j ). 

In the framework of a grand canonical ensemble the particle number density of species i is given 
by a Boltzmann distribution ^ 



V 



(2^ 



-{Ei-tlBBi-flsSr)/T 



(51) 



Ei = ^/k^ + mf is the single particle energy of hadron species i, T is the temperature, /i^ and jj^s are 
the baryon chemical potential and strangeness chemical potential. Bi and Si are the baryon charge 
and strange charge and Qi the degeneracy factor of hadron i. V is the volume of the system. The 
characteristics of the fireball can be described by means of four independent parameters, such as volume 
V, total energy E, net baryon density pb and net strangeness density ps- These determine the fireball 
temperature T, the strangeness and baryon chemical potentials ps and pb- Since the total strangeness 
is zero the number of free parameters reduces to three, namely, V, E and ps = Nb/V. If one is only 
interested in particle ratios instead of total yields the volume dependence drops out and the number 
of parameters is reduced to two. They are usually chosen as the temperature and the baryon chemical 
potential and are determined by fitting experimental hadron yields. 

Such a grand canonical description fails at SIS energies. However, as shown by Cleymans et al. |149j 
the statistical model works also at SIS if one accounts for exact strangeness conservation which leads 
to the following modification of Eq. ()51|) for kaons and antikaons 
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*In most actual calculations the more accurate Fermi and Bose distributions are used. The usage of Boltzmann 
distributions makes the present discussion, however, more transparent. 
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Figure 16: Kaon (left) and antikaon (right) production dNc/dt and annihilation dNi/dt rates in a 
central Ni+Ni reaction at 1.0 AGeV. The solid line shows the time evolution of the central density. The 
figure is taken from |lllj . 



With this modification the statistical model delivers an accurate description of particle yields and 
ratios also at SIS energies (except of the r] yield which is systematically underestimated) |149t I15U] . 
The assumption of common freeze-out parameters for all hadron species stands, however, in clear 
contradiction to the dynamical transport calculations. In particular a freeze-out density of about po/4 
|149| I150j as obtained by the statistical model fit to the hadron multiplicities is much lower than the 
corresponding freeze-out density obtained from the dynamical models. 

The statistical model reproduces nonetheless several basic experimental features of the strangeness 
production at SIS. The canonical treatment leads to a quadratic Apart dependence of the kaon multiplic- 
ities Nx+ ~ y ^ ~ Apart which is close to the observed Apart dependence [21 1141j . Since pions are treated 
grand-canonical they scale linearly with V and Apart- Hence, the observed centrality dependence of the 
K^/tt~^ ratio [2] is obtained for free. The same holds for the Apart dependence of the K~ / ratio 
which is predicted to be flat in agreement with the Ni-|-Ni data. A fine tuning can be achieved by a 
variation of the temperature as a function of centrality |151j . This enforces, however, in a slight decrease 
of T for central reactions which contradicts the intuition of a larger energy deposit and a hotter fireball. 
In the larger Au-|-Au system the K~ / ratio has even be found to slightly decrease with centrality 
|152j which is at variance with a statistical interpretation but easy to understand from a microscopic 
point of view due to the large K~ absorption cross section. 

However, here appears another contradiction to the dynamical approaches: In contrast to the sta- 
tistical model based on bare hadron masses, the transport models need to introduce in-medium kaon 
masses in order to reproduce the experimental K~ /K^ ratio |154[ 1155] . This apparent contradiction 
between dynamical and statistical approaches has been addressed by Brown et al. |156j . From eqs. 
()52|) and ()53|) the / K~ ratio is approximately given by 

N^. ^ (rn,^ i e-("^-^^)/^ _ ^^^^ 



Nk- yrriK J e~^K-/'^ 

The experimental value ^ of Nk+/Nk- = 30 in Ni-|-Ni at 1.8 AGeV |157j enforces a chemical potential 
which is significantly smaller than the A mass in order to punish hyperon production relative to that 
of K~'s. A dropping K~ mass = rriK — aps compensates the /is dependence in (jH^ to large 

§In Refs. | 15Ulll5l] the earlier and smaller ratio Nk+ /N^- — 21±9 from | 141| has been used. 
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extent |156j and enables one to reproduce the experimental ratio with much larger freeze-out densities 
p ~ 1 2po, depending on the strength of the mass reduction. Such values for the freeze-out density 
are now in agreement with the estimates from transport models. 

However, the problem is not completely resolved since a more recent analysis |158j comes to the 
conclusion that the usage of in-medium masses obtained from coupled channel G-matrix calculations 
requires extremely low values for p and T, i.e. T < 34 MeV and p < 0.02po- Such freeze-out parameters 
seem, however, to be far from reality. When the same in-medium potentials are, on the other hand, 
applied within dynamical transport calculations [121] the K~ / ratio is fairly well reproduced. From 
this analysis one has to conclude that the statistical approach is at SIS energies not applicable to this 
observable, at least if one accounts also for in-medium effects. 



3.3.4 Thermal freeze-out 



In contrast to the chemical kaon freeze-out which takes place early, the kaon momentum distributions 
are strongly influenced by final state interactions. Both, elastic KN scattering as well as the propagation 
in a repulsive mean field from Coulomb and strong interactions makes the spectra harder |122l I124j and 
infiuences the collective in-plane li^ IM} I159| ll6Uj and out-of-plane fiow |123U16lj . In the case of the 
K~^s thermal and chemical freeze-out are not clearly separated since elastic and strangeness exchange 
cross sections are of the same order of magnitude. Propagation in an attractive mean field tends here 
to make the spectra softer. 

The infiuence of the mean field on the shape of the spectra has e.g. been studied in Refs. 
I12H 1124] . Both, a repulsive potential as well as an attractive K~ potential lead to deviations from 
a thermal spectrum which can be characterized by a radial fiow component. In |124j the transverse 
mass (niT = a/pt^ + ^k) spectra from QMD transport calculations were fitted by a radially expanding 
thermal source 

~ e-^--^\^^E - ^aT{— + 1) + {aTf — }^ ^ (55) 



d(j)dymTdmT p^ p 



with a common radial kaon velocity /3 = v/c. In ()55p E = rriT coshy, p = \/ "Pt^ + m^-^ sinh^y, 
a = 'jPp/T and 7 is the Lorentz factor. The repulsive potential leads to a depletion of the low mass rriT 
spectrum and creates a characteristic 'shoulder- arm' shape which gives rise to a radial fiow component. 
E.g. in typical Au-|-Au reactions at 1 AGeV the collective motion of the kaons due to the mean field has 
been found to be about 20% of the thermal motion (/5coii. — 0.1) |124j . K~ mesons exhibit an anlogous 
collective motion in the radial direction. However, the attractive potential leads here to a characteristic 
'concave' structure in the transverse mass spectrum |124[I121] which was attributed to a 'virtual' radial 
fiow in |124j . Thus a collective motion could in principle be distinguished from the thermal motion 
which would, however, requires to measure low mass spectra with high precision. 

Experimental evidence for different kinetic freeze-out of kaons and antikaons has been reported in 
|152j where in the same reaction (Au+Au at 1.5 AGeV) K~ spectra were found to be systematically 
softer than those of the K~^^s (see Fig. ITTj) . Also the angular distributions are qualitatively different, 
i.e. strongly forward-backward peaked for and much more isotropic for K~ |152j . Both facts 
indicate that kaons decouple early from the reactions dynamics while antikaons suffer strong final state 
interactions. 
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Figure 17: Inverse slope parameters T of K"^ spectra in Au+Au reactions at 1.5 AGeV measured by 
KaoS jl52j . The Figure is taken from |153j . 

4 Strangeness production in transport models 
4.1 Mean field dynamics for strange particles 

Semi-classical transport equations of a Boltzmann type can be derived from first principle quantum field 
theory. The starting point is the non-equilibrium real time Green's functions formalism |162j which leads 
after truncation of the Green's functions hierarchy at the two-body level, to the famous Kadanoff-Baym 
equations |163j for the one-body Green function. After a gradient expansion of Wigner transforms 
one arrives finally at the semi-classical kinetic equation. The kinetic equation is composed by a drift 
term, the Vlasov part, which propagates the one-body density and a collision term which accounts for 
two-body scattering. Derivations of the complete kinetic equation in full beauty and complexity are 
involved and can be found in the literature [201 12H1 11^31 11641 \lbb\ I166j for both, fermions and bosons. 
To set the context the relevant steps are briefly reviewed: 

Starting from the real time Schwinger-Keldysh formalism one obtains a matrix G of Green func- 
tions, (anti-) chronologically ordered Green functions and correlation functions G^'^ . The correlation 
functions are defined as G'>(1,1') = -i < ^(l)^(l') > and G'<(1,1') = i < "^(l)^(l') > using the 
notation 1 = (ti,xi). The quantity of interest is because in the limit ti = t'^ it corresponds to a 
density. The Dyson equation in non-equilibrium is given as a matrix equation 



with D(l, 1') = 5(1 - l')(n/.5f - m) for spin-1/2 fermions and D{1, V) = -5(1 - V){di^d'^ + m^) for 
bosons. The matrix of self-energies in Eq. ()56|) contains all higher order correlations originating 
from the higher order Green functions of the Schwinger-Keldysh hierarchy. When the hierarchy is 
truncated at the two-body level this leads to the T-matrix approximation summing all two-body ladder 
correlations. The kinetic equation for the correlation function is obtained by subtracting from the 
Dyson equation (j^UI) its adjoint. A Wigner transformation then allows to represent the kinetic equation 
in phase space, i.e. x — /c-space, rather than in coordinate space. An essential step is the truncation 
of the gradient expansion of the Wigner transform of products retaining only terms of first order in 
which neglect memory terms. The self-energy, Eq. (j56|) . is decomposed into scalar and vector parts 



^In this subsection fermions or bosons are not distinguished in the notation for their self-energy. Otherwise a bosonic 
self-energy is denoted by IT. 




(56) 



(57) 
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and the real part of S is used to define effective masses and kinetic momenta 

m* = M + ReEt{x,k) , /c^ = A;^ + /?eS+(x, A;) (58) 

of the dressed particles in the nuclear medium. S+ is the retarded self-energy constructed by the 
difference of the corresponding correlation functions 6*^(1,1') = 0{±{ti — ti')) — G^(l,l')] 

|167j . The Dirac structure of the correlation functions G^'^ can be separated off by a decomposition 
into a scalar spectral function a, a scalar distribution function / and the projector onto positive energy 
states A"*" 



G<{x,k) oc iA'^ a{x,k)f{x,k) (59) 
G>{x,k) oc -iA-^ a{x,k)[lT f{x,k)] . (60) 

The =F sign in Eq. ()60|) stands for fermions or bosons, respectivly. In the case of nucleons the projector 
is given by = {/c* + m*) /2m*. In an essential approximation which concerns the the spectral 
properties of the hadrons is the quasiparticle approximation. It is valid in the limit of a small imaginary 
part of the self energy (JmS+ << i?eS+). The spectral function then reduces to the mass shell 
constraint a{x,k) = 27r5 (/c*^ — m*^)20(A;Q) which sets the energy on the mass shell k^ = E*{x,\<.) = 
Vk*^ _l_ j^*2_ Thus, the number of variables of the distribution function f{x, k) is reduced from eight to 
seven 

a{x, k)f{x, k) = 2TT6[k*^ - m*2]2e(A;*)/(x, k) (61) 

which simplifies considerably practical implementations. As can be seen from Fig. El the quasiparticle 
approximation is well justified for kaons but less applicable to antikaons. 

In many cases the kinetic equation is further treated in the Hartree approximation which implies to 
neglect the explicit momentum dependence of the mean field, i.e. i^eS"*" = i?eE^(a;). Then the resulting 
kinetic equation can completely be formulated in terms of kinetic momenta instead of canonical momenta 

m\ 

[k*^d; + {klF^^ + m*9>*) d^;] {af){x, F) 

1 f ^2 k^ ki\. . , \ / 7 \ / 7 \ / 7\TTr/TJl7J\ 

X (27r) V (k + k^-ks- k,) [fix, ks)f{x, k,) (1 T fix, k)) (1 t fix, k^)) 

- fix,k)fix,k2)ilTfix,ks))ilTfix,k,))] . (62) 

Eq. (|62|) resembles the well known transport equation of a Boltzmann-Uehling-Uhlenbeck type. The 
left hand side is a drift term driven by the mean field via the kinetic momenta k*, the field strength 
tensor F'^^(x) = d^ReE'^'^ix) - d^ReEj/ix), and the effective mass m*. The right hand side is a 
collision integral which contains the transition rate W or equally the in-medium cross section given by 
(fc* + k2)'^da/dflik,k2) = W ikk2\k3k4/) . The colhsion term contains Pauli-blocking, respectively Bose 
enhancement factors (1 =F /) for the final states. 

As soon as one has to deal with different particle species which are interacting, e.g. via production 
and absorption processes, the transport problem becomes a coupled channel problem. This means that 
one has to solve a seperate transport equation for each degree of freedom, but these are coupled to each 
other by their collisions intergrals and the mean fields. 



4.1.1 Transport equation for kaons 

A full transport equation for kaons can in principle be derived as outlined above. However, since we are 
only interested in the structure of the relativistic equations of motion, we restrict the discussion in the 
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following to the mean field level and the left hand side of the transport equation (jU^ . i.e. the Vlasov 
equation. 

The Vlasov equation can easily be derived from the Klein-Gordon equation for the kaon field. For 
simplicity we restrict the discussion to mesons since the derivation for K~'s is fully analogous. 
Starting from the Klein-Gordon equation ()23|) the field equations for the kaon field and its adjoint read 



[{d^ + iV^y + mf] 0k(xi) = (63) 
[{d,-iV,f + n4]ct>],{x2) = Q . (64) 



In mean field approximation to the chiral 5'?7(3) Lagrangian as discussed in Chap. 2 a vector potential 
with alternating sign for kaons and antikaons, i.e. = ±^j„ for and an attractive scalar part 

TV 

occur. The latter can be absorbed into the effective mass = — ^p^Ps + V^iV^ which is equal 

for kaons and antikaons. Here one has already made use of the quasiparticle approximation since Eqs. 
contain only the real part of the kaon self-energy (j4H) . Multiplying equations (j63p and (j64j) 



by the respective adjoint field and subtracting them, yields 

+V^{X2)V^'{X2) + m*i{xi) - m*^{x2)] 0(xi)0*(x2) = . (65) 
Introducing center-of-mass and relative coordinates 

X = -[Xi + X2) , r = Xi- X2 



one has 



The fields are now expanded to first order in gradients 

V,ix^)d^^ + V,ix2)d^,^ = cosh ■ 9.) V,ix)d!i + 2 sinh ■ 9.) V,ix)di^ 

^ V^{x)d^ + (r ■ d,y^ix)) di! (66) 

where r ■ dx denotes a four- vector product. Similar one obtains 

"^K (^i) ~ "^K (^2) ~ 2m^(x) (r ■ 9a;m^(x)) 
V^{xi)V^{xi)-V^{x2)V^{x2) ^ 2V^{x){r-dRV''{x)) . (67) 

Inserting Eqs. (jUUj) . into (jHH) gives 

+ ^{v,ix)d^ + (r ■ dxV,{x)) d^^ + mUx) (r ■ d^mUx)) 
-V,{x) (r ■ dxV^ix))] (x + 0* (x - = . (68) 

Integration of Eq. ()68|1 over J d'^re^^'^ yields the Vlasov equation 

[k*^d; + {Kd^V^ - V^d^V^ + mld^^ml) 93 {af){x, k) = (69) 

for the scalar phase-space distribution f{x,k). The phase-space distribution is defined as the Wigner 
transform of the density matrix 00* 

ifc.. , / r\ ,^ / r 



fix, k) = I d're''-^<P [x + - j 0* i^x - -) . (70) 
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Introducing the field strength tensor F^^ = d^V^ — d'^V^ and kinetic momenta k*^ = V^, Eq. (jUUI) 
can now be written as 

[k^^dl + iklF^^ + m],d^,ml)d';]f{xX) = ^ ■ (71) 

The corresponding equation for K~ mesons is obtained by replacing in Eq. (fTTj) the field strength tensor 
F^'" with F''^ = -F'"''. 

Equation (|7T|l has exactly the same structure as the relativistic Vlasov equation for nucleons which is 
obtained from the projection onto the scalar component of the Wigner density matrix for spin-| fermions 
in the spinor representation. Thus, in the semi-classical limit the information of the fermionic/bosonic 
character of the particles is lost at the mean field level. In the collision integral differences are retained, 
e.g. fermionic Pauli-blocking factors are replaced by Bose-enhancement factors. 

4.1.2 Equations of motion 

From the Vlasov eq. ()7H) classical equations of motion are obtained in the quasi-particle limit which 
implies to put the phase space distribution on mass shell 

/(x,F) = /(q,t;k*)5(fc*-F*(q,k)) , (72) 

with 



E* = -^k*2 + m*^ . (73) 
Inserting ()72p into ()71|) yields the covariant equations of motion 

-f- = = u^' (74) 

dr mj^ 

^ = ^F^^ + d^m*^ . (75) 

Again the classical equations of motion are identical to those for nucleons, obtained e.g., in relativistic 
mean field theory [JIlJ. Only the different structure of the effective mass reflects the bosonic character 
of the kaons. 

Alternatively the equations of motion can also be obtained from the dispersion relation E = ko = 
E* ±Vo through the Hamilton equations 

dq_dE dk__dE 

ltt~~dk ' ~dl~ ~~dq ' ^ ^ 



Thus one obtains 



rfq _ k* 
'dt ~ W 



(77) 



dt E* dq ^ dq E* dq ' ^ ' 

Eq. ()78|) provide the non-covariant version of (fTSj) . The term proportional to the space- like component 
of the vector potential gives rise to a momentum dependence in Eq. (fTHj) which can be attributed to a 
Lorentz force term. This structure becomes more evident when Eq. (|78|l is rewritten in terms of kinetic 
momenta 

dk* E* dm*^ dv^ ^ k* ( d \ dUK ^ dV, 

= — ± — 7^ ^ = — K — ivj— — (79) 

dt oq oq E* \dq J dq dq 
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Figure 18: Action of the Lorentz force between nucleons {qqq) and kaons. The force acts opposite for 
K~{sq) and K~^{sq). 



with V = k*/E* the kaon velocity. The appearance of the velocity dependent (v = k*/E*) Lorentz 
force in Eqs. fj75|) -()79 |) is a genuine feature of the relativistic kaon dynamics as soon as a vector field is 
involved. In addition to the trivial |k| momentum dependence of the optical potential the k dependent 
part of the interaction appears in reference frames where the spatial components V and k do not vanish. 
Hence, in nuclear matter at rest or p+A reactions no Lorentz force is present. The importance of this 
interaction for the nucleon dynamics in heavy ion collisions is well known since a long time, see e.g. 
|l(i8j . For the kaon dynamics it was first discussed by Fuchs et al. [46^ . The physical picture behind is 
completely analogous to electrodynamics: assuming that the interaction between kaons and nucleons is 
determined by the non-strange quark content of the hadrons, two parallel currents are attracted, two 
anti-parallel currents are repelled. Hence the Lorentz force acts opposite for K~{sq) and K^{sq) as 
illustrated in FiglT^ In heavy ion reactions kaons are primordially produced in the early non-equilibrium 
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Figure 19: Transverse flow in 1.93 AGeV ^^Ni + ^^Ni reactions at impact parameter b< 4 fm. 
The calculations are performed without and including the kaon mean field dynamics. In the latter case 
results without and including the contribution from the Lorentz-force (LF) are shown. 

phase. In contrast to p+A reactions, there the baryon vector current is in general locally non-zero even 
in the cm. frame of the colliding nuclei. Thus a Lorentz force will be acting on the kaon. As discussed 
in jin] the Lorentz force counterbalances to large extent the infiuence of the time-like Vq component of 
the vector field acting on the mesons. This cancellation can easily be understood from Eq. ()78|). 
The vector field is proportional to the baryon current = {pb,upb) where u denotes the streaming 
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velocity of the surrounding nucleons. If u is locally constant, then the total contribution of the vector 
field in Eq. (I78|) can be written as =Fo72 (1 ~ |v||u| cosG) dpB/d(\. Now the angle G between the kaon 
and the baryon streaming velocities determines the influence of the Lorentz force. For = one obtains 
an attraction (repulsion) for (-^ ) and the opposite in the case of = 180°. 

Such dynamics have e.g. consequences for the kaon transverse flow. This is illustrated in Fig. IT^ 
where the transverse flow in Ni+Ni reactions at 1.93 AGeV is shown for various cases. It is seen 
that the strongly repulsive static potential Uk in (fT9|l tends to push the kaons away from the 
spectator matter, leading to a strong anti-flow around mid-rapidity. The Lorentz force pulls the kaons 
back to the spectator matter, resulting finally in a flow pattern which is close to that obtained without 
any in-medium potential. However, the magnitude of the cancellation effect and the final flow pattern 
depends also on the strength of the vector potential and the interplay between repulsive vector and 
the attractive scalar potential. In the early works on kaon flow [33 I159j the Lorentz force has been 
disregarded which led to an overestimation of the vector repulsion. One has, however, to keep in mind 
that correlations beyond mean field lead to an explicit momentum dependence of the scalar and vector 
potentials. Such an explicit momentum dependence reduces e.g. in the case of the nucleons the size 
of the scalar and vector potentials and, correspondingly the Lorentz force which softens the nucleon 
optical potential. A detailed discussion of the kaon flow and a comparison to data will be performed in 
chapter 5. 

4.2 OfF-shell transport 

Essential for the validity of the classical equations of motion is the quasi-particle approximation (QPA) 
which assumes that the spectral strength of a hadron is concentrated around its quasi-particle pole. 
Since kaons have practically zero width this condition is readily fulfilled in the vacuum. Particle widths 
can, however, dramatically change in a dense hadronic environment. To first order in density the in- 
medium width of a hadron in nuclear matter can be estimated by the collision width r*°* = F™'^ + F™", 

r-ii = 7wps , (80) 

with V the hadron velocity relative to the surrounding matter and a the total hadron-nucleon cross 
section. Due to their long mean free path the collisional broadening of kaons is small. The kaon- 
nucleon cross section is dominated by its elastic part which is about 10 mb and the majority of kaons is 
produced with relatively small velocity. Thus for i^'+'s in-medium widths can be expected to be small 
and the quasi-particle-approximation appears to be justified. For antikaons the situation is completely 
different. The strong coupling of the K~N system to the A(1405) and S(1385) hyperons leads not only 
to a shift of the quasi-particle pole but also to a strong distortion of the K~ spectral function [HIl EZ] 
(see Fig. Ej). The application of mean field dynamics is therefore questionable. 

First attempts to account for in-medium spectral properties in a more consistent way within the 
framework of semi-classical transport models were already made some time ago. A fully consistent 
treatment of the off-shell dynamics, i.e. a solution of the full Kadanoff-Baym equations has up to now 
only been performed for toy models and simplified geometries. Danielewicz jTHI was the first to solve 
the quantum evolution equations for the correlation functions C^'^ (|59I60|) for an infinite cylindrically 
symmetric system using a simplified two-body potential of Gaussian form. The work of Kohler |169j 
who solved the Kadanoff-Baym equations under the same conditions on a Cartesian momentum space 
grid can be viewed as the first step towards a realization of quantum transport on a lattice. To develop 
a lattice quantum transport for non-uniform systems and realistic interactions will be one of the future 
challenges in theoretical heavy ion physics. 

On the other hand, substantial progress has been made in the recent years to map part of the off- 
shell dynamics on a modified test-particle formalism I166t I172t [T7^ . This allows to apply off-shell 
dynamics, although in a simplified form, to the complex space time evolution of a heavy ion reaction. 
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Figure 20: Left: momentum distribution in nuclear matter at po from a self-consistent off-shell transport 
(solid) is compared to a many-body calculation |17()j (dotted) and to data from electron scattering |171j . 
The figure is taken from |172j . Right: Off-shell dependence of the elastic in-medium np cross section 
at Po- Off-shell DBHF matrix elements are folded over nucleon spectral functions of different width. 
Results are from [173j. 

That such type of approach is able to describe essential features of nuclear correlations beyond mean 
field has e.g. been demonstrated in Ref . |172j : in the left part of Fig. EO the nucleon phase space 
distribution /(k) in nuclear matter resulting from a self-consistent off-shell calculation is compared 
to a microscopic many-body calculation |170j and to electron scattering data |171j . Although these 
results were obtained with adjusted on-shell scattering matrix elements the agreement with the realistic 
distribution function is remarkably good. 

The present knowledge of off-shell matrix elements is rather limited. Except of particular cases of 
interest, e.g nY ^ K^N matrix elements [HSlinSl; theoretical investigations are scarce. In Ref. |173j 
in-medium half-off-shell matrix elements for elastic nucleon-nucleon scattering have been determined 
from relativistic DBHF many-body calculations based on realistic meson exchange potentials. The 
results indicate that a smooth transition from on-shell to off-shell matrix elements is possible and the 
usage of on-shell cross sections is justified for A^A^ scattering. The right part of Fig. 1201 displays the 
elastic off-shell in- medium neutron-proton cross section at po, folded over nucleon spectral functions of 
different width, as a function of the cm. momentum. The mean A^A^ cross section is in general only 
little affected which justifies the standard on-shell transport approach for the overall reaction dynamics. 

The question to what degree a depletion of the Fermi surface due to particle-hole excitations and 
the high momentum tails of the nuclear spectral functions will affect subthreshold particle production 
is, however, not so obvious to answer. The high momentum tails correspond to deeply bound states 
which are off-shell and to treat such states in a standard transport approach like on-shell quasi-particles 
would violate energy-momentum conservation. The contribution of the nuclear short-range correlations 
to subthreshold production in p+A reactions have been estimated in (144J, based on the spectral 
distributions of Benhar et al. ,170^ shown in Fig. 1201 The removal energy for a high momentum 
state compensates the naively expected energy gain and the short-range correlations do therefore not 
significantly contribute to subthreshold particle production |144j . The situation changes, however, 
when the medium is heated up and high momentum particles become on-shell or when the spectral 
distributions of the produced hadrons themselves are broadened. 

Possible implications of off-shell dynamics connected to a spectral broadening of the kaons have first 
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Figure 21: Influence of liadron off-sliell propagation on spectra in Ni+Ni reactions at 1.8 AGeV. 
On-sliell cross sections fiave been used for tfie strangeness production and excfiange reactions in tlie 
RBUU calculation of [TTIj . 



been studied by Gassing and Juchem |174U12lj . These calculations were based on on-shell cross sections 
and hadronic spectral functions {N,^, A, , K~) which are exclusively determined by the collisional 
broadening (jHUjl of these particles. Nucleons and pions were found to be practically unaffected by the 
off-shell propagation while the high momentum tails of the spectra are enhanced. But even at 
subthreshold energies the effect on the yield was found to be moderate which justifies the standard 
quasi-particle approach and is in line with the observations of Ref. [144] . As expected, the off-shell 
effects were found to be much stronger for antikaons. K~ 's aquire a large spectral width in the medium 
and it is energetically favorable to produce them off-shell in the low momentum tails of their spectral 
functions. During the propagation through the medium these particles become on-shell as they reach 
the detector. Obviously this transition must take place in a controlled way which is guaranteed within 
consistent off-shell dynamics jl74l I12t| . As can be seen from Fig. 1^ where the corresponding 
spectra in Ni-|-Ni reactions at 1.8 AGeV are shown, the collisional broadening leads to an enhancement 
of the K~ yield by about a factor of 2-3 compared to the on-shell treatment. However, in particular for 
K~ the off-shell behavior of the in-medium strangeness exchange cross sections ttY < — > NK~ plays an 
important role |121j . Results based on off-shell cross sections will be discussed in the following chapter. 



4.3 Collisions 

Transport models treat particle production at subthreshold energies usually in a perturbative way. In 
each individual hadron-hadron collision the incident energy ^/s must be above the production threshold. 
However, such processes are rare and the corresponding cross sections lie in the vicinity of thresholds 
usually several orders of magnitude below the total cross sections. Since Monte-Garlo methods which 
are applied to solve the collision integral of the kinetic equation |175j select the actual reaction channel 
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i with probability 



tot ( g^jj channels 

meson production is at subthreshold energies an extremely rare process. To obtain nevertheless rea- 
sonable statistics, for such reactions a perturbative treatment is applied \l7i}\ \m\ I178j . In short, this 
means that the meson is produced if kinematically allowed. Phase space and Pauli-blocking factors for 
the final states are determined and the produced meson k gets the probability 

Y^K+X 

assigned. Y denotes the initial state and the sum runs over all possible final states. From the set of 
possible final states the actual reaction channel k, i.e. Y — > K~^Xk, is selected by Monte-Carlo. Thus, 
e.g. the total multiplicity is given 



Nk+ = J2 Pk . (83) 



all K+ 

Finally the initial particles are reset to their initial values and the cascade continues as if no such 
reaction would have taken place. The perturbatively produced mesons are propagated as 'virtual' 
particles parallel to the ongoing reaction dynamics in such a way that they are influenced by the bulk of 
particles, however, without any feed back to the global reaction dynamics. The perturbative treatment 
is justified as long as the considered particles are rare enough not to influence the average reaction 
dynamics. 



4.3.1 Elementary cross sections 

As discussed in the previous chapter, the elementary production mechanism for mesons can be 
divided into baryon induced reactions BB — > BYK^ {B stands either for a nucleon or a A-resonance 
and Y for a A or a S hyperon, respectively) and processes nB — > YK^ induced by pion absorption. 
Here and in the previous works of the Tiibingen group jlEl CUDl 1123 CSni EZHl IMIl CHIl for 
pion induced reactions the elementary cross sections of Tsushima et al. |llUj have been used. The 
cross sections are derived within the resonance model in Born approximation, including all baryonic 
resonances with masses below 2 GeV as intermediate states. In the meantime they are standardly and 
applied in most transport calculations [Iini CHI 123 EH CMl EHl HM EH] • 

For baryon induced channels there are presently several different parameterizations on the market: 
Those from Li and Ko |154j used by the Texas and Stony Brook groups; those given by Sibirtsev 
P7j which are used by the Giessen group 123 and have been applied in some of our previous 
works |lU9t 11591 1123| 1^: and more recent calculations within the resonance model by Tsushima et al. 
|lUm IIUH ll(J2j . In principle it would be desirable to base the cross sections for both, pion and baryon 
induced reactions on the same model |limilU2j . However, the resonance model |lU2j under-predicts the 
COSY-11 data [HZl for pp — > pAK^ at threshold. Thus, here and in |181l 1182] the cross sections from 
|97j were used for NN — > NA{Ti)K^ and those of ref. |lU2j for reactions involving nucleon resonances 
[NN — > AYK+, NA — > BYK+ and AA — > BYK+). Figure |221 compares the model predictions 
for pp and vrp channels to available data. In the vicinity of the threshold the nucleon-hyperon final 
state interaction is attractive and strong which tends to increase the cross section |188j . Since both, 
the OBE as well as the resonance model calculations do not account for FSI effects, the COSY data are 
under-predicted near threshold [57]. As pointed out in |189j the modification of the energy dependence 
of the cross section due to FSI is determined by the on-shell T-matrix. Since in the semi-classical 
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Figure 22: Elementary cross section for kaon production. Resonance model calculations from Tsushima 
et al. [niIllIII2] and OBE calculations from Sibirtsev ^ are compared to data CHSl EHEl CHII • 



transport models on-sliell potentials enter into the evaluation of the scattering amplitude, e.g. by shifts 
of thresholds, some FSI contribution is effectively included which justifies the usage of the cross sections 
from P7] . 

In contrast to the A^A^ channel, cross sections with an incoming A resonance NA — > BYK^ are 
not constrained by data. Here the largest differences between various parameterizations occur. Fig. 
ESI compares the resonance model results |lU2j with parameterizations based on isotopic relations for 
PP 1^ IMj . The latter cross sections turn out to be about one order of magnitude smaller. This 
difference is crucial since the A^A — > K^X is one of the primary production channels for kaons. 
It explains why there existed a significant inconsistency between transport calculations from different 
groups concerning the total kaon multiplicity in heavy ion reactions. Kaon yields obtained from the 
Giessen group [23 IHH turned out to be systematically smaller than those of the Texas and Stony 
Brook I154j . the Nantes and the Tiibingen |18H 1181?] groups. The source of this discrepancy, 
which is crucial with respect of an interpretation of corresponding data in terms of kaon in-medium 
potentials, could finally be traced back to the usage of the different cross sections for kaon production 
in A^A collisions |96j . 

Reactions with multi-particle final states, such as A^A^ — > NnYK, NN — > NttttYK, are subdom- 
inant at subthreshold energies. If taken into account they can either be treated by direct parameteriza- 
tions of experimental data |154j or by two-step processes via resonance production A^A^ — > RYK and 
subsequent decay R — > N{n ■ vr). Fig. ESI compares resonance model calculations |lU2j to data. The 
experimental values in Fig. |2Slare normalized to the iso-spin factors of the corresponding A — > Nn 
decays. The present calculations include reactions with a A resonance in the final state with cross 
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Figure 23: Cross sections involving A resonances. Upper part: for the experimentally unknown pA — > 
pYK^ cross sections resonance model calculations |l()2j are compared to the parameterizations of |25j . 
Lower part: Resonance model calculations |lU2j for pp — > ^YK are compared to pp — > NttYK data 



sections taken from |lU2j . 



4.3.2 Shift of thresholds 

The cross sections are obtained for free scattering. The incorporation of medium effects in the scattering 
process is a peculiar question, in particular since at SIS energies most kaons are created near threshold. A 
shift of the thresholds by in-medium potentials affects the production mechanism significantly. However, 
the treatment of the threshold conditions and the determination of the phase space of final states is 
in the medium a subtle problem. Since there exists up to now no unified description of relativistic 
and non-relativistic approaches we will discuss this point in more detail. In this context problems 
arise generally due to on-shell scattering of quasi-particles. In the presence of the medium mass-shell 
conditions are modified and phase space relations valid in free space have to modified as well. 

Let us for example consider baryon-baryon induced production BB — > BYK^ . In free space 
the momenta of the outgoing particles are distributed according to the 3-body phase space 

d$3(v^, m^, my, rriK) = (i$2(v^, m^, M)dM^$2(Ai', my, m^^) (84) 

with ^/s the center-of-mass energy of the initial baryons. The two-body phase space in Eq. (|5^ has 
the well known form 

^2Ws,mi,m2) = j= (85) 
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where 

is the momentum of the particles 1 and 2 in their center-of-mass frame. 

Concerning the in-medium description differences arise also between non-relativistic and relativistic 
approaches. In the relativistic case the mean field is generally composed by scalar and vector parts, a 
feature which is absent in purely non-relativistic approaches. However, also the kaon mean field is of a 
scalar-vector type structure which makes a non-relativistic treatment somewhat ambiguous. The reason 
lies in the fact that in most non-relativistic approaches canonical momenta fc^ are propagated whereas 
in relativistic approaches usually kinetic momenta A;* are used. To illustrate this effect we consider first 
nucleon-nucleon scattering. 

In relativistic dynamics, e.g. given by the au model of Quantum Hadron Dynamics jl7] the nucleon 
mean field is also composed by a scalar and a vector part. Like for kaons the vector field enters 
into the kinetic momenta k*^ = k^j, — and the scalar field into the effective mass m* = m -|- Eg of the 
particles. The dressed quasi-particles fulfill the mass shell condition 



kf - = , E* = Vk*2 + m*2 (87) 

and thus transport models are usually formulated in terms of these kinetic quasi-particle quantities. If 
two nucleons are scattered, energy-momentum conservation requires 

/ciM + k2f, = k[^ + 4^ (88) 

which is equivalent to the conservation of the kinetic quantities 

kifj, + = ki^ + (89) 

as long as the fields are density dependent but do not explicitely depend on momenta, i.e. S^^5 = 
^fi,s{p)- Then both conditions (jHHjl and (jHIIjl can be fulfilled simultaneously. This means that the phase 
space relations (jHiHHH|) can be used by replacing the bare quantities s,mB,my,mx by the effective 
quantities s*,m^,my,m^. 

However, mesons are created by associated strangeness production which leads automatically 
to a shift of the corresponding production thresholds in the medium, even when no modifications of 
the kaon properties are taken into account. The reason is the creation of the associated hyperon. Due 
to its reduced non-strange quark content the mean field of the hyperon should scale - at least in a 
simple SU{3) fiavor picture - with about a factor of 2/3 compared to the nucleon field. Such a scaling 
is in rough agreement with mean field calculations for hyper nuclei |19m I19H I192j . Now the baryon 
fields are no more conserved and therefore it is no more possible to conserve both, kinetic and canonical 
momenta simultaneously. From the derivation of kinetic equations fHl 1201 it is clear that in this case 
single particle energies 

E = kQ = E* + (90) 

and canonical momenta A;^ = {E, k) are the quantities which have to be conserved. Consequently, 
the usage of on-shell phase space relations (jSlHSHj) with the quantities s* ,m*Q,mY^m*j^ will lead to a 
violation of energy conservation. 

To overcome this problem and to make nonetheless use of on-shell relations it is useful to formulate 
the mass-shell conditions in terms of canonical momenta. This can be achieved with the help of the 
optical potential defined in ()29j) which allows to rewrite the in-medium dispersion relation as = 
k*i^ — m^ = kf^ — m^ — 2mKf/opt- Since Uopt is a Lorentz scalar it can be absorbed into an newly defined 
effective mass rfiK 



rriKip, k) = V mf^ + 2mKt/opt(p, k) (91) 
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which sets the canonical momenta on mass-shell 

= kf_m*^ = kl-ml . (92) 
The single particle energy follows from the dispersion relation ()27|1 written now as 

E = ^k2 + m| . (93) 

The threshold condition for production in baryon induced reactions reads now 

y/s > fflB + rhy + rfiK (94) 

with y/s the center-of-mass energy of the colliding baryons. The momenta of the outgoing particles 
are distributed according to the 3-body phase space 

d^3{y/s,rhB,rhY,mK) = d^2{Vs,rhB,M)dM^^2{M,mY,rh . (95) 

The integration over the mass distribution of the mB-,ffiY-,fhK system in Eqs. ()86|) and has to be 
performed numerically. The introduction of mx is thereby of practical use. However, in contrast to the 
quasi-particle mass which is in mean field approximation only density dependent, mx depends on 
density and momentum. This means that 3-body-phase-space of the final states and the corresponding 
final state fh masses for the final state hadrons have to determined self-consistently by iteration. Thus 
the production threshold depend also on the final momentum distribution. The usage of a simply shifted 
in-medium mass mx = frtK ± ap in eq. ()93|) neglects such a momentum dependence. 

The potential of the final state hyperons is usually determined under the assumption of SU(3) 
symmetry. The scaling with the non-strange quarks leads than to the reduction by 2/3 compared 
to the nucleon potential. Such a scaling has been found to be in reasonable agreement with the A 
dynamics in heavy ion reactions, in particular the A flow |179l IHHl and is also close to the value 
extracted from hyper nuclei |190t 11911 1192j . Thus one obtains an additional shift of the thresholds 
hy Uy — Ub = —\Ub- In the QMD calculations discussed in the work the shifts of the thresholds 
due to the different initial and final state in-medium potentials have been treated within a relativistic 
framework, i.e. Uy — Ub = ~|(Ss ~ ^o)- This allows well defined Lorentz transformations of the kaon 
and baryon mean fields concerning their scalar- vector structure. In correspondence with the soft/hard 
Skyrme forces a soft/hard version of the non-linear acj-model [133] with K=200 MeV and K=380 MeV 
are used to determine the reaction threshold. 

Besides the scaling of the hyperon optical potential also the momentum dependence of baryon 
optical potentials influences the production thresholds. The nucleon optical potential is repulsive at 
high momenta which are necessary for the incoming states to overcome the production threshold. 
The momenta of the final states are significantly smaller. Due to the less repulsive optical potential 
experienced by the final states additional energy is gained to overcome the production threshold. At 
subthreshold energies all these effects are essential and have to be taken into account in the transport 
calculations. 

4.3.3 Angular distributions 

For most reactions elementary cross sections are assumed to be isotropic. In some cases there exists 
experimental evidence for anisotropic angular distributions. However, from heavy ion spectra it is 
difficult to disentangle the possible sources for anisotropic kaon emission. One source are final state 
interactions such as KN rescattering and the mean field propagation [122| I123| IM] . For a reproduction 
of the strongly forward-backward peaked pattern reported in ^193, these sources turned out to 
be not sufficient and an additional forward-backward anisotropy was introduced into the pion induced 
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Figure 24: Inclusive spectra at 6'iab = 32°, 40° in C+C reactions at 2.0 AGeV. The calculations are 
performed with different BB — > BYK^ final state 3-body phase space distributions and compared to 
KaoS data |194j . The dotted curves refer to an isotropic 3-body phase space while the solid curves are 
obtained using the parameterization of Eq. (19) with an additional empirical cm. angular anisotropy 
(see text). The figure is taken from |182j . 



vrA^ — > YK^ channel in |123j . Instead of an ideal 3-body final state phase space in the baryon induced 
reactions A^A^ — > NY Li et al. jHj proposed an empirical parameterization of the form 

c/$3(v/s, m^, my, ffiK) = dWxiVs, m^, my, MK)dM]^^2{Vs - Mk, rhy, rfiB) (96) 

where the kaon momentum p is distributed according to 



'K 



k W k 



k \ k 

"'max / \ "'max 



1-7— ' (97) 



with fcmax = P*{y/s,'fnB + fiT'Yi'fnK) the maximal kaon momentum in the BB cm. frame. Mk = 
\/k^^~+~A^ in Eq. ()9(ij) . The parameterization of Eq. has been motivated by analyzing correspond- 
ing pp pAK^ data |195j . Compared to the ideal 3-body phase space it shifts the kaon spectrum 
towards lower momenta. The influence of different treatments is illustrated in Fig. 1^ where spectra 
in C-l-C collisions at 2.0 AGeV are shown. KaoS measured inclusive spectra at various cm. angles 
with high precision |194j . From the comparison to data one sees that an isotropic 3-body phase space 
in the BB BYK^ channel shifts the spectra to too high momenta. The empirical parameterization 
of Eq. (j^JTj) improves the situation but is still not fully sufficient to account for the angular asymme- 
try of the data |194j . A relatively good agreement can be achieved introducing an empirical angular 
dependence 

da (X {1 + a cos^ Oc.m.)d cos 9c.m. (98) 

in the elementary cross sections. An asymmetry parameter of a = 1.2 leads to slightly forward/backward 
peaked elementary BB — * BYK^ cross sections and the corresponding spectra shown in Fig. |^ are 
then well reproduced. 



4.4 Comparison of different transport models 

This subsection addresses the question how consistent the results of present transport models are. 
Differences may occur due to the use of different physical input such as e.g. elementary cross sections. 
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The type of the model, i.e. BUU or QMD should not be of relevance. However, the corresponding 
simulation codes are complex and sometimes based on different numerical and methodical solution 
techniques. Although physical observables should be independent on such technical questions one has 
to exclude them as possible sources of uncertainties as far as possible. This was the major goal of two 
workshops held in Trento 2001 and 2003 where all major groups doing transport model calculations in 
the SIS energy range participated. In a first round of homeworks the default versions of the codes were 
compared, in a second round further specifications were made for the comparison. The results of the 
second round have been published in |19(ij . 




Y-Y Y-Y Y-Y 

cm. cm. cm. 

Figure 25: rapidity distributions in central (b=l fm) Au+Au reactions at 0.96 and 1.48 AGeV and 
Ni+Ni reactions at 1.93 AGeV from various transport models: RBUU (Texas, open circles), RBUU 
(Giessen, full circles), IQMD (Nantes, full squares) and QMD (Tiibingen, full triangles). 

In the following I compare results for pion and kaon production from four state of the art transport 
models which participated in that benchmark test. This selection represents only a subset of the 
models compared in |196j and is motivated the fact that most of the publications on subthreshold 
production over the last five to ten years are based on one of these codes. Of course each of the models 
experienced a steady devellopment during this period but the selection represents well established 
transport models used to describe production at subthreshold energies: the RBUU model of the 
Texas A&M group [HSl IM|, the RBUU model of the Giessen group [23, the IQMD model of the 
Nantes group and the QMD model of the Tiibingen group |18H \1H'J\ . Additional reasons why the 
present discussion is restricted to the four mentioned approaches instead of all models presented in |196j 
are the following ones: the BUU model of Danielewicz [198J does not contain explicit kaonic degrees of 
freedom, UrQMD (Reiter) [21] treats kaon production non-perturbatively through resonance excitations 
and is not well suited for an application at subthreshold energies, the Munich/Catania RBUU model 
(Gaitanos) |199j uses the Tiibingen kaon package and gives therefore similar results. The status of the 
Rossendorf /Budapest BUU model (Barz/Wolf) is hard to estimate since there exist no publications 
on production based on this model from recent years. Finally the Giessen BUU model (Larionov) 
included for the benchmark test all baryonic resonances with masses below 2 GeV under the assumption 
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that the excited A* and A^* resonances contribute to the production with the same cross sections 
as the A and the nucleon, respectively. Hence a comparision to the four models which include only the 
A(1232) and A'^*(1440) for kaon production can be missleading, in particular since the contributions 
from higher excitated nucleon resonances to the kaon production will be suppressed by a quenching 
of the elementary resonance production cross sections at finite density. The resonance quenching is 
taken into account in |2UUj where the Giessen BUU model shows a much better agreement with the 
four selected models than the benchmark test from |196j would suggest. For the benchmark results of 
this model we refer the reader to However, the model is included in Fig. ITfl where the default 

versions of the various codes are compared. 

The bechmark tests were performed for three different systems, Au+Au at 0.96 and 1.48 AGeV 
and Ni+Ni at 1.93 AGeV, all at impact parameter b=l fm. For this comparison a soft nuclear EOS 
with momentum dependent forces was applied. Except of the Giessen RBUU model" the results were 
obtained with a constant A width Fa = 120 MeV, respectively a constant lifetime r = 1/Fa. A 
constant resonance life time is unphysical and not used in the default versions of the codes but it 
simplifies such a comparison to some extent, in particular when predictions for pion production are 
considered. In QMD** and IQMD only the A(1232) and A^*(1440) resonances have been included, 
while the Giessen RBUU includes additionally the A^*(1535). Fig. 1^ shows the resulting vr" and vr"*" 
rapidity distributions: For Au+Au at 0.96 and 1.48 AGeV the agreement between the different codes 
is quite satisfactory, i.e. within a 5-10% level. Only at the highest considered energy of 1.93 stronger 
deviations are visuable since the Giessen calculations includes one higher lying resonance. In this context 
it should be mentioned that the vr^ yields from Texas RBBU shown in Fig. |2S1 differ from those of |196j . 
There isospin averaged rapidity distributions are shown while here the isospin dependence is taken into 
account by the isobar model. The default versions of the codes use different descriptions of the energy 
dependence of the resonance life times. As discussed in [196 the two RBUU models under consideration 
use r oc l/rA(|k|) with an energy dependent width. The QMD [20] and IQMD |201j models apply the 
time delay description |202j where the lifetime is obtained from the phase shift r = 2dS{E)/dE which 
results in a Breit-Wigner form r(/i) = 4:'7ifidW/\{fi)/dfi'^ with 

'i^'-M = -r^^^¥#fT^ (99) 

In Eq. fl99|) mA is the resonance pole mass and fi the running mass. The pion yields obtained by such 
a description are lower than those shown in Fig. 1221(15-20% for QMD and about 30% for IQMD) while 
they do not much change when the constant width is replaced by an energy dependent width in the 
RBUU calculations. 

Though pion yields agree well the situation is not yet as clear concerning the kaon production. Fig. 
|2ni shows the rapidity distributions obtained without in- medium effects and including an in-medium 
kaon potential. First we will discuss the results without medium effects: 

The results are based on comparable, however, not completely identical input. The Giessen RBUU, 
IQMD and QMD calculations are based on exactly the same set of elementary cross sections, i.e. the 
parameterizations of [HZI for A^A^ — > NYK^ and those of |102j for A^A — > NYK^ (in contrast to 
previous publications of the Giessen group). Pion induced reactions are based on the cross sections of 
|110j . Hence the RBUU results from Giessen agree relatively well with those of the Tiibingen QMD, i.e. 
on a 10-15% level. Compared to these two models the Nantes IQMD leads to an about 30% higher kaon 
yield, in particular at higher energies. This discrepancy can probably not be traced back to the level 
of elementary input and further clarification is needed. The Texas RBUU yields lie significantly below 



"Private communication with E.L. Bratkovskaya. 
**An extended version of the QMD model which includes all nucleon resonances with masses below 2 GeV was used 
for the description of dilepton production at SIS energies [2^1 ■ 
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Figure 26: rapidity distributions in central (b=l fm) Au+Au reactions at 0.96 and 1.48 AGeV 
and Ni-|-Ni reactions at 1.93 AGeV from various transport models |203j : RBUU (Texas, open circles), 
RBUU (Giessen, full circles), IQMD (Nantes, full squares) and QMD (Tiibingen, full triangles). The 
upper curves show results without kaon in-medium potentials while the lower curves include potentials. 

the other ones. This discrepancy is, however, understandable since the Tsushima cross sections |102j 
have been used for the A^A^ — ^ NYK^ channel which are at threshold about one order of magnitude 
smaller than those of [HZl, see Fig. [221 Consequently, the suppression of the RBUU Texas yields is 
most pronounced at the lowest energy of 0.96 AGeV. Here the impact of the elementary input is clearly 
reflected. 

All models show the same qualitative in-medium effect, namely a sizable reduction of the yield. 
One has, however, to be aware that the different models use repulsive in-medium potentials of different 
strength and treat their momentum in different ways. When the kaon potential is parameterized in the 
form of an effective mass 

"^K = "^k(1 + ap/po) (100) 

the coefficients for the results shown in Fig. |2niare a = 0.07 (QMD; MFT ChPT-|-corr. potential), 0.075 
(IQMD), 0.04 (RBUU Giessen) and 0.04 (RBUU Texas) This leads automatically to different 

predictions for the kaon yields. According to the weakest potential the relative suppression of the kaon 
yields is least pronounced in the RBUU calculations. There exist, however, not only differences in the 
size of the potentials but also in the way how the kaon mean field is treated. The present QMD ^[ I181j 
calculations (and those of |200j ) are the only ones which account for the full covariant structure including 
the Lorentz force from the spatial components of the vector field ([7^1) . RBUU Texas and IQMD neglect 
these contributions and in the RBUU Giessen model the total kaon potential accroding to expression 

(unni). 

From the above comparison it seems surprising that all four models are able to fit the experimental 
K+ data for Ni+Ni at 1.93 AGeV [M EHl 1203 IMI (see also Fig. M)- Therefore in Fig. |27| the 
predictions of the default versions of the codes are shown, again for Ni+Ni at 1.93 AGeV at impact 
parameter b=l fm. This figure contains in addition the predictions of the BUU model of Larionov 
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Figure 27: rapidity distributions in central (b=l fm) Ni+Ni reactions at 1.93 AGeV from various 
transport models in their default versions: RBUU (Texas, open circles), RBUU (Giessen, full circles), 
IQMD (Nantes, full squares), QMD (Tiibingen, full triangles) and BUU (Giessen). The left figure shows 
results without kaon in-medium potentials while the right one includes potentials. 



and Mosel (BUU Giessen) which is also able to reproduce experimental kaon multiplicities when a 
quenching of higher nucleon resonances is taken into account 200 . This model uses exactly the same 
in-medium potential as is used in the present QMD calculations. From Fig. EZI it becomes clear 
that the unphysical constraint of a constant A width (which violates detailed balance in the codes 
if pion absorption is not modified accordingly) has a strong influence on the kaon yields. Similar to 
the pion yields the kaon yields are smaller in the QMD/IQMD default versions, in particular when in- 
medium potentials are used. For the Texas RBUU the behavior is opposite, i.e. the kaon yields are now 
enhanced. Considering the case without potential one sees that the models group into two fractions: 
(IQMD, BUU (Giessen)) and (QMD, RBUU (Texas), RBUU (Giessen)). The latter ones yield an about 
30-40% smaller kaon yield. The reason for this discrepancy is still an open question which has to be 
settled in future. 

When in- medium potentials are taken into account, the QMD and BUU rapidity distribution are 
slightly broader than the other ones. This is due to the momentum dependence originating from the 
relativistic Lorentz force which is only accounted for in these two models. The relative potential effect is 
now larger for the Giessen RBUU compared to Fig. EDI since the stronger chiral RHA potential of j204j 
is applied, corresponding to a ^ 0.055 in the above parameterization. When in-medium potentials are 
applied, the uncertainty of the present model calculations, is of similar magnitude than experimental 
error bars. However, as mentioned above, the agreement of the codes in their bare versions is still not 
completely satisfactory and requires further efforts to improve on this. 

For a compilation of the predictions for K~ production, which turned out to be still burdened with 
much higher uncertainties, we refer the reader to Ref. |196j . 
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5 Probing in-medium kaon potentials 



In-medium kaon potentials shift the single particle energies (I26p and the thresholds for kaon produc- 
tion (|9HI94p . The measurement of kaon multiplicities should therefore provide a direct access to such 
in-medium potentials and, consequently probe an expected partial restoration of chiral symmetry at 
supra- normal nuclear densities |122j . For mesons the repulsive mean field reduces the yields while 
the situation is opposite for K~ where the yields should be significantly enhanced by the attractive 
potential. The hope is to extract information about the existence and size of such potentials from 
heavy ion reactions by transport calculations. However, the situation turned out to be more complex 
than originally expected: 

Despite the fact that the kaon mean field has indeed a strong influence on the multiplicities their 
absolute values depend as well on the elementary cross sections. Uncertainties due to an incomplete 
knowledge of the cross sections turned out to be of the same order as the potential effects. However, in 
the meantime a more or less consistent picture has emerged concerning the mesons. For antikaons 
the situation is less satisfying since cross sections are known with less precision and the K~ yield itself 
is strongly coupled to the production rate via strangeness exchange reactions. 

As a way out of this dilemma there were also strong attempts to extract the potential from dynam- 
ical observables which are altered by the in-medium optical potential while uncertainties in the total 
production rates drop out. Promising probes are, e.g., collective flow pattern. 



5.1 Total Yields 

For mesons the various transport models now provide a relatively consistent picture what concerns 
the net potential effect on the kaon multiplicities. The repulsive mean field leads to a reduction of the 
yields by about 30-50%, depending on the actual strength of the potential, the system size and the 
energy of the reaction. The magnitude of the reduction within different transport model realizations 
can be read off from Fig. 123 Figures. EHl to IHUl demonstrate the energy and system size dependence and 
the conclusion which can be drawn from the comparison to data. Fig. |2H1 shows the potential effect in 
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Figure 28: Influence of in-medium potential on the kaon yields. Multiplicities obtained with and 
w/o in-medium potential are shown for central (b=0) Au+Au and C+C reactions. 
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central Au+Au and C+C reactions as a function of beam energy. Throughout this work the calculations 
which include an in-medium potential are based on the mean field proposed by Brown and Rho , 
denoted in Figure Q as MFT ChPT+corr., which has been derived from ChPT and includes effectively 
higher order corrections beyond mean field. The reduction of the yield due to the repulsive potential 
is, as expected, slightly larger in heavy systems than in light systems and most pronounced at energies 
far below threshold. The excitation functions of the cross sections in inclusive Au+Au and C+C 




Figure 29: The excitation functions in Au+Au (scaled by 10~^) and C+C reactions are compared 
to data from KaoS |2U5| 1142] . Calculations include an in-medium kaon potential. For C+C also results 
w/o in-medium potential are shown. 

are shown in Fig 1221 Calculations were performed with fomax = 11 fm for Au+Au and 6max = 5 fm for 
C+C and are normalized to the experimental reaction cross sections. The comparison to data from 
KaoS |142| l2U5j clearly supports the existence of such a repulsive potential. 

The same conclusion has been obtained by other groups, e.g. the Texas and Stony Brook groups 
I154[ IT^ l2U6j and the Nantes group jHEj • For some time these findings were in contradiction to 
calculations from the Giessen group jM] I116t 123 1121j where the in-medium suppression of the yield 
was found to be of similar size but corresponding data where then underpredicted. The reason for this 
discrepancy could be traced back to the usage of different elementary cross sections, in particular for 
the A^A 1-^ NYK^ channels (see also Fig. The cross sections used by the Giessen group have been 
derived from pp pAK^ by isotopic relations while the calculations of the present work (Tiibingen 
group) are based on the cross sections of Tsushima et al. |2U7j for the NA channel. The latter ones 
have a better theoretical foundation and are now standardly used in transport calculations. However, 
in earlier calculations the uncertainty due to this channel (not constrained by data) was of the same 
order as the net potential effect. When model calculations are based on comparable sets of elementary 
cross sections the results of the various transport models are rather consistent, as can be seen from 
Figs. EiandEOl 

Fig. IHUl compares the rapidity distributions in Ni+Ni reactions at 1.93 AGeV to data from FOPI 
|135j and KaoS [157] *. Again the description of the data requires the repulsive mean field. The same 
conclusions are obtained from IQMD and the RBUU calculations of Mishra et al. |2U4j , now with 

*y(o) denotes the center-of-mass rapidity normalized to the projectile rapidity Y^°^ = (F/Fproj)c.m. = 2Yiab/yboam — 1- 
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the A^A i-H> XK^ cross sections of Tsushima et al. j2DTj . The RBUU calculations shown here are based 
on a chiral mean field evaluated in relativistic Hartree approximation with Skn = 450 MeV |2U4j . The 
range term is included and thus the mean field is close that one used in the present QMD calculations. 
The IQMD calculations are based on the RMF kaon optical potential of Schaffner et al. ISH] which is also 
of similar strength than the chiral mean field. The suppression of the kaon production by the repulsive 




Figure 30: rapidity distributions in Ni+Ni reactions at 1.93 AGeV obtained including and w/o an 
in-medium kaon potential. QMD, IQMD [^0] and RBUU (Giessen from |204j . Texas from [197J ) and 
calculations are compared to data from FOPI |135j and KaoS |157j . 

in-medium potential is most pronounced at mid-rapidity which is understandable from kinematical 
reasons. Since subthreshold particle production takes dominantly place around mid-rapidity the kaons 
are produced close to threshold where in- medium shifts have maximal impact. As already refiected in 
Fig. [2nithe in- medium effect is slightly stronger in the IQMD calculations [96j compared to QMD and 
RBUU. Both RBUU calculations use a slightly weaker in-medium potential (RBUU Giessen is based on 
the chiral RHA potential of ^204 and RBUU Texas on the empirical potential of |197j ). In the QMD 
calculation the in-medium potential leads not only to a suppression but also to a slight broadening 
of the rapidity distribution which is due to the covariant dynamics including the Lorentz force. This 
contribution is absent in the IQMD and RBUU calculations. Although the theoretical descriptions show 
still same variance they allow nevertheless to distinguish clearly between the two scenarios with and 
without in-medium effects. Data support the first one. 

Another parameter which has infiuence on the kaon multiplicities is the stiffness of the nuclear 
equation of state. The EOS dependence will in detail be discussed in the next chapter. The usage of a 
stiffer nuclear EOS reduces the yields to some extent, in particular in heavy systems and thus one 
might be worried that the combination of a hard EOS w/o kaon potential might also be able to match 
with existing data. But here one can take advantage from the fact that light systems like C+C show 
almost no nuclear EOS dependence while the influence of the kaon potential is still present. Hence the 
two effects can be disentangled and the scenario without in-medium kaon potential can be ruled out 
from the light reaction systems. 

Concerning the K~^s mesons the situation is more complex and less clear. This is to large part due 
to the interplay of two different productions mechanisms, i.e. strangeness production and strangeness 
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exchange reactions. Here transport model calculations do not yet deliver a consistent picture. On a 
qualitative level the in-medium energy shifts can easily be understood within the mean field picture. In 
K^K~ pair production reactions the vector potentials ±Vo cancel due to alternating signs for and 
K~ and only the attractive scalar parts S5 fj30p lead to a shift AQ of the production thresholds 

BB — > BBK+K- , Ag = -2S5 . (101) 

The same holds for the combination of strangeness production and strangeness exchange which involves 
the excitation of an intermediate hyperon 



BB — ^ BYK+ , Ag = ^0 - S5 
TiY — > BK- , Ag = -Vq - Ss 



Ag = -2Ss . (102) 



The latter can be viewed as a three-body process ttBB — > BBK^K^ which is energetically favored 
compared to the two-body process BB — > BBK'^ K~ . Exactly the same arguments can be applied to 
the pion induced pair production -kB — > BK^K~^ and effective three-body process irirB — > BK^K~ 
which runs over strangeness exchange via an intermediate hyperon. If there are enough pions in the 
system these processes are the dominant sources for K~ production |116|I12()| . From these considerations 
one would conclude that the K~ production threshold is generally lowered by Ag = — 2E5 and the 
corresponding K~ should be significantly enhanced by the presence of the in-medium potentials. 
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Figure 31: rapidity distributions in Ni+Ni reactions at 1.93 AGeV obtained including and without 
in-medium kaon and antikaon potentials. RBUU calculations (Giessen from j204j . Texas from |197j ) 
are compared to data from KaoS |157j . 

However, these arguments are only valid if strangeness production and strangeness exchange reac- 
tions take place at equal nuclear densities where the corresponding potentials are of equal size. This 
is generally not the case (see also Chap. 3). While strangeness production takes predominantly place 
in the early high density phase, strangeness exchange reactions are the driving processes for K~ pro- 
duction and absorption at later stages and lower nuclear densities jllH 1120] . Now the balance between 
attractive and repulsive potential shifts is violated. The same holds when less attractive K~ potentials, 
e.g. from coupled channel calculations are used instead of mean fields potentials. Hence the net effect of 
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the in-medium potentials on the K~ yield may be small |12m 11211 1204j . Also the data situation is much 
less satisfying than for . Some older transport calculations explain the available K~ yields better 
using in-medium potentials |1541 1^ . The comparison to old Ni+Ni in 1997 data clearly needed strong 
K~ potentials. In more recent measurements the K~ yield was found to be about a factor of two lower 
|157j . Fig. [^compares recent RBUU results from Giessen |2U4j and Texas |197j to K~ multiplicities 
measured by KaoS in semi-central Ni+Ni reactions at 1.93 AGeV. Both models support the K~ in- 
medium scenario, but on the basis of a qualitatively different behavior. While the results of Chen at al. 
|197j follow the argumentation above, leading to an enhanced K~ yield when potentials are switched 
on, the results from Mishra et al. |2()4j show the opposite behavior. The K~ yields depend not only 
on the strength of the attractive K~ potential but, due to strangeness exchange, also on the strength 
of the repulsive potential. Around threshold the number of mesons and, correspondingly, that 
of A's in the system is about one order of magnitude larger than that of primary K~ mesons. Small 
relative changes in the A abundances can have large impact on the final K~ yields. This complicated 
interplay can even lead to a reversed potential dependence of the yield as has first been pointed 
out by Hartnack et al. |12Uj . The results from Ref . |2U4j shown in Fig. |^ are based on the chiral 
RHA potential, i.e. the same model as applied in Fig. EOlfor , which is able to reproduce both sets 
of data. Slightly different mean fields which reproduce the data as well were found to fail for K~ 
|2U4j . Hence, the measured K~ yields allow at present no definite conclusions on the strength of the 
attractive potential. The situation is further complicated by a possible strong medium dependence 
of the strangeness exchange reactions tiY < — > NK~ which is, however, theoretically not yet completely 
settled. The predictions obtained within coupled channel calculations range from a moderate enhance- 
ment close to threshold [HHl EZI to a strong suppression |121j . A consistent treatment of these effects 



requires in any case to account for the off-shell dynamics within the transport approach. The off-shell 
calculations of |121j support the scenario of in-medium potential shifts for K~ . 
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Figure 32: K~ / ratio as a function of rapidity in Ru+Ru reactions at 1.69 AGeV and Ni+Ni reactions 
at 1.93 AGeV. RBUU calculations (thick lines: Ref. thin lines: Ref. [206 ) are compared to FOPI 
data |155j . In both cases solid lines include in-medium potentials, dashed lines refer to calculations 
without in-medium potentials. The figure is taken from 



The situation should become clear if one considers the / ratio, in particular its phase space 
dependence. This is done in Fig. 021 where FOPI data |155j for the K~ / ratio as a function of 
rapidity are compared to transport results from |2U6j and [221 • Without further in- medium effects the 
distributions are predicted to be flat as also expected within a statistical approach. The presence of the 
potentials pushes the kaons outwards to higher rapidities while the attractive antikaon potential binds 
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K~^s at mid-rapidity. Both effects lead to an increase of the K~ / ratio around mid-rapidity as also 
seen in the data. Supplemenatry data from KaoS ^157. show that K~ / ratio reaches in Ni-|-Ni at 1.93 
AGeV a value of about 0.04 at mid-rapidity which is in good agreement with the predictions from [22] 
(with pot.) but in contrast to those of |206j were the in-medium effects are over estimated. However, 
more recent calculations which include in-medium modifications of the pion induced K~ productions 
cross sections ttF — > NK~ and the corresponding absorption cross sections do no more deliver such a 
clear picture |68t l62] . The K~ chemistry and the freeze-out time depends crucially on the magnitude of 
the strangeness exchange cross sections and this seems also to be reflected in the corresponding rapidity 
distributions. As already mentioned, the medium modifications of the cross sections are still a matter 
of current debate. 

5.2 Dynamical Observables 

Dynamical observables such as collective flow patterns are to large extent free from uncertainties in the 
total production rates. They depend on the phase space pattern of the primordial sources and the final 
state interaction. For mesons the final state interaction is well under control since only elastic (and 
charge exchange) reactions occur, the total elastic cross section is of the order of 10 -^ 15 mb. 

5.2.1 In-plane flow 




• 


FORI, new 






reflected 




A 


w/o kaon pot. 




T 


with kaon pot 




Ea 


FORI, old 





Y 





(0) 



1.6 



Figure 33: Transverse flow in 1.93 AGeV ^^Ni + ^^Ni reactions at impact parameter b< 4 fm. 
Calculations with and w/o kaon in-medium potential are compared to FOPI data [3j (old) and |2()8j 
(new). The MFT ChPT-|-corr. potential is used brown96b and the Lorentz force is included. 

The transverse or in-plane flow of mesons is a particularly attractive observable. It was first 
proposed by Ko and Li j33 that the kaon flow pattern should provide a sensitive probe for the in- 
medium kaon potential. The repulsive potential should push the kaons away from the nuclear matter 
and produce slight anti-flow at spectator rapidities and a zero flow signal around mid-rapidity. This 
was found to be consistent with the first available flow data from FOPI [S]. Other theoretical studies 
predicted similar features for the kaon flow jHH 11831 1159| lltiOj . However, as pointed out by Fuchs et al. 
|46j the scalar- vector type structure of the kaonic mean field implies the occurrence of a Lorentz- force 
(LF) in moving frames which has been disregarded in the previous investigations. The Lorentz-force 
from the vector field counterbalances the influence of the time-like vector potential on the in-plane 
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flow to large extent (see discussion in Chap. 4) which makes it more difficult to draw definite conclusions 
from transverse fiow pattern. 

In Refs. jl821 l2UUj the fiow in Ni+Ni reactions was re-investigated since in the meantime FOPI 
data with improved statistics became available |2U8j . As can be seen from Fig. ESI it is difficult to 
distinguish between the scenarios w/o in-medium potentials and full co variant in-medium dynamics 
around mid-rapidity. However, at spectator rapidities clear differences appear and the data favor again 
the in-medium scenario. The QMD calculations are based on the MFT ChPT+corr. potential brown96b 
including the LF. Very similar results have been obtained in |2()()j where the same potential (including 
LF) has been used. In both, QMD |182j and BUU |2()()j . the data are best described by the relatively 
strong MFT ChPT+corr. mean field while the weaker MFT ChPT potential or no potential at all lead 
to too strong fiow at spectator rapidities. The dependence of the kaon fiow on the nuclear EOS and the 
Coulomb force is quite weak |159| I182j . The present calculations are based on a soft nuclear EOS and 
the Coulomb force is included. 

Since the magnitude of the Lorentz-force depends on the size of the vector potential an explicit 
momentum dependence of these potentials beyond mean field can reduce the Lorentz force. Such a 
momentum dependence reduces the nucleon fiow in relativistic approaches compared to a mean field 
description, e.g. within QHD. An explicit momentum dependence is necessary in order to comply with 
the empirical optical nucleon- nucleus potential and nucleon fiow data above 1 AGeV |2(J9t 12101 I^Hj - 
Similarly, spectra (piab, Vt and m-r) imply that the KN interaction is less repulsive at high pt which 
might be an indication for an explicit momentum dependence counterbalancing the Lorentz force to 
some extent. Slopes obtained - in particular in central - Au+Au reactions are too hard while C+C 
spectra spectra are well described |2UUj (see also Fig. 1^^ . The FOPI Collaboration measured also the 
Pt dependence of Vi in Ni+Ni and Ru+Ru reactions at spectator rapidities where a transition from 
anti-fiow to fiow with rising px was observed |212j . Also these data require a repulsive in-medium 
potential as found in RBUU |212j . Recent BUU calculations |2UUj indicate that the Lorentz force is 
thereby essential not to overestimate the data. Hence, for a precise determination of the density and 
momentum dependence of the potential more theoretical efforts are needed. 



E 

X 

Q. 

V 



-0.2 



Ni-hNi, 1.93 AGeV 
p/m>0.5 




© FOPI 

▲ all protons 

V kaon producing sources 



0.5 



,(0) 



Figure 34: Transverse flow of the kaon production sources in L93AGeV ^^Ni + ^^Ni reactions, 
total proton flow is compared to FOPI data |^ . 



The 



In this context one might worry about conclusions based on the small kaonic flow signal (about 10~^ 
of the initial kinetic energy) if it is not completely clear how well the flow of the primordial sources is 
under control. The proton flow, shown in Fig. OH is reasonably well described by present transport 
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models but there exist still deviations from data which are generally on the 10% level. The situation 
is, however, less severe since the final proton and kaon fiow pattern are only very loosely connected 
|184j . The production sources carry a large in-plane fiow while the kaons themselves carry a much 
smaller fiow fraction since, at given rapidity, they are produced from two baryons originating from very 
different rapidity regions. Therefore baryon sources with positive and negative px add up to an almost 
vanishing net flow (see the results given by open down triangles in Fig. |^ A's which are produced in 
association with K~^^s have very similar flow pattern as protons, i.e. they show almost the same Px/m 
scaling |2()8j . In |18()j the A flow was investigated within the present model and the data were best 
reproduced including the hyperon mean fleld according to SU(3) scaling f/^^. = ff/J^t- also observed 
in |16Uj the primordial A flow is moderate but strongly enhanced by the AA^ flnal state interactions, 
i.e. rescattering and the A mean fleld. 



5.2.2 Out-oflF-plane flow 
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Figure 35: azimuthal angular distributions in semi-central Au+Au reactions at 1.0 AGeV and Ni+Ni 
reactions at 1.93 AGeV. QMD, BUU |200j and RBUU |213j calculations without and with in-medium 
potential are compared to data from KaoS |ltiH I2l4j . 



The phenomenom of collective flow can generally be characterized in terms of anisotropics of the 
azimuthal emission pattern, expressed in terms of a Fourier series 

dN 

— (0) ocl + 2vi cosU) + 2v2 cos(2</)) + . . . (103) 
d(p 

which allows a transparent interpretation of the coefficients vi and V2. The dipole term vi arises from 
a collective sideward deflection of the particles in the reaction plane and characterizes the transverse 
flow in the reaction plane. The second harmonics describes the emission pattern perpendicular to the 
reaction plane. For negative V2 one has a preferential out-of-plane emission, called squeeze-out. Pions 
exhibit a clear out-of-plane preference j215| I216j which is due to shadowing by spectator nucleons. The 
short mean free path of the pions hinders pions produced in the central reaction zone to traverse the 
spectator matter located in the reaction plane. Therefore it is easier for them to escape perpendicular 



59 



to the reaction plane |136j . Since the K~ mean free path is comparable to that of the pions one might 
expect the same phenomenon for K~ while the mean free path of mesons is large and no squeeze- 
out signal should be observed. These arguments hold when the final state interaction is exclusively 
determined by scattering and absorption processes. However, the presence of a mean field changes the 
dynamics. Azimuthal anisotropics are therefore considered as a suitable tool to study medium effects. 

The azimuthal asymmetry of the production in heavy ion reactions has been first studied by Li 
et al. |217j . First data from KaoS for semi-central Au+Au at 1 AGeV showed a clear squeeze- 
out signal for midrapidity kaons. In corresponding transport calculations from the Texas/Stony Brook 
group [2r7l I161j and the Tiibingen group |179j the data could only be reproduced by the presence of 
the repulsive mean field. Elastic rescattering of mesons was found to be too weak to create 
the observed squeeze-out signal. However, if the repulsive potential is taken into account, the kaons are 
driven by potential gradients preferentially out-of-plane since gradients are larger perpendicular than 
parallel to the reaction plane. With other words, in the reaction plane the kaons are repelled by the 
spectator matter. Thus, the potential leads to an additional dynamical focusing out of the reaction 
plane. 

Fig. IHKl shows the azimuthal distributions for semi-central Au+Au reactions at 1 AGeV and Ni+Ni 
at 1.93 AGeV |214j . In both cases a transverse momentum cut of 0.2 < pt < 0.8 GeV/c and a mid- 
rapidity cut of < 0.2 (Au+Au) and \Y^^^ < 0.4 (Ni+Ni) has been applied. The corresponding 
impact parameters ranges are 5 < 6 < 10 fm (Au+Au) and 3.8 < b < 6.5 fm (Ni+Ni). In the Au+Au 
case we compare QMD calculations and recent BUU calculations from Larionov et al. |200j to the 
KaoS data |161j . The results confirm the findings that the in- medium potential is needed in order to 
explain the experimental squeeze-out signal. Another interesting observation is the fact that the Lorentz 
force, present in covariant dynamics, has only a small influence on the out-of-plane flow, contrary to 
the in-plane flow discussed above |200j ^. The BUU calculations are based on the same kaon mean 
field as QMD, i.e. on Ref. [31]. For the Ni+Ni system shown in the right panel of Fig. EHl also full 
experimental filter cuts are apphed (0.267 < piab < 1.182 GeV/c and 28 < ©lab < 54). The calculations 
with in-medium potential include the full covariant dynamics. For completeness the right panel of Fig. 
EHl shows also recent RBUU/HSD calculations |213j based on the chiral RHA potential |204j which 
agree well with QMD. In summary, there exists a convergence of the various transport models on the 
conclusion that the azimuthal emission pattern require a repulsive mean field. 

Turning to K~ the situation is less clear. The first predictions for K~ out-of-plane emission pattern 
were made by Wang et al. |179j and are shown in Fig. EH The energy of 1.8 AGeV has been 
chosen threshold equivalent to at 1 GeV. The results shown in Fig. 1211 were obtained at mid- 
rapidity (|F(°)| < 0.2) for a semi-central (b=8 fm) reaction. Standard strangeness exchange reactions 
ttY <-i> K~N (Y = A, S) have been taken into account and the MFT ChPT mean field (see Fig. 
has been used. These calculations concentrated on the dynamical mean field effect in the particle 
propagation and thus shifts of the and K~ production threshold by the in-medium potentials have 
been disregarded. 

Under these assumptions a very transparent picture was predicted: due to the short mean free 
path K~ behave similar like pions, i.e. they show a clear squeeze-out signal caused by absorption and 
rescattering. By the presence of a strongly attractive potential this signal is destroyed. In summary: 
without in-medium potential no squeeze-out signal for and a strong squeeze-out signal for K~ , 
with in-medium potentials a sizable squeeze-out signal for and an isotropic emission pattern for 
K~ . However, already in |179j a freeze-out time or correlation has been observed, reflected in the 
Pt dependence of the -Rout/in ratio. -Rout/in quantifies the strength of the azimuthal asymmetry and is 
defined by the ratio of the particle multiplicity emitted perpendicular to that emitted parallel to the 



'I'ln contrast to the earlier statement made in |179j we agree with | 2UU| on this point. 
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Figure 36: K~ azimuthal angular distribution and the -Rout/in ratio as a function of pt as predicted 
in Ref. |179j for semi-central Au+Au reactions at 1.8. The calculations have been performed with 
in-medium K~ potential. A mid-rapidity cut of \Y^'^^ \ < 0.2 has been applied. 

reaction plane 



-Rout/in>l corresponds to a preferrential out-of-plane emission. The calculation without K~ potential 
shows a steady rise of -Rout/in with pt, reflected in the emission pattern (Fig. EHl left part) where a 
Pt > 0.5GeV/c cut has been applied. When the potential is switched on -Rout/in is close to unity but 
high Pt particles which are assumed to freeze out early carry still some squeeze. The low pt K~ mesons 
are, on the other hand, equilibrated and have flat emission pattern or even a slight in-plane flow. Such 
a transition from in-plane to out-of-plane flow has recently been observed by KaoS in Au+Au at 1.5 
AGeV |218j but the experimental signal is much larger than that obtained in |179j . 

First data on K~ azimuthal emission pattern have only recently been delivered by KaoS for the 
Ni+Ni system at 1.93 AGeV |214j . The data agree with none of the two predicted scenarios but a 
preferred in-plane emission of the K~ mesons at mid-rapidity has been observed. Although IQMD 
transport calculations from the Nantes group match the data when the in-medium potential is used it 
is not really obvious how such in-plane flow can develop. In |214j it is argued that it might be due to 
late emission times caused by strangeness exchange reactions K~N — > Yn — > K~ N . However, the 
same holds for pions which undergo several absorption cycles vrA^ — > A — > ttN, have late emission 
times but show a clear squeeze-out signal. The main difference between the 1997 QMD calculations 
from Wang et al. |179j and recent IQMD and RBUU results shown in Fig. EH lies in the fact that 
in ^179j free cross sections without in-medium shifts, both for and K~ production and absorption, 
have been used. Shifts of the thresholds change the and K~ multiplicities which, at a first glance, 
one would not expect to affect flow pattern significantly. However, reduced multiplicities go in line 
with a reduced number of hyperons in the system which reduces the K~ absorption rate and the rate 
of strangeness exchange reactions Yir — > K~N . 

Thus the K~ flow pattern are determined by the interplay between mean fleld and in-medium cross 
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Figure 37: K~ azimuthal angular distribution in Ni+Ni at 1.93 AGeV. RBUU calculations ^213J without 
and using in-medium potentials of different type are shown for 47r- acceptance (left) and compared 
to the KaoS data |214j with the corresponding acceptance cuts (right). 



sections. The medium dependence of the cross sections is significant, in particular when off-shell effects 
are included. The An pattern from recent RBUU calculations |213j demonstrate the following: Without 
any medium effects a clear squeeze-out signal is observed, as predicted in |179j . Using off-shell dynamics 
based coupled channel G-matrix potentials and in-medium cross sections of Tolos et al. jl58^ I121| this 
signal is even enhanced. This potential see Fig. 01 is significantly weaker than the chiral RHA mean field 
but the K~ absorption cross sections are strongly enhanced at threshold. This means that the shadowing 
wins against the attractive potential in this scenario which leads to an even stronger squeeze-out than 
without in-medium effects. When, on the other hand, pion dressing is taken into account (G-matrix 
with TT-dressing) the situation changes and the azimuthal pattern are almost flat. The corresponding 
potential is slightly less attractive than in the previous case but the absorption cross sections are now 
strongly reduced. The last case, namely the strongly attractive chiral RHA potential with free cross 
sections leads to an anti-squeeze-out which is also seen in the data. However, the KaoS acceptance cuts 
lead to a strong distortion of the flow pattern and none of the discussed scenarios is able to explain the 
observed distribution shown in the right panel of Fig. EH Thus the observed K~ emission pattern are 
at present not really understood and further experimental and theoretical efforts are needed to clarify 
the picture. 



5.3 Consistency of the results 

5.3.1 Consistency between transport predictions 

Now the question arises up to which degree a consistent picture has emerged after more than ten 
years of intensive experimental and theoretical efforts to understand kaon production in heavy ion 
reactions at intermediate energies. Since one of the major goals was to extract information on the 
existence and size of in- medium potentials we summarize in Table El the answers which are provided by 
present transport calculations. This summary is organized in the following way: for the most important 
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Table 3: Comparison of various transport model calculations with existing data: the calculations 
are from ^Tiibingen (f besides the given references also results from the present work are included), 
^Nantes, ^Giessen (RBUU/HSD), ^Texas/Stony Brook, ^Giessen(Berkeley). Symbols denote: •= good 
description of available data within error bars; 0= clear failure to describe available data within error 
bars; -k= situation unclear. The first symbol corresponds to calculations based on kaon in-medium 
potentials while the second one corresponds to the case w/o in-medium potential. 

measured observables the agreement or disagreement of transport calculations of different groups, using 
independent simulation codes and models, is classified by three classes: 

A filled bullet (•) denotes a good agreement within error bars with the bulk of existing data for this 
observable, an open bullet (O), on the other hand, denotes a clear disagreement. When the situation 
is unclear, i.e. when calculations match part of the data and fail for other parts, this is indicated by a 
star (^). In each case two symbols are shown. The first one corresponds to the in-medium scenario, i.e. 
the calculations using in-medium potentials, while the second one corresponds to calculations without 
kaon potentials. Hence a combination (•O) means that the corresponding observable allows to clearly 
distinguish between the two scenarios and data support the existence of in-medium potentials, (O*) 
would mean the opposite. 

It is clear that such a classification is rough since the data situation strongly differs from observable to 
observable, the transport models differ at least partially in their input and the level of sophistication and 
do not in each case compare to the complete set of available data in one observable class. Nevertheless, 
there exists the necessity to bring some systematics into the theoretical predictions. 

As already discussed in Chapters 3 and 4 the various transport approaches differ partialy in the 
elementary input and in the technical realizations. Table 0] summarizes the most relevant differences. 
These lie in usage of different paramterizations for the strangeness production cross sections NN ^-^ 
NYK~^^, NA I— >■ NYK^ and the medium dependence of the strangeness exchange cross sections Ntt ^ 
YK~ . The production cross sections are either based on one-boson-exchange (OBE) or resonance 
model (R) calculations where parameters are fixed by the measured pp ^ pAK^ reaction. For the pion 
induced production all approaches apply the cross section of |110j . (P) denotes parameterizations 
of experimental cross sections which are in particular used for the well constrained strangeness exchange 
reactions. However, here the intrinsic medium dependence of these cross sections is crucial and has been 
explored in chiral coupled channel (CC) [HHl and coupled channel G-matrix (CCG) [63j calculations. 
The medium dependence due to shifts of the threshold by the attractive/repulsive mean field are 
usually taken into account (mass shifts). Concerning the dynamics one has to distinguish between a 
non-relativisitc treatment of the mean field (nonrel.) and full covariant dynamics (covariant) which 
includes the vector Lorentz force according to ()79|) and the off-shell dynamics of |28| I121j . 

By far the best measured quantity are multiplicities. High precision data exist for a broad range 
of energies and different mass systems. For the dynamical observables flow and squeeze, data are 
in the meantime also precise enough to constrain the models. For mesons one can summarize the 

■fin QMD different parameterizations are applied for the A and S channel which is indicated in Tab. ^ 
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Table 4: Elementary input, realization of the medium dependence of the cross sections and mean field 
dynamics in the transport models summarized in Table El 
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situation as follows: there exist no data which contradict the in-medium scenario. In contrast, most 
observables can only be described within the in-medium scenario. For K~ the data situation as well as 
the theoretical situation is much less satisfying. No consistent picture has yet emerged which allows to 
discriminate between the two scenarios. 

5.3.2 p + A reactions 

The picture was recently complemented by measurements of the production in proton-nucleus 
reactions |219| |HHI • Although such reactions test only subnormal nuclear densities they are much easier 
to handle than the complicated dynamical evolution of heavy ion reactions. A particularly sensitive 
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Figure 38: Ratios of production cross sections for p+Au/p+C at Tp = 2.3 GeV as a function of the 
kaon momentum, a) Transport calculations including only the Coulomb potential (dash-dotted) and in 
addition a kaon potential of different strength, i.e. 20 MeV (dotted) and 40 MeV (dashed-double-dotted) 
at po- The broken line corresponds to simulations without Coulomb and nuclear kaon potentials. In all 
cases considered here, rescattering in the nucleus has been taken into account, b) The open circles 
are the experimental data. The solid line shows the result of RBUU transport calculations starting 
from the dotted line in the top figure with a baryon potential added. The figure is taken from |219j . 

observable was based on the measurement of the kaon momenta in p+Au and p-|-C reactions. A repulsive 
potential leads to a momentum shift which, in the absence of further final state interactions, results in 
minimal possible momenta 

Pmin = V2mK{Vc{r) + Voir)) . (105) 

This feature is, e.g. well known from the suppression of (3^ emission in heavy nuclei at low positron 
momenta due to the Coulomb shift. Naturally, such a low momentum shift is more pronounced in a 
heavy nucleus. Building now the ratio between the heavy and the light system one is sensitive to density 
effects and thus to the size of the potential. 
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This is demonstrated in Fig. EH where the ratio of the cross section in p+Au/p+C reactions is 
shown as a function of the kaon momentum. Taking the Coulomb potential Vc into account one can 
determine the strength of Vq from the peak in the data. This has been done in |219l ISH] using transport 
calculations and a repulsive kaon potential of Vq ~ 20 ± 5 MeV at po has been extracted. Hence the 
result is consistent with heavy ion reactions and the magnitude of Vq as predicted by effective chiral 
Lagrangiens. Such a potential was also found to be consistent with the measured spectra in p+A 
reactions at subthreshold energies [85] . 

6 Probing the nuclear equation of state 

Heavy ion reactions provide the only possibility to reach nuclear matter densities beyond saturation 
density po — 0.16 fm~'^. Transport calculations indicate that in the intermediate energy range -Eiab — 1 
AGeV nuclear densities between 2 ^ 3po are accessible while the highest baryon densities (~ 8po) will 
probably be reached in the energy range of the future GSI facility FAIR |22Uj between 20-^30 AGeV. At 
even higher incident energies transparency sets in and the matter becomes less baryon rich due to the 
dominance of meson production. Since the knowledge of the nuclear equation-of-state (EOS) at supra- 
normal densities is essential for our understanding of the nuclear forces as well as for astrophysical 
purposes, the determination of the EOS was already one of the primary goals when first relativistic 
heavy ion beams started to operate in the beginning of the SOties |221j . In the following we will briefly 
report the knowledge on the nuclear EOS from a theoretical point of view, then turn to the compression 
phase in heavy ion reactions, give a short review on possible observables and finally discuss the recent 
progress achieved by the kaon measurements. 

6.1 Modeling the nuclear EOS 
6.1.1 Predictions for the nuclear EOS 

Models which make predictions on the nuclear EOS can roughly be divided into three classes: 

1. Phenomenological density functionals: These are models based on effective density depen- 
dent interactions such as Gogny or Skyrme forces |222j or relativistic mean field (RMF) models 
[57j . The number of parameters which are fine tuned to the nuclear chart is usually larger than 
six and less than 15. This type of models allows the most precise description of finite nuclear 
properties. 

2. Effective field theory approaches: Models where the effective interaction is determined within 
the spirit of effective field theory (EFT) become recently more and more popular. Such ap- 
proaches lead to a more systematic expansion of the EOS in powers of density, respectively the 
Fermi momentum kp- They can be based on density functional theory |223t l224j or e.g. on chiral 
perturbation theory |225l 1226^ . The advantage of EFT is the small number of free parameters 
and a correspondingly higher predictive power. However, when high precision fits to finite nuclei 
are intended this is presently only possible by the price of fine tuning through additional param- 
eters. Then functionals based on EFT have approximately the same number of parameters as 
phenomenological density functionals. 

3. Ab initio approaches: Based on high precision free space nucleon-nucleon interactions, the 
nuclear many-body problem is treated microscopically. Predictions for the nuclear EOS are pa- 
rameter free. Examples are variational calculations |:227t 122^ . relativistic |167| 12291 12301 or 
non- relativistic Brueckner calculations and Greens functions Monte-Carlo approaches |232j . 
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Phenomenological models as well as EFT contain parameters which have to be fixed by nuclear prop- 
erties around or below saturation density which makes the extrapolation to supra-normal densities 
somewhat questionable. However, in the EFT case such an extrapolation is safer due to a systematic 
density expansion. One has, nevertheless, to keep in mind that EFT approaches are in general based 
on low density expansions. Many-body calculations, on the other hand, have to rely on the summation 
of relevant diagram classes and are still too involved for systematic applications to finite nuclei. In the 
following we will restrict the discussion mainly to the prediction from many-body calculations. 

In the relativistic Brueckner approach the nucleon inside the nuclear medium is viewed as a dressed 
particle in consequence of its two-body interaction with the surrounding nucleons. The in-medium 
interaction of the nucleons is treated in the ladder approximation of the relativistic Bethe-Salpeter 
(BS) equation 



where T denotes the T-matrix, while V is the bare nucleon-nucleon interaction. The intermediate off- 
shell nucleons in the scattering equation are described by a two-particle propagator iGG. The Pauli 
operator Q accounts for the influence of the medium by the Pauli-principle and projects the intermediate 
scattering states out of the Fermi sea. The Green's function G fulfllls the Dyson equation 



Go denotes the free nucleon propagator while the influence of the surrounding nucleons is expressed 
by the nucleon self-energy S. In Brueckner theory this self-energy is determined by summing up the 
interaction with all the nucleons inside the Fermi sea in Hartree-Fock approximation 



The coupled set of equations ()106p - ()108|) represents a self-consistency problem. 

Fig. inni compares the saturation points of nuclear matter obtained by relativistic Dirac-Brueckner- 
Hartree-Fock (DBHF) calculations using the Bonn potentials |233j as bare NN interactions to non- 
relativistic Brueckner-Hartree-Fock calculations for various NN interactions. The DBHF results are 
taken from Ref. |167j (BM) and more recent calculations based on improved techniques are from 
|229j (Tiibingen). In these calculations the Bonn A interactions matches with the empirical region of 
saturation, however, still at a slightly too high density. For various A^A^ potentials the non-relativistic 
results lie on the so called Coester line which misses the empirical region of saturation. Only by the 
inclusion of 3-body-forces (shown by the dashed square in Fig. the situation can be improved 
|234[ 1235] . The contributions from 3-body-forces (3-BFs) are in total repulsive which makes the EOS 
harder and non-relativistic calculations come close to their relativistic counterparts when 3-BFs are 
included. The same effect is observed in variational calculations |228j . The variational approach shown 
contains relativistic boost corrections to the potential which lead to additional repulsion |228j . Both, 
the BHF calculations from [234' and the variational calculations from f228! are based on the latest AVis 
version of the Argonne potential. In both cases phenomenological 3-body-forces are used, the Tucson- 
Melbourne 3-BF in [234j and the Urbana IX 3-BF (using boost corrections the repulsive contributions of 
the UIX interaction are reduced by about 40% compared to the original ones) in |228j . Except of ChPT 
|236j there exists no systematic generation of 3-BF contributions *. On the other hand, contributions 
from 3-body-forces are to large extent canceled by box diagrams containing resonance excitations and/or 
are partially effectively included in the relativistic approach (see e.g. the discussion in |25nil23Zj)- This 
fact should make an application of DBHF at supra-normal densities more reliable. Fig. HUl compares 

*Next to leading order all 3-BFs cancel while non- vanishing contributions appear at NNLO (236) . 
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Figure 39: Nuclear matter saturation points from relativistic (full symbols) and non-relativistic (open 
symbols) Brueckner-Hartree-Fock calculations based on different nucleon-nucleon forces. The diamonds 
show results from variational calculations. Shaded symbols denote calculations which include 3-body 
forces. The shaded area is the empirical region of saturation. 

the equations of state from the different approaches, i.e. relativistic DBHF from Ref. |229j based on the 
Bonn A interaction^^ |233j , non-relativistic BHF |234j and variational calculations |228j . The latter ones 
are based on the Argonne AVig potential and include 3-body forces. All the approaches use modern 
high precision NN interactions and represent state of the art calculations. The corresponding EOSs 
can also be compared to phenomenological parameterizations. 

In Fig. 1301 two Skyrme EOSs are shown which correspond to the limiting cases of a soft (K=200 
MeV) and a hard (K=380 MeV) EOS. One can conclude from Fig. llHlthat ab initio calculations predict 
throughout a soft EOS in the density range relevant for heavy ion reactions at SIS energies, i.e. up to 
about three times po- There seems to be no way to obtain an EOS as stiff as the hard Skyrme force 
shown in Fig. HUl This observation stands somehow in contrast to the observations made from fits 
to finite nuclei. When density functionals are fine tuned to the nuclear chart, e.g. in RMF theory, 
the corresponding EOS turns out to be relatively stiff [HZ]. The same observation can be made within 
EFT. E.g. the EOS obtained from chiral pion-nucleon dynamics by Finelli et al. |226j is rather soft 
but when phenomenological correction terms are added in order to improve the description of the finite 
nuclei this results in much stiffer EOS. However, one has to keep in mind that finite nuclei constrain 
the interaction at saturation density and below. The predictive power of such density functionals at 
supra-normal densities is therefore restricted. 

6.1.2 Skyrme forces in QMD 

Like for QMD calculations shown in the foregoing chapter the following investigations are based on the 
soft and hard Skyrme parameterizations. These forces are easy to handle, cover the range of uncertainty 
concerning the EOS for isospin symmetric nuclear matter and are therefore widely used in transport 
calculations for heavy ion collisions. 

^The high density behavior of the EOS obtained with different interaction, e.g. Bonn B or C is very similar. |229| 
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Figure 40: The nuclear matter EOS from soft and hard Skyrme forces are compared to the predictions 
from microscopic ab initio calculations, i.e. relativistic DBHF |229j . non-relativistic BHF |234j and 
variational 12281 calculations. 



The QMD N-particle Hamiltonian is given by [23 121 12SH1 



(109) 



The Hamiltonian (|109p contains 2-body interactions which are determined as classical expectation 
values from local Skyrme forces V^^^ supplemented by a phenomenological momentum dependence, an 
effective Coulomb interaction V^^°"' and a Yukawa- type potential V^^^. The Yukawa potential mainly 
serves to improve the surface properties and the stability of the initialized nuclei when used in heavy 
ion collisions. The individual nucleons are described by Gaussian wave packets with fixed width 2\/L. 
This leads to a one-particle Wigner density 



/i(q,k,t) 



.g-(q-q.W)^2/L g-(k-k,(i))2L/2 



(110) 



The Skyrme interaction contains an attractive local two-body part, an effective density dependent 
repulsive two-body part and an nonlocal momentum dependent two-body part. The elementary two- 
body potentials entering into p09|) read then 
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In the case 7 = 2 the density dependent interaction can be derived from a local three-body interaction. 
The parameterizations of 1211 1238j treat 7 as a phenomenological parameter. The folding over the 
Wigner distributions (jllOj) yields the expectation values 
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where pij is an interaction density 



A, = T-^e-^'^'-''^^ (114) 
(47rL) 2 

which arises due to the folding of the two Gaussian wave packets with fixed width 2\/L in coordinate 
space. In analogous way the expectation values of the Yukawa and Coulomb potentials are obtained. 
The parameters 7,5, e in Eq. pi3|) are fitted to the saturation point (po = 0.16 fm~^, Eb = —16 
MeV) and the momentum dependence of the real part of the nucleon-nucleus optical potential. The 
linear density dependence of V^^ is obtained from the point-like 2-body interaction in Eq. ()111|) while 
the nonlinear density dependence is motivated by point-like 3-body interactions"'" With this Hamiltonian 
(jl(J9|) the EOS of isospin saturated nuclear matter, i.e. the binding energy per particle, is of the simple 
forming § 

---f-li-f(^)-T|:^(^y-5''>M^|-|.il-. (n5) 

In contrast to the Skyrme functional (jll5|) where the high density behavior is fixed by the compression 
modulus, in microscopic approaches the compression modulus is only loosely connected to the curvature 
at saturation density. Below 3po, e.g. the DBHF EOS with K=230 MeV is close to the soft Skyrme 
EOS but becomes significantly stiffer at higher densities. 



6.2 Particle production and the compression phase in HICs 
6.2.1 Pions 

With the start of the first relativistic heavy ion programs the hope was that particle production would 
provide a direct experimental access to nuclear EOS |239j . Without additional compression two times 
saturation density should be reached in the participant zone of the reactions where the difference between 
the soft and hard Skyrme EOS is about 13 MeV in binding energy. If the matter could be compressed 
up to 3po the difference is already ~ 55 MeV. It was expected that the compressional energy should be 
released into the creation of new particles, primarily pions, when the matter expands |239j . However, 
pions have large absorption cross sections and they turned out not to be suitable messengers of the 
compression phase. They undergo several absorption cycles through nucleon resonances |136| I137| ITT] 
and freeze out at final stages of the reaction and at low densities. Hence pions loose most of their 
knowledge on the compression phase and are not very sensitive probes for stiffness of the EOS. This 
fact is illustrated in Fig. |^ which shows the tt^ excitation function in minimal bias Au+Au and C+C 
reactions and compares to data from KaoS |2U5|ll42j ^. For C+C there exists practically no dependence 
of the pion yield on EOS, in Au+Au the dependence is moderate, i.e. of the order of ~ 15 — 20%. 

A second and important observation, in particular with respect to the production discussed 
below, is the fact that the enhancement of the pion multiplicities in Au+Au when using a soft compared 
to a hard EOS is almost completely independent of the beam energy. This feature becomes even more 
transparent from the lower part of Fig. |^ which shows the ratio of the total pion multiplicities in 
Au + Au over C + C reactions scaled by the corresponding mass numbers. For this observable the data 
from KaoS |205( I142j are qualitatively well reproduced. The usage of different nuclear forces leads to 
a small shift of the theoretical curves but does not change their slope. The fraction is generally below 

parameter 7 ^ 2 is purely phenomenological. 
^The small Coulomb and Yuka"wa terms as "well as relativistic corrections to the kinetic energy corrections have been 
suppressed. 

^Details on the treatment of pions in the present QMD calculations can be found in |27| EOl 
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Figure 41: Excitation function of the vr^ production cross section in Au+Au (scaled by 10~^) and C+C 
reactions (top), and ratio of the total pion multiplicities in Au+Au over C+C. The calculations are 
performed using a hard/soft nuclear EOS and are compared to data from KaoS |2U5| \1A'J\ . 

unity indicating the larger absorption rate in Au + Au compared to C + C. With rising energy the 
pion suppression in the heavy compared to the light system becomes smaller as also seen in the data. 
The dependence of this observable on the nuclear EOS is rather moderate and would not allow to draw 
some definite conclusions from the model calculations. 

The other observation from Fig. ^2 namely that experimental pion multiplicities are in the SIS 
energy range systematically overestimated in large systems like Au+Au is a common (and well known) 
feature of transport models, independent if they are of BUU or QMD type [133 CSZl 12111 123 12121 • It 
demonstrates that the pion production and/or absorption mechanisms are theoretically not yet fully 
understood. The pion as a Goldstone boson is the lightest meson with a large Compton wave length 
and quantum aspects beyond the semi-classical treatment in transport models could be of particular 
relevance. However, several attempts to go beyond the standard treatments and to account for particular 
quantum effects such as the coupling to AN~^ excitations |243[ 1244} I241j or non-localities, i.e. memory 
effects and time delays in resonance decays |2451 1246j did not really improve on this discrepancy and 
solve the problem. Since pions are a source for kaon production through pion induced reactions one could 
consider this fact as a severe uncertainty. Fortunately, the situation is less worse since the overestimation 
of the experimental pion multiplicities is restricted to low momentum pions which dominate the yields. 
These pions are, however, not energetic enough to contribute to subthreshold kaon production and the 
relevant high momentum tails of the pion spectra are generally well described by the transport models 
[miTTTl lTTj. 

However, for a fair comparison of experiment and theory it should be mentioned that also the 
experimental situation, in particular the Apart dependence of the multiplicities, is not yet completely 
settled for the SIS/BEVALAC energy range. FigJ^ shows a compilation of published pion multiplicities 
(vr" + 7r° + 7r+) for different mass systems as a function of Apart !134j at 1 AGeV incident energy. 
The data set consists of measurements with the streamer chamber at LBL |247j and with the TAPS 
[21H1 1213 l2Sni, KaoS |2ni 1203 and FOPI j2S21 spectrometers at the GSI. In none of the experiments 
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Figure 42: Pion multiplicities as a function of Apart for various mass systems at 1 AGeV incident energy. 
QMD calculations for central reactions (b=0 fm) are compared to a compilation of experimental data 
taken form |134j . Recent preliminary results from FOPI |24Uj are denoted as FOPI new. 



all pion charges have been measured simultaneously and thus each data set contains extrapolations 
in isospin space which are based on the isotopic relations for the A production and decay |252j . The 
largest experimental uncertainty lies probably in the determination of Apart- On the one side this 
quantity is biased by the impact parameter resolution but it depends on model assumptions as well. " 
The corresponding QMD calculations are performed for central reactions (b=0 fm) with Apart = A + A. 
The experimental yields are well described in light and intermediate mass systems but are significantly 
overestimated in heavy systems with Apart > 120. The La+La data point |247j . which falls off the 
general systematics, is compatible with the transport result within error bars. Interesting is, however, 
that the FOPI Collaboration re- measured tt^ yields and a new analysis, denoted in Fig 021 as FOPI new, 
is compatible with the streamer chamber systematics and also with the transport model predictions. 

6.2.2 Kaons - historical overview 

After pions turned out to fail as suitable messengers, mesons were suggested as promising tools 
to probe the nuclear EOS. This idea was first put forward by Aichelin and Ko almost 20 years ago 
|132j . At subthreshold energies mesons are produced in the high density phase and due to the 
absence of absorption reactions they have a long mean free path and leave the matter undistorted by 
strong final state interactions. Moreover, at subthreshold energies nucleons have to accumulate energy 
by multiple scattering processes in order to overcome the threshold for kaon production and therefore 
these processes should be particularly sensitive to collective effects. 

Already in the first theoretical investigations by transport models it was noticed that the yield 
reacts rather sensitive on the EOS |176| 1^ 1^ I133j . Both, in non-relativistic QMD calculations based 
on soft/hard Skyrme forces |^ 1^ and in RBUU |176| 113^1 1^3^ with soft/hard versions of the (non- 
linear) crcu-model for the nuclear mean field it turned out that the yield is about a factor 2-3 

" For the determination of the overlapping volume of the two interpenetrating nuclei as function of the impact parameter. 
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Figure 43: QMD and RBUU transport model calculations from Nantes, Texas and Giessen are compared 
to the first KaoS data [2] for Au + Au. The figure is taken from 

larger when a soft EOS is applied compared to a hard EOS. At that time the available data favored 
a soft equation of state. This fact is illustrated in Fig. |331 which compares QMD calculations from 
the Nantes [12] and RBUU calculations from the Texas/Stony Brook |133j and Giessen jHSj groups 
to the first KaoS data for Au + Au j2] • However, at that stage the theoretical calculations were still 
burdened with large uncertainties. First of all, it was noticed [201122] that the influence of the repulsive 
momentum dependent part of the nuclear interaction, Eq. leads to a strong suppression of the 

kaon abundances which made a quantitative description of the available data more difficult. Moreover, 
at that time the pion induced reaction channels ttB — > YK^ have not yet been taken into account 
explicitely. In |176j the frozen A approximation has been used which includes pionic degrees of freedom 
implicitly. The importance of this channel was first pointed out by Fuchs et al. [109 . These additional 
channels which contribute up to 30-i-50% to the total yield allowed to explain the measured yields with 
realistic momentum dependent interactions |1U91 IM] . However, the dependence of the total yield 
on the nuclear EOS turned now out to be much smaller than originally expected, i.e. in the order of 
15-20%. 

A further breakthrough was achieved when the COSY-11 Collaboration measured the pp — pK~^A 
reactions at threshold |HZ1 • The strangeness production cross sections A^A^ — > NK^Y [HZl I102j which 
are nowadays in use are based on these data and are in particular close to threshold three orders of 
magnitude smaller than the parameterizations of Randrup and Ko jHEj which were used in the early 
QMD and RBUU/Texas calculations shown in Fig. 

6.3 The ratio Au+Au/C+C 

Within the last decade the KaoS Collaboration has performed systematic measurements of the 
production far below threshold [21 114H I2U5| 11401 1142j . Based on the new data situation, in Ref . |181j 
the question if valuable information on the nuclear EOS can be extracted has been revisited and it 
has been shown that subthreshold production provides indeed a suitable and reliable tool for 
this purpose. These results have been confirmed by the Nantes group later on ||96j. In subsequent 
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Figure 44: EOS dependence of the multiplicities as a function of Apart in Au+Au, Ni+Ni, C+Au 
and C+C reactions at 1 AGeV. For Au+Au and Ni+Ni data from KaoS |134[ I141j are shown as well. 

publications these findings were worked out in more detail |2531 1254j . Here we summarize and complete 
these investigations. If not denoted differently, throughout this section all model calculations contain a 
repulsive in-medium potential as discussed in Chap. 5. 

In Chap. 3 we discussed already the Apart dependence of the yield on a qualitative bases. The 
calculations shown in Fig. IT3l were based on a soft EOS. Fig. demonstrates the interplay between 
Apart, system size and the nuclear EOS. It shows the multiplicities as a function Apart in Au+Au, 
Ni+Ni, C+Au and C+C reactions at 1 AGeV. The multiplicities are normalized to the mean Apart: 
(Apart) = Amax/2 with A^ax = A + A for Symmetric systems and Amax = 56 for C+Au. A significant 
dependence of the kaon multiplicities on the nuclear EOS requires a large amount of collectivity which is 
easiest reached in central reactions of heavy mass systems. Consequently, the EOS dependence is most 
pronounced in central Au+Au reactions. Also in Ni+Ni effects are still sizable while the small C+C 
system is completely insensitive on the nuclear EOS even in most central reactions. The available data 
for Au+Au and Ni+Ni jl34[ I141j support the soft EOS. Interesting is in this context the asymmetric 
C+Au system: Though in central C+Au reactions the number of participants is comparable to Ni+Ni 
the yield does not depend on the EOS. This indicates again that a sensitivity on the EOS is not only 
a question of Apart but of the compression which can be reached by the colliding system. Remarkable 
is the saturation of the kaon yield as a function of Apart predicted by the transport calculations. It 
stands in clear contradiction to the Apart dependence predicted by the thermal model [ISUI I151j and a 
measurement of the quantity would allow to distinguish between these two approaches. In Fig. E31the 
QMD calculations for inclusive reactions are compared to KaoS data from |255] . The data are generally 
described within error bars by the soft EOS. Compared to Fig. 031 in inclusive reactions the EOS 
dependence survives for large mass systems although the large difference between the soft and hard - 
seen in most central Au+Au reactions - is washed out to some extent. In minimal bias reactions the 
bulk of kaons originates from semi-central reactions 6 ~ 5 fm, corresponding to Apart /Amax ~ 0.7. The 
fact that the multiplicities in C+Au are significantly smaller than in Ni+Ni reactions although both 
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Figure 45: multiplicities in inclusive C+C, Ni+Ni, Au+Au and C+Au reactions at 1 AGeV. QMD 
calculations using a hard/soft nuclear EOS are compared to KaoS data |255j . The figure is taken from 



systems correspond to comparable mean Apart is a strong experimental evidence for a strong sensitivity 
of the kaon production on the compression. As also discussed in |255j even in central reactions the 
carbon projectile is simply to small in order to achieve a high compression of the gold target. The 
situation is different in central Ni+Ni reactions. Hence the mass dependence of the kaon multiplicities 
seems not primordially to be an Apart effect but a compression effect. This is also clearly reflected by 
the transport calculations. The next and natural step is to consider the energy dependence of the EOS 
effect. It is expected to be most pronounced most far below threshold because there the highest degree 
of collectivity, reflected in multi-step collisions, is necessary to overcome the production thresholds 
(see also discussion in chapter 3.2). The calculations for the excitation function shown in Fig. EHlwere 
obtained for a soft EOS and are performed under minimal bias conditions with femax = 11 fm for Au + Au 
and 6rnax = 5 fm for C + C and normalized to the experimental reaction cross sections jl42[ I205j . For 
both systems the agreement with the KaoS data is quite good. 

The effects become even more evident when the ratio R of the kaon multiplicities obtained in Au+Au 
over C+C reactions (normalized to the corresponding mass numbers) is built |18HI142] . Such a ratio has 
moreover the advantage that possible uncertainties which might still exist in the theoretical calculations 
should cancel out to large extent. 

This ratio is shown in Fig. 0^1 Both, soft and hard EOS, show an increase of R with decreasing 
energy down to 1.0 AGeV. However, this increase is much less pronounced when the stiff EOS is 
employed. In the latter case R even decreases at 0.8 AGeV whereas the soft EOS leads to an unrelieved 
increase oi R. At 1.5 AGeV which is already very close to threshold the differences between the two 
models become small. The strong increase of R can be directly related to higher compressible nuclear 
matter. The comparison to the experimental data from KaoS |142j . where the increase of R is even 
more pronounced, strongly favors a soft equation of state. 

6.3.1 Phase space dependence 

To obtain a quantitative picture of the explored density effects in Fig. |^the baryon densities are shown 
at which the kaons are created. The energy is chosen most below threshold, i.e. at 0.8 AGeV and only 
central collisions are considered where the effects are maximal. dMx+/dp is defined as in Eq. (j50|) . For 
the comparison of the two systems the curves are normalized to the corresponding mass numbers. 
Fig. im illustrates several features: Only in the case of a soft EOS the mean densities at which kaons 
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Figure 46: Excitation function of the ratio R of multiplicities obtained in inclusive Au+Au over 
C+C reactions. The calculations are performed with a hard/soft nuclear EOS and compared to the 
data from the KaoS Collaboration |142j . 

are created differ significantly for the two different reaction systems, i.e. < p/po >=1.46/1.40 for C+C 
and 1.47/1.57 for Au+Au using the hard/soft EOS. Generally, in C+C reactions densities above 2po 
are rarely reached whereas in Au+Au the kaons are created at densities up to three times saturation 
density. Furthermore, for C+C the density distributions are weakly dependent on the nuclear EOS. The 
situation changes completely in Au+Au. Here the density profile shows a pronounced EOS dependence 
|133j . Moreover, the excess of kaons obtained with the soft EOS originates almost exclusively from high 
density matter. A second quantitative measure for the collectivity is the average number of collisions 
for those hadrons which were involved in the production displayed in Fig. ^\ Again only central 
collisions are considered where the effects are maximal. < Nq > is defined as in Fig. (jll|) . In average the 
particles undergo about twice as much relevant collisions in the heavy compared to the light system. 
Furthermore, the collectivity, i.e. the accumulation of energy by multiple scattering, increases with 
decreasing energy. Thus one can conclude that the increase of R is not due to a trivial phase space 
effect, namely the fact that far below threshold the C+C system is simply too small to provide enough 
collectivity for the kaon production. If such a scenario - which in principle also explain the rise of R 
seen in the KaoS data - would be true, < Nc > would have to saturate for C+C at low energies. This 
demonstrates that production far below threshold always requires a certain amount of collectivity 
which can be provided also in a very small colliding system, though such processes are rare. There is, 
however, no sharp limit were such collision histories become impossible. Thus trivial phase space effects 
can be excluded for an explanation of the increase of R. In [142] a similar argument was based on the 
measurement of high energy pions which can test the phase space available for particle production. 

6.3.2 Stability of the EOS dependence 

Now the question arises, how firm the conclusions on the nuclear EOS are. The influence of the repulsive 
in-medium potential has been discussed in \181\ and can be seen from Fig. 1^ which shows the ratio R 
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Figure 47: Phase space dependence of the production. Left: nuclear density at the production of 
mesons. Right: number of coUisions which the particles encountered prior to production. In 
both cases the EOS dependence in central Au+Au and C+C reactions is studied. 

in central reactions from simulations with and without kaon potential. It is remarkable and at a first 
glance surprising that R shows for both cases qualitatively the same behavior. The in-medium kaon 
potential acts opposite to the EOS effect: a higher compression increases the kaon yield but also the 
value of the in-medium kaon mass which, on the other hand, tends to lower the yield again. However, the 
increase of the in-medium mass goes linear with density whereas the collision rate per volume increases 
approximately with p^. E.g. in central Au+Au reactions at 0.8 AGeV the average density < p > at 
kaon production is enhanced from 1.47 to 1.57 po switching from the hard to the soft EOS. This leads 
to an average shift of the in- medium mass (j91|) compared to the vacuum value of 55/61 MeV using the 
hard/soft EOS, i.e. a relative shift of 6 MeV between soft and hard. However, collective effects like the 
accumulation of energy by multiple scattering show a higher sensitivity on the compression resulting in 
an enhancement of the available energy < -^i >= 90 MeV applying the soft EOS. 

The next question concerns the knowledge of the elementary reaction cross sections. As discussed 
in Chapter 4 the NN and ttN cross sections are well under control since these channels are constrained 
by data. The reactions which involve nucleon resonances in the initial states [i = NA, vrA, AA) are 
less secure due to the lack of data and one has to rely on model assumptions. The cross sections which 
have been used in the present transport calculations are based on the effective Lagrangian model of 
Refs. |102[ I110[ I207j . The isospin dependence of the cross sections was determined by isotopic relations 
assuming isospin independent matrix elements. 

In summary, some uncertainty in the transport calculations is still existing due to the fact that 
elementary production channels involving A resonances are not constrained by data ** . The ratio built 
from different mass systems should, however, be robust against such uncertainties: 



Changes of the production cross sections shift absolute yields but considering the ratio possible 



**Elcmcntary reactions with a N*{1AA0) resonance in the initial states are included in the present calculations. For 
these reactions the same cross sections as for nucleons are used. Higher lying nucleon resonances can be neglected at 
subthreshold energies. 
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Figure 48: Excitation function of the ratio of multiplicities obtained in central (b=0 fm) Au+Au 
over C+C reactions. The calculations are performed with/without in-medium kaon potential and using 
a hard/soft nuclear EOS. The figure is taken from |181j . 

errors drop out in leading order. 

• Conclusions are based on the slope of this ratio as a function of energy. It is extremely unlikely 
that an incomplete knowledge of cross sections, i.e. an unknown isospin dependence, can create 
the observed energy dependence. The systematics of spurious contributions should rather be flat 
as a function of energy. 

To illustrate these arguments the influence of the different elementary channels is shown in Fig. 
EHl There the ratios Ri are built separately for the production channels with initial states i = 
NN,7rN, NA,7rA, AA. The shape of R is not strongly influenced by the A^A , vrA channels which 
are the most insecure ones. The excitation function for the A^A contribution varies only little as a func- 
tion of energy and is similar using the different EOSs. The contribution of the vrA channel is decreasing 
for both, a hard and a soft EOS. The shape of R is to most extent determined by the A^A^ and vrA^ 
contributions which are well under control. 

These findings are generally confirmed by independent transport calculations of the Nantes group 
using the IQMD transport model [SH] shown in Fig. EU] together with the present QMD results. The 
IQMD calculations include an in-medium kaon potential derived in relativistic mean field theory (RMF) 
[HB] which is somewhat less repulsive than that one used in the present calculations. For the soft EOS the 
IQMD calculations coincide almost with the present results |181j . The two sets of transport calculations 
show a good overall agreement and both rule out the hard EOS from the comparison with data. The 
shaded area in Fig. EOl can be taken as the existing range of uncertainty in the theoretical model 
description of the considered observable. 

Moreover, the IQMD calculations were also repeated with an alternative set of A^A; A A NYK^ 
cross sections taken from |2S1 which are almost one order of magnitude smaller than those from Tsushima 
et al. [102J (see Fig. 1221) • The ratio R is almost completely independent on this change in elementary 
cross sections and also total yields. An even more systematic study of possible uncertainties concern- 
ing this observable has been performed by Hartnack in ^23Bj . A possible medium dependence of the 
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Figure 49: Dependence of the excitation function of R on the various K'^ production channels. Central 
(b=0 fm) Au+Au and C+C reactions are considered. The results are taken from |253j . 

elementary cross sections, i.e. a density dependent reduction, as well as a general scaling has been 
investigated. Also in these extreme cases the EOS dependence of the ratio R survived and conclusions 
stayed stable. This demonstrates once more the robustness of this observable. 

6.4 Constraints from other sources 
6.4.1 Nucleon flow 

Concerning the nuclear equation of state one has to confront the information from subthreshold pro- 
duction with the knowledge obtained from other sources: At intermediate energies heavy ion reactions 
test the density range between one and two, maximally three times nuclear density. The information 
from kaon production implies that in this density range the EOS shows a soft behavior. One has of 
course to be aware that the adopted Skyrme forces are simplified interactions which are easy to handle 
but must not be very realistic. 

Another observable which helps to constrain the nuclear mean field and the underlying EOS at 
supra- normal densities is the collective nucleon flow |258j . The transverse flow vi has been found to be 
sensitive to the EOS and, in particular in peripheral reactions, to the momentum dependence of the 
mean field |198t I209ll259j . The elliptic flow f2, in addition, is very sensitive to the maximal compression 
reached in the early phase of a heavy ion reaction. The cross over from preferential in-plane flow f 2 < 
to preferential out-off-plane flow f 2 > around 4-6 AGeV has also led to speculations about a phase 
transition in this energy region which goes along with a softening of the EOS ^260j . However, the 
corresponding AGS data can also be explained conventionally 261 . 

The present situation can be summarized as follows: Flow data at SIS energies are consistent with 
a soft EOS I2IIII CnH]. The full flow excitation function, ranging from low SIS {E^^^, ~ 0.2 ^ 2 AGeV) 
up to top AGS energies (-Eiab — 2 11 AGeV), has been studied in |257t 1261] . The conclusion from 
Ref. j257j was that, both, super-soft equations of state (K=167 MeV) as well as hard EOSs (K>300 
MeV) are ruled out by data. Fig. ED displays the pressure- density area which, according to the analysis 
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Figure 50: Excitation function of the ratio R of multiplicities obtained in inclusive Au+Au over 
C+C reactions. Our results are compared to IQMD calculations [96; . The shaded area indicates thereby 
the range of uncertainty in the theoretical models. In addition IQMD results based on an alternative 
set of elementary K'^ production cross sections are shown. 



of |257j . is consistent with heavy ion flow data. The soft Skyrme EOS is in agreement with flow data. 
The boundaries of Fig. |^ are the result of a compilation from the analysis of sideward and elliptic 
anisotropies. In the models used by Danielewicz et al. |198| I257j sideward flow favors indeed a rather 
soft EOS with K=210 MeV while the development of the elliptic flow requires slightly higher pressures. 
The BHF and variational calculations including 3-body-forces^^ fit well into the constrained area up to 
4po- At higher densities the microscopic EOSs, also DBHF, tend to be too repulsive. 

However, conclusions from flow data are generally complicated by the interplay of the compressional 
part of the nuclear EOS and the momentum dependence of the nuclear forces. A detailed comparison to 
Vi and V2 data below 1 AGeV from FOPI \2b'2\ and KaoS |263j favors again a relatively soft EOS with a 
momentum dependence close to that obtained from microscopic DBHF calculations |198t I209| I264j . In 
Fig. I^the microscopic DBHF EOS (K=230 MeV) lies at the upper edge of the boundary, but is still 
consistent in the density range tested at SIS energies, i.e. up to maximally 3 po- This fact is further 
consistent with the findings of Gaitanos et al. |2091 1264j where a good description of Vi and f 2 data at 
energies between 0.2 and 0.8 AGeV has been found in RBUU calculations based on DBHF mean fields. 
However, as pointed out in Refs. |209t ITHl 12641 1211j is thereby essential to account for non-equilibrium 
effects and the momentum dependence of the forces which softens the EOS compared to the equilibrium 
case shown in Fig. |^ 

In summary, production and nucleon flow provide a consistent picture so far that stiff equations- 
of-state are ruled out. Details concerning the interplay between density and momentum dependence 
have still to be settled and require future efforts. 



^^^^For the BHF + 3-BF calculation the pressure shown m Fig. [^has been determined from the parameterization given 
in |235| which is based on the Urbana IX 3-BF different to that used in |234| . 
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Figure 51: Constraints on the nuclear EOS from heavy ion flow data. The shaded area shows the 
pressure- density which is compatible with heavy ion flow data according the analysis on |257j . The 
equations-of-state from the models shown in Fig. HUl are displayed. 

6.4.2 Neutron stars and symmetry energy 

Models for neutron stars are constrained by the lower limit of the maximal neutron star mass. This 
means that any equation-of-state must be stiff enough to produce a neutron star of mass greater than 
1.44 solar masses, the largest neutron star in the PSR 1913+16 system. The upper limit for a nuclear 
EOS is thereby obtained by the conventional neutron star consisting of neutrons and protons. The 
occurrence of additional degrees of freedom such as pion, K~ or H-dibaryon condensation or the exci- 
tation of hyperons in hybrid stars softens the EOS and reduces the maximal neutron star mass. The 
same is true when quark cores or cores of strange matter are considered. Conventional neutron stars 
put only weak constraints on the nuclear EOS. Super-soft EOSs with a compressibility less than K< 
120 MeV can be ruled out The microscopic models discussed above yield maximal neutron star 
mass above two solar mass which are all very close: DBHF (Bonn A) [231J gives Mmax = 2.26M0, 
BHF + 3-BF Mmax = 2.3Mq 1235" and the variational calculations with 3-BFs + boost correction gives 
Mmax = 2.2IM0 |228j . A soft EOS comparable to the soft Skyrme force, i.e. a chiral relativistic mean 
field model with K=194 MeV used in jEj yields Mmax = 2.OM0. 

Since strangeness is not conserved in weak interactions K~ condensation can occur in neutron stars 
at densities above 3 -r- 5po [121 UHl UH UHl ISS]; depending on the strength of the K~N interaction and 
the nuclear EOS. The K~ condensate introduces additional negative charge which enhances the proton 
fraction in the star and makes the EOS softer. Its influence depends therefore strongly on the high 
density behavior of the symmetry energy. It has been found that the maximal mass is reduced by about 
20 -7- 25% by a K~ condensate E] which would then rule out an EOS softer than K< 180 MeV. 
Neutron stars constrain the isospin symmetric EOS from below. 

Of particular interest is in this context the symmetry energy which characterizes the isospin de- 
pendence of the nuclear forces and determines the proton fraction inside a neutron star. In isospin 
asymmetric matter the binding energy is a functional of the proton and neutron densities. Asymmetric 
matter is characterized by the asymmetry parameter f3 = Yn — Yp which is the difference of the neu- 
tron and proton fraction Yi = pi/ p ,i = n,p. The isospin dependence of the energy functional can be 
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Figure 52: Left: symmetry energy of the various models as a function of density. Right: K'^ jK^ ratio 
in nuclear matter at density 2.5 po and temperature T = 60 MeV as a function of the proton-neutron 
asymmetry /3 obtained from RBUU calculations with periodic boundary conditions. The results are 
taken from |199j . 



expanded in terms of /3 

E{p, (5) = E{p) + Esym(p)/5' + 0{f3') + ■■■ . (116) 

The high density behavior of the symmetry energy is at present largely unconstrained and kaons could 
again turn out to be a suitable tool to derive constraints from heavy ion collisions. In this case the 
/ ratio has been suggested as a promising observable by the Catania group |199j . The isospin 
dependence of the elementary production cross sections reflects the asymmetry of the emitting source 
which is connected to the stiffness of symmetry energy and thus translated into the corresponding 
and yields. RBUU calculations with periodic boundary conditions for equilibrated nuclear matter 
at 2.5 times saturation density and a temperature of T = 60 MeV, shwon in Fig. indicate a strong 
sensitivity on -Egym and the isospin dependence of the nuclear forces. The calculations are based on non- 
linear Walecka mean field including either the vector isovetor p-meson (NLp) or in addition the scalar 
isovector 5- meson {NLp6). The NL model contains only isoscalar a and u mesons and has therefore 
Egym = 0. The measurement of K~^'^ mesons in symmetric and highly isospin asymmetric colliding 
systems might therefore provide experimental access to the isospin dependent part of the nuclear EOS. 



7 Summary and outlook 

Kaon and antikaon production in nucleus-nucleus collisions around threshold energies opens the pos- 
sibility to attack a variety of physics questions which have important implications in nuclear physics, 
QCD and astrophysics. The present article tried to summarize the status of the field and to point out 
which problems could be settled and which questions are still open. 
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Theory predicts strong modifications of tlie kaon and antikaon properties in a dense liadronic envi- 
ronment. Mean field models as well as chiral perturbation theory predict a repulsive potential of 
about Vk+ — +(20 -r- 30) MeV at nuclear saturation density. Such a value is in agreement with em- 
pirical kaon-nucleon scattering. The i^~-nucleon interaction, in contrast, is resonant around threshold 
and requires non-perturbative approaches. The strength of the K~ potential is still an open question. 
The depth of the attractive antikaon-nucleon potential ranges from Vk- — —(50 -h 100) MeV, obtained 
within chiral coupled channel dynamics, to Vk- — — (lOO-i-200) MeV predicted by mean field approaches 
and the analysis of kaonic atoms. In contrast to the mesons the K~ mesons develop complicated 
spectral properties in the medium. 

In heavy ion collisions at intermediate energies, i.e. at energies around the threshold region, 
strangeness is generally produced in the early and high density phase of the reaction. However, the 
freeze-out conditions for kaons and antikaons are completely different. Due to strangeness conservation 
mesons cannot be reabsorbed by the surrounding nucleons and their chemical freeze-out takes place 
early. Final state interactions, i.e. elastic scattering or charge exchange reactions and the influence 
of the optical kaon-nucleon potential change their dynamical pattern but not the abundances. This 
makes mesons to a suitable 'penetrating' probe to study the dense fireball created in a heavy ion 
reaction. Antikaons, in contrast, are strongly coupled to the environment through strangeness exchange 
reactions. This leads to a late freeze-out and a loss of memory on the early reaction stages. 

The link of the underlying physics to the heavy ion experiments must be provided by dynamical 
transport models. The precision of such models depends thereby crucially on the elementary input, 
i.e. the knowledge of the elementary reaction cross sections. This input is much better constrained by 
data for the mesons than for K~ . In addition, the quasiparticle picture which underlies all types of 
semiclassical transport approaches, is much better justified for the kaons than for the antikaons. Hence 
state of the art transport calculations have reached a reasonable degree of consistency concerning 
production and dynamics. The comparison to experiment concerning both, total yields as well as 
dynamical observables such as in-plane and out-off-plane flow pattern, supports the existence of a 
slightly repulsive in-medium potential as predicted by chiral dynamics. This picture is complemented 
by data from proton-nucleus reactions. 

Concerning the antikaons the situation is less satisfactory, both from the experimental as well as 
from the theoretical side. To settle the question of a strongly attractive /^"-nucleon potential is one of 
the major challenges of this field in future. If the i^~-nucleon potential is strong enough, this can lead 
to deeply bound K~ states in nuclei |265j and even to light kaonic nuclear clusters . Such states 
would be much stronger bound than the 7r~ states observed in pionic atoms [51j. Some calculations 
predict even a collapse of the nuclear wave functions to densities significantly above saturation density 
in if~-nuclear clusters which would allow to access the nuclear forces at very short distances. The 
existence of such molecule states would thus open a completely new field in hadron physics with strong 
implications on nuclear structure. Heavy ion reactions as well as proton-nucleus reactions are the tools 
to clarify the preconditions for this hypothesis. 

To draw firm conclusions on the in-medium antikaon properties will, however, require significant 
efforts to control their off-shell dynamics. The same holds for other hadrons, e.g. vector mesons, 
where dramatic changes of their spectral properties are expected. First attempts towards a quantum 
transport have been made but an exact treatment of the quantum evolution equations requires a better 
knowledge of off-shell transition elements and in-medium spectral functions of all the involved hadron 
species. To develop a consistent quantum transport is the major challenge of future theoretical heavy 
ion physics. The same holds for proton-nucleus reactions. Also there the quantitative understanding of 
high precision meson production data, e.g. from GSI, COSY or KEK, requires to control the off-shell 
dynamics of such processes. 

Finally the kaons turned out to provide a suitable tool to attack another longstanding question, 
namely the stiffness of the nuclear equation-of-state. The high density behavior of the EOS has severe 
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astrophysical consequences since it determines e.g. the maximal mass and the radii of neutron stars. 
Although heavy ion reactions test mainly isospin symmetric matter they put constraints on theoretical 
models which are also applied to neutron stars. The systematic measurement of the excitation 
function in heavy and light systems down to energies far below threshold can here be considered as a 
breakthrough. Extreme subthreshold energies exclude distortions from surface effects and ensure that 
the mesons originate from supra-normal nuclear densities. To overcome the production thresholds 
requires a high degree of collectivity which, on the other hand, introduces the sensitivity on the com- 
pression achieved in the reaction. The comparison with data strongly supports an EOS which is 'soft' 
in a density regime between 1-3 times saturation density. Such a behavior is consistent with the pre- 
dictions from microscopic many-body calculations and the constraints obtained from nucleon flow data 
in heavy ion reactions. Based on similar arguments the / ratio has recently been proposed as a 
tool to access the isospin dependence of the nuclear EOS and to constrain the symmetry energy. If this 
turns out to be true, subthreshold kaon production can probably be considered as the most successful 
observable to constrain nuclear forces at high densities. 
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